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Section 1 


SUMMARY * 

This report provides user instructions and software descriptions 
for the Base Program of the Coupled Rotor/Airframe Vibration 
Analysis of Reference 1. The Base Program, SIMVIB, was developed 
to simulate the dynamics of a coupled helicopter rotor and air- 
frame. The program can be used as a design tool for predicting 
helicopter vibrations. Parametric studies can be performed during 
the preliminary design phase of an aircraft. The analysis is 
capable of representing various dynamic configurations, including 
fixed system and rotating system vibration absorbers and vibration 
isolation devices . 

Information is provided on the functional capabilities and pro- 
cedures for running the SIMVIB program. Interfaces with external 
programs are discussed. The method of synthesizing a dynamic 
system and the various solution and operational modes are de- 
scribed. A detailed description of the program input data and 
output results is presented. 

The section "Overall Program Structure" provides detailed informa- 
tion on the overlay segmentation structure, flow diagrams, Fortran 
subroutines descriptions. Common blocks, data transfer files, 
storage size of the computer program, and necessary additions and 
revisions for incorporating a new component. 

Sample test case results are provided and discussed. The complete 
input data and output results are included in this manual for five 
representative cases . System responses are plotted to demonstrate 
the plot capability of the SIMVIB computer program. The Job 
Control Language (JCL) required to execute the computer program on 
the IBM 370/168 is included. Program installation instructions 
for the CDC system complete this user's manual. 

Sample test case listings are provided for key test cases in their 
entirety in Volume II of the Program Manual. 


*The research effort which led to the results of this report 
was financially supported by the Structures Laboratory, 
USARTL, (AVRADCOM). 




Section 2 


INTRODUCTION 

The purpose of computer program SIMVIB is to simulate the coupled 
dynamics of a helicopter rotor and airframe. The program can be 
used as a design tool for predicting helicopter vibrations . The 
analysis can predict vibrations for parametric studies conducted 
in preliminary design- User features provided in the program 
include minimum core storage, cost effective operation, and 
interactive input/output options. The SIMVIB program has a 
flexible capability for representing various dynamic configura- 
tions, including fixed system and rotating system vibration 
absorbers and vibration isolation devices. It is an advanced 
design tool for coupled helicopter rotor/airframe vibrations, 
which overcomes the deficiencies of other current programs. The 
basis of the analysis is described in Reference 1. 

This manual provides user's and programmer's instructions. The 
instructions enable a user to prepare input to synthesize config- 
urations for vibration analysis. Sample outputs are interpreted. 
Program structure and subroutine functions are defined to aid the 
development of modifications. 

Section 3 

COMPUTER PROGRAM FUNCTIONAL CAPABILITIES 


The Coupled Rotor/Fuselage Vibration Analysis is a package of 
programs consisting of a Base Program, called SIMVIB, communicat- 
ing with a set of External Programs to achieve a vibration model- 
ing capability for helicopters. A substructure assembly method is 
used in the Base Program to assemble a dynamic model of the 
helicopter from physical helicopter components contained in the 
Base Progreim and External Program models of the rotor system. The 
Base Program solution modes enable the user to obtain steady state 
forced responses, eigensolutions , time history solutions, and 
parameter variations yielding the effects of variations in system 
properties. After the rotor data are transferred to the Base 
Program, the Base Program functions in a stand-alone mode to yield 
rapidly the effects of parameter variations on vibrations. 
Stand-alone operation is achieved with a small core requirement. 
The data flows between the Base Program and the External Programs 
provide the user with wide latitude to select applications tai- 
lored to user interest, aptitude, and convenience. Applications 
range from simple use of the Base Program to obtain the effects of 
parameter variations to more advanced configurations involving 
External Programs. Input to the Base Program is designed to 
provide an economical format requiring input only for components 
selected by the user. Results are displayed for only those 
components comprising the final system configuration. 
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Section 4 


COMPUTER PROGRAM EXECUTION 

The Base Program SIMVIB employs a substructure method for assemb- 
ling a dynamic configuration. The substructure assembly method is 
discussed in detail in Reference 1. The user specifies the number 
and kinds of components or substructures making up the dynamic 
system and the manner in which the separate components interact 
with each other through connection nodes . The computer program 
reads the input data for the system components, converts param- 
eters to the internal units (foot, pound, second, and radian), and 
prints out the component input data. The user may elect to read 
data for the rotor model from external sources such as the "Rotor 
Aeroelastic Stability Analysis" (E927 program) described in 
Reference 5 or the "Rotor Aeroelastic Analysis" (G400 program) 
described in Reference 3. The transformation matrix which relates 
the set of dependent coordinates for the component assembly to the 
final set of independent coordinates used to obtain the solution 
to the equations of motion is developed. After forming the labels 
needed to identify the output system coordinates, the computer 
program proceeds to calculate the mass, damping, and stiffness 
matrices of the individual components and assembles the various 
matrices together. The system generalized forces are specified by 
the user. The final dynamic system equations can now be solved by 
one of several methods selected by the user: eigensolutions (real 

or complex), steady state forced response, or time history. The 
final component responses are printed out and saved in a data file 
for execution of the "Plot Package Program" . Time history solu- 
tion results can be stored in a separate data file for restart 
conditions. 

Section 5 

BASE PROGRAM INTERFACING WITH EXTERNAL PROGRAMS 

The Base Program SIMVIB communicates with four external programs, 
as illustrated by the block diagram in Figure 1: G400 (Reference 

3), E927 (Reference 5), RIEVA (Reference 6) , and the Plot Package 
(Reference 7). The transfer of data is done through file units 
which are identified in Figure 2. 

5 . 1 G400 Rotor Aeroelastic Analysis 

The G400 Rotor Aeroelastic Analysis calculates the time history 
response of conventional articulated, non-articulated, teetering 
and composite bearingless rotors . The theoretical development is 
presented in Reference 2. Instructions to execute the program are 
found in Reference 3 and in this report in Appendix A, "G400/F389 
Coupled Program". The options available are uniform or variable 
inflow, steady or unsteady aerodynamic airfoil data, radial flow, 
blade mounted pendulum absorbers, higher harmonic control (HHC) of 
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blade pitch, and coupling with the F389 Rotor Induced Inflow 
Program. The program calculates the hub excitation vector for a 
fixed hub and the rotor impedance matrix relating 6 hub forces and 
moments to harmonic perturbations of the hub. The data are 
written to unit 11 and then loaded into the Base Program. The 
rotor hub excitation vector and impedance matrix are assembled 
with the mass, damping, and stiffness matrices of the internal 
components by activating component RE3 . The steady state forced 
response mode, which is invoked through component FRl, is used to 
calculate the coupled system responses. The parametric variation 
option, PVl, can be executed if changes in design characteristics 
need to be investigated. 

5.2 E927 Rotor Aeroelastic Stability Analysis 

The E927 Rotor Aeroelastic Stability Analysis prior to modifica- 
tions made for the contract provided a description of the sta- 
bility of coupled helicopter main or tail rotor/airframe systems. 
The analysis obtained linear stability characteristics for axial 
flow or hover conditions from calculations of the eigenvalues and 
eigenvectors of the coupled system. The analysis allows a maximum 
of 24 blade degrees of freedom: up to four blade coupled flat- 
wise/edgewise bending modes and two torsional modes plus blade 
rigid flapping and lead lag motion. Each blade has a collective 
and two cyclic rotor mode coordinates associated with it. The 
rotor hub is represented by five displacements: longitudinal, 
lateral, vertical, pitch, and roll. The trim rotor conditions may 
derive from G400 or a similar analysis. The theoretical develop- 
ment of the analysis is presented in Reference 4. The instruc- 
tions for executing the E927 program are available from Reference 
5, and Appendix B of the present report. 

For this contract the E927 program was modified to calculate only 
mass, damping, and stiffness matrices associated with rotor and 
hub coordinates. This specialization enabled the E927 model to be 
represented as a substructure and to be coupled to other sub- 
structures by means of the assembly method contained in the Base 
Program. The mass, damping, and stiffness matrices are trans- 
mitted to the Base Program through unit 10. The rotor matrices 
are then assembled with the matrices from the other system com- 
ponents by activating the RE2 component and the eigensolution EG2 
for systems with damping. The system stability is then determined 
from an inspection of the calculated eigenvalues. Steady state 
forced response can also be calculated if a rotor hub force and 
moment excitation vector is input to the Base Program by activat- 
ing components FRl and GFl. The excitation data would be provided 
by the G400 program or by any other suitable analysis. 

5.3 RIEVA - Empennage Excitation Analysis 


The RIEVA Empennage Excitation Analysis from Reference 6 calcu- 
lates the N/Rev (N is the number of blades) unsteady aerodynamic 
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loads induced by a rotor wake at an empennage which may be a 
vertical or a horizontal surface. The RIEVA utilizes rotor 
induced velocities from the F389 program through unit 15. Lifts 
and moments which affect the empennage are evaluated and then are 
input to the Base Program as harmonic coefficients through the 
substructure component GFl. The steady state forced response mode 
(FRl) is used to provide a solution for the final coupled system. 

5.4 Base Program Plot Package 

Results from the Base Program are stored in unit 2 for processing 
by the Base Program Plot Package. Execution of the plot program 
is discussed in detail in Reference 7. The program can display 
results for the different solution methods (eigensolutions, forced 
response, and time history). Examples of output plots are avail- 
able in this report in section 10, "Test Cases Descriptions and 
Results". 

5.5 Refinement of External Rotor Model 

The rotor model provided to the Base Program from the G400 program 
can be refined by executing the F389 Rotor Induced Flow Program 
and the WABAT (Wing and Body Aerodynamic Technique) Program. 

5.5.1 F389 Rotor Induced Flow Program 

This program calculates the , velocities induced by the 
rotor at the rotor disc or an empennage. It can be 
executed in an iterative manner with the G400 analysis 
to obtain the variable inflow field. Convergence is 
obtained when the blade circulations and induced veloc- 
ities are consistent from the two programs. Wake models 
which can be utilized are classical (skewed helix) or 
distorted wakes. The G400 and F389 programs communicate 
with each other by transfer of data through units 13, 16 
and 23 (G400 to F389) and unit 14 (F389 to G400). More 
details are in section 5.6, "Execution of the G400/F389 
Coupled Program". The F389 program is described in 
Reference 8. The F389 Program may also be used to 
calculate the rotor induced velocities at an empennage 
as input to the RIEVA Program. The data transfer is 
accomplished through unit 7 . No iteration between the 
F389 and RIEVA programs is performed. 

5.5.2 WABAT (Wing And Body Aerodynamic Technique) Program 

This program provides the fuselage induced velocities at 
the rotor disc to the F389 program through unit 15. No 
iteration between the two programs is performed. The 
WABAT geometry module is run before the potential flow 
program to define the fuselage geometry, as indicated in 
Figure 1. The WABAT plot package is executed to verify 


5 



that the fuselage geometry input is correct. Additional 
information on the WABAT program is available in Ref- 
erences 9 and 10. 

5.6 Execution of the G400/F389 Coupled Program 

The G400 and F389 programs have been coupled to allow the user 
three execution options: 

1. G400 only 

2. F389 only 

3 . G400/F389/G400 cycling operation 

The options are controlled through an additional input quantity 
specified on the first input card. A value of -2 specifies 
execution of the F389 program alone , -1 or zero are used to run 
G400 alone, while a positive number indicates the number of cycles 
of G400/F389 to perform. Appropriate input for the G400 or the 
F389 program (if run alone) is then loaded. In the uncoupled 
mode, the G400 .program calculates the time history response with 
uniform or variable inflow (calculated from a previous F389 run). 
Input needed to run the F389 program by itself can be punched out; 
this includes blade physical data, airfoil section character- 
istics, and non-induced axial velocities. If the coupled mode is 
exercised, the above input is passed through automatically from 
G400 to F389 using units 13, 16, and 23 as shown in Figure 2. 

The F389 program calculates the rotor induced velocities which are 
passed to the G400 program through unit 14 in the coupled mode or 
can be punched out in cards for input to G400 when run alone. 
When the coupled version is activated, a number of cycles can be 
specified; one cycle consists of execution of G400, F389, and 

G400. The accuracy of the rotor induced flow increases as more 
iterations are performed, however, the computer time also increas- 
es. One cycle normally requires 2 minutes of CPU (central proces- 
sing unit) time, while two cycles require about 3 and one-half 
minutes of execution time on the IBM 370/168 system. A trade-off 
between variable inflow accuracy and computer time usage and cost 
must be considered. For additional information on the execution 
of the G400/F389 Coupled Program refer to Appendix A, "G400/F389 
Coupled Program" . 


Section 6 

OVERVIEW OF BASE PROGRAM COMPONENTS 

This section provides an overview of the ten substructures pre- 
sently available in the Base Program: BFl, BF2, BMl, CNl, FAl, 

GFl, I SI, MSI, RE2, and RE3 . The four solution methods (EGl, EG2 , 
FRl, and THl) and the two operational modes (GEN and PVl) are also 
outlined. A summary of the attributes of the SIMVIB program 
components is provided in the last part. 
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6 . 1 Substructures Attributes 


The substructures available with the Base Program are discussed in 
detail in Reference 1. A summary of the substructures and a brief 
description is presented in Table 1 . Each substructure or com- 

ponent is identified by a three-character name and is assigned by 
the user a unique element number. The Base Program requires input 
data for only those components identified by the three-character 
name and element number. Several components may be present in the 
synthesis of a dynamic system. One component may be connected to 
another component by assigning the same connection node number 
among the two components. When components are connected, the 
displacements of the components at the point of connection are 
equal. The specification of connections enables a system to be 
synthesized from components . A component may be repeated several 
times. For example, a dynamic structure may be represented by 
five normal mode shapes (MSI). However, each MSI component must 
have a unique element number. The user should refer to section 7 
entitled "Base Program Component User Information" for examples of 
input data. 

6 ,2 Solution and Operational Modes Attributes 

The Base Program can employ four methods of solution to calculate 
the dynamic system responses. The methods are eigensolutions 
(real and complex), forced response, and time history. Two 
operational modes are available also. These are parametric 
variation and general control mode. The solution and operational 
modes are summarized in Table 2 which includes a brief description 
of each mode. If additional information is desired by the user. 
Reference 1 should be consulted. 

6 .3 Summary of Component Attributes 

A summary of the attributes of the Base Program components (10 
substructures, 4 solution modes, and 2 operational modes) is shown 
in Table 3. • This table provides the three-character name, the 
nvimber (NDLOAD) of component-dedicated input locations (see 
sections 9.7 and 9.8), the number (NWORK) of component-dedicated 
work locations (see sections 9.7 and 9.8), the number of connec- 
tion nodes which may be used with each substructure, the total 
number of coordinates, and the Fortran names associated with each 
substructure output coordinates. 

Section 7 

BASE PROGRAM INPUT DESCRIPTION 


This section contains information on the input data required to 
execute the SIMVIB program. A description of the input case may 
be provided through the use of up to four title cards. The input 
data for the individual components then follow. Each component is 
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at first identified by a three-character name and a unique element 
number. Component names and element numbers may be located in any 
order; the computer program reads the data, checks for duplication 
of element numbers, and processes the data to calculate the system 
response. The individual components are coupled to each other 
through the connection nodes . 

Each connection node has associated with it three Euleir angles 
defining the angular orientation of a substructure relative to a 
global axis system. When these angles are defined, one substruc- 
ture may be connected to another with an arbitrary angular orien- 
tation. All connection node displacements are resolved to global 
axes and the system solution is obtained in terms of coordinates 
resolved to global axes (see Reference 1). Figure 3 illustrates 
the angular disposition of any local axis system relative to the 
global system. For small Euler angle displacements, the Euler 
angles may be identified with pitch, roll, and yaw displacements. 
For large angles these names are not conventional. However, for 
the sake of convenience the Euler angles 0, 4>, and *1' are referred 
to as pitch, roll, and yaw angles, respectively. After the system 
response has been obtained, sxibstructure responses referred to 
substructure local axes are recovered and are displayed as line 
printer output or plot output. Coordinates which are not con- 
nected - e.g. absorber mass displacement - are invariant to 
rotations and have the same values in all axis systems . These 
properties of the coordinates are discussed in Reference 1 (see 
Substructure Assembly Method). 

The component data are loaded three items per card. The addi- 
tional space on the input cards is reserved for comments to 
provide a brief description of the input data. The input data for 
the various substructures, solution modes and operational modes 
employ the same format. 

Any solution mode may be invoked for a dynamic configuration. The 
only exception occurs when the rotor model component, RE3, from 
the G400 analysis is used. In this case, the only solution method 
applicable is the forced response solution (FRl). The operational 
mode, GEN, can be used for all cases. The parametric variation 
mode, PVl, cannot be activated when a time history response 
solution (THl) is performed. 

Several cases may be executed in succession (this is discussed in 
greater detail in the following section). Only revised or addi- 
tional input data for the present components need to be loaded. 
New components may not be introduced after the first case since 
computer storage for components data is allocated for the first 
case only. The final case must have the word STOP on the last 
card. 
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7.1 Input Data Format 


The card input format for the Base Program components is illus- 
trated in Figure 4. Three different loader formats are used for 
describing titles, component specification/multiple cases con- 
trol/end of run, and component characteristics as described below. 

1. Titles (card 1) 

a. Up to four titles can be specified. 

b. Each title is coded in columns 1 through 40. 

c. Read format is 4A10 for CDC and 10A4 for IBM. 

d. Titles are only read once for a complete run. 

e. If results are to be plotted, then the first 3 
titles will comprise the heading for each plot 
while the fourth title will describe the abscissa. 
For plotting purposes, 4 titles must be coded (the 
second and third title may be blank, cards). Column 
1 must be blank for all 4 titles when the Tektronix 
plot package is used. The fourth title is limited 
to 15 characters (columns 2 - 16). 

f. The exact image of the input titles is printed out 
before the component input data images are printed. 

2. Component Specification (card 2) 

a. The component name, NAMl, is coded in columns 1 
through 4, left adjusted. Format is A4. 

b. Each component is assigned a unique element number, 
NBl. It is coded in columns 10 through 12, right 
adjusted. Format is 13. 

c. If input data are to be repeated for a component, 
then COPY is coded in columns 18 through 21 and the 
previous element number, NB2 , is specified in 
columns 25 through 27 (right adjusted - 13 format). 
If no data are to be copied from a previous com- 
ponent, COPY and NB2 are not coded. 

d. Comments may be added in columns 40 through 79. 
Format is 4A10 for CDC and 10A4 for IBM. 

3. Component Characteristics (card 3) 

a. The input location of the first data item on the 
card, LOC, is located in columns 1 through 3, right 
adjusted. Format is 13. The last input data card 
for a given component has a minus sign in columns 1 
or 2. 

b. A delimeter comment, DEL, may be coded in columns 4 
through 6 to separate the location number and the 
input data. Format is A3 . 
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c. Up to 3 input values may be loaded per card. The 
format field is 3G10.3. 

d. Comments may be added in columns 40 through 79, 

4. . Multiple Cards Control (card 4) 

a. Card loader format is similar to that described in 
section "Component Specification" above. 

b. To run multiple cases, the control flags END and 
CASE are coded exactly as shown on card 4 in Figure 
4 after the last input data card for the case. 

c. The case number, NCS, is optional. 

5. End of Run (card 5) 

The last input data card contains the control flag STOP 

coded in columns 1 through 4. 

The loader is designed to allow the user to supersede data in a 
card by inserting an additional card which contains new data. 
This card must be placed further down in the input stream than the 
card whose data are to be replaced, and it must be placed in a 
data region applicable to the component whose data are to be 
modified. One or more cards may be modified. More than one card 
may be used to provide a succession of modifications in the input 
stream. This modification capability is designed to provide 
flexibility in input preparation. A restriction is that if a new 
value is zero, this new value must always be placed in columns 7 
through 16 (first data field, Figure 4). Two examples of loader 
data are presented in Figures 5 and 6 . 

The first run in Figure 5, represents a forced response study 
(FRl) of a system including a fixed system absorber (FAl), and a 
single normal mode (MSI) being excited by a longitudinal (x) force 
(GFl) at a frequency range from 1.59 to 9.55 hz (PVl). The output 
printout is not suppressed and the variable stored for three- 
dimensional plotting is the damping ratio of the fixed system 
absorber (GEN). For plots, the abscissa is labelled "FREQUENCY 
(HZ)", as indicated by the fourth title card. Only one case is 
executed in this example. 

The running of a second case with a new value of damping ratio in 
FAl location 3 and GEN location 10 would cause new results to be 
stored for the second damping ratio. The specification of a 
damping value in GEN location 10 causes this value to be stored 
for plotting. Successive cases similar to those described may be 
run to store results corresponding to additional damping ratios. 
A three-dimensional plot generated by the plot program consists of 
two axes of abscissas and one axis of ordinates . In the sample 
problem, the first axis of abscissas is the axis of frequencies. 
The second axis of abscissas is the axis of dampings . The axis of 
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ordinates is the axis of system responses. A curve of variation 
of system responses with frequency is displayed for each value of 
damping. These curves are disposed on the second axis of ab- 
scissas at points defined in GEN location 10. If the value in GEN 
location 10 is not loaded (blank string), the plot package de- 
faults to displaying points on the second axis of abscissas 
equidistantly. 

The second run in Figure 6 calculates the time history response 
(THl) for the system described above employing a second normal 
mode (MSI) introduced by the COPY command. Three cases are 
executed using the END CASE commands. In the second case, the 
fixed system absorber mass has been changed from 10 to 15 slugs; 
the mass is increased to 20 slugs for the third case. The print- 
out of the output results is eliminated by the control switch in 
GEN. For plots, the abscissa is labelled "TIME (SEC)". 

Additional inputs are needed if elastic rotor representations 
(component RE2 from the E927 program and component RE3 from the 
G400 program) are utilized in a computer run. The E927 rotor 
mass, damping, and stiffness matrices are written to unit 10 from 
an independent execution of the E927 program. The Base Program 
reads the rotor data from unit 10 in subroutine LK2RE2 and, after 
conversion to the proper internal units, they are printed out in 
subroutine ASMMCK, if the debug switch is activated. A sample 
output of the E927 matrices appears in Figure 7 for a rotor/ 
fuselage system with 17 degrees of freedom. The rotor impedance, 
higher harmonic control (HHC), and blade stress matrices, and the 
force vector from the G400 program are input to the Base Program 
through unit 11 . Unit 12 may be used to transmit user input data 
defining the impedance properties and force vector of a general 
substructure. A sample input from unit 12 is shown in Figure 8. 
The data are read in subroutine LK2RE3 and printed in PRTRE3 and 
PRTSTR. A sample output is shown in Figure 9. No controls are 
provided to suppress the input force excitation vector and rotor 
impedance matrix display. 

7.2 Component Descriptions, Input Data Specifications, and 

Output Coordinate Labels 

Detailed descriptions are given in this section of the functions, 
assumptions, and coordinates employed by components available in 
the Base Program. Components comprise substructures, solution 
modes, and operational modes. Input specifications are provided 
for each component. Appendix B in Reference 1 defines the equa- 
tions of motion for each substructure. Reference 1 also describes 
the mathematical basis of the solution modes (see section Solution 
Modes ) . 
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7.2.1 


Horizontal Linear Bifilar BFl 


Horizontal linear bifilar BFl is a bifilar vibration 
absorber designed to reduce N per rev vibrations in the 
plane of rotation of a helicopter rotor. The bifilar 
absorber is assumed to consist of N arms equal in number 
to the rotor blades. Each arm holds a pendular mass 
free to lead or lag. The arms are assumed to be equally 
disposed and each holds an equal absorber mass. The 
bifilar absorber rotates at the helicopter rotor speed. 
Figure 10 is a schematic of one arm of the bifilar. 
Each bifilar mass lead-lag displacement may be expressed 
as the sum of bifilar lead-lag collective and cyclic 
displacements . 

The BFl substructure employs 8 coordinates consisting of 
bifilar mass sine (y ) and cosine (y ) cyclic amplitudes 
and 6 hub displacemlnts . The reduced equation set has 
constant coefficients. Restrictions on the configura- 
tions that may be treated are that they must have three 
or more equally disposed and equal masses . Pendulum 
perturbations are assumed to be small enough to permit 
linearization of the equations of motion. Solutions are 
obtained in terms of the cyclic amplitudes y and y . 
Solution modes that may be used with BFl are ‘trI, EG?, 
EG2, and THl. Initial values of bifilar mass cyclic 
mode amplitudes may be specified for transient analysis. 

HORIZONTAL LINEAR BIFILAR BFl 
INPUT DATA 


Location Fortran 

No. Name Input Item Description Units 


1 

NCNl 

Connection node number 

nd 

2 

N 

Total number of bifilar masses 

nd 

3 

MB 

Bifilar mass 

slugs 

4 

RR 

Distance from center of bifilar 
tracking hole to center of rotation 

in 

5 

R 

Equivalent pendulum arm 

in 

6 

2ETGAM 

Bifilar damping ratio 

nd 

7 

OMEGA 

Bifilar rotational speed 

rad/sec 

8 

THETA 

Euler pitch angle at node (0 ) 

deg 

9 

PHI 

Euler roll angle at node (0 ^ 

deg 

10 

XSI 

Euler yaw angle at node (0 ^ 

deg 

11 

GAMMAC 

Bifilar mass cosine cyclic^initial 
amplitude 

deg 

12 

DGAMMC 

Bifilar mass cosine cyclic initial 
velocity 

deg/sec 

13 

GAMMAS 

Bifilar mass sine cyclic initial 
amplitude velocity 

deg/sec 

14 

DGAMMS 

Bifilar mass sine cyclic initial 
velocity 

deg/sec 


\ 
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HORIZONTAL LINEAR BIFILAR BFl 
COORDINATE LABELS 


Coordinate Description 

Symbol Code Output 


VC 

GAMMAC 

GAMC 

Bifilar mass cosine cyclic lead 
amplitude 

vs 

GAMMAS 

GAMS 

Bifilar mass sine cyclic lead 
amplitude 

u 

X 

X 

Hvib 

X linear displacement 

V 

Y 

Y 

Hub 

y linear displacement 

w 

Z 

Z 

Hub 

z linear displacement 


THTX 

THTX 

Hub 

roll displacement 

®2 

THTY 

THTY 

Hub 

pitch displacement 


THTZ 

THTZ 

Hub 

yaw displacement 

7.2.2 

Vertical 

Linear Bifilar 

BF2 



Vertical linear bifilar BF2 is a bifilar vibration 
absorber designed to reduce N per rev vibrations per- 
pendicular to the plane of rotation of a helicopter 
rotor. The bifilar absorber is assumed to consist of N 
arms equal in ntunber to the rotor blades. Each arm 
holds a pendular mass free to flap up or down. The arms 
are assumed to be equally disposed and each holds an 
equal absorber mass. The bifilar absorber rotates at 
the helicopter rotor speed. Figure 11 is a schematic of 
one arm of the bifilar. Substructure BF2 employs 9 
coordinates consisting of bifilar mass collective (p ), 
sine cyclic (p ), and cosine cyclic (P_) amplitudes, and 
6 hub displacements. Restrictions on the configurations 
that may be treated are that they must have three or 
more equally disposed and equal masses . Solution modes 
that may be used with BF2 are FRl , EGl , EG2 , and THl . 
Initial values of bifilar mass collective and cyclic 
amplitudes may be specified for transient analysis. 
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VERTICAL LINEAR BIFILAR BF2 
INPUT DATA 


Location 

No. 


Fortran 

Name 


Input Item Description 


Units 


1 

NCNl 

Connection node number 

nd 

2 

N 

Total number of bifilar masses 

nd 

3 

MB 

Bifilar mass 

slugs 

4 

RR 

Distance from center of bifilar 
tracking hole to center of rotation 

in 

5 

R 

Equivalent pendulum arm 

in 

6 

ZETBET 

Bifilar damping ratio 

nd 

7 

OMEGA 

Bifilar rotational speed 

rad/sec 

8 

THETA 

Euler pitch angle at node (0 ) 

deg 

9 

PHI 

Euler roll angle at node (6 Y 

deg 

10 

XSI 

Euler yaw angle at node (0 ^ 

deg 

11 

BETAO 

Bifilar mass collective initial 
amplitude 

deg 

12 

DBETAO 

Bifilar mass collective initial 
velocity 

deg/sec 

13 

BETAC 

Bifilar mass cosine initial amplitude 

deg 

14 

DBETAC 

Bifilar mass cosine initial velocity 

deg/sec 

15 

BETAS 

Bifilar mass sine initial amplitude 

deg 

16 

DBETAS 

Bifilar mass sine initial velocity 

deg/sec 


Coordinate 
Figure Code 

VERTICAL LINEAR BIFILAR BF2 
COORDINATE LABELS 

Description 

Output 

P.O 

BETAO 

BETO 

Bifilar mass collective 
flapping amplitude 

Ps 

BETAC 

BETC 

Bifilar mass cosine cyclic 
flapping amplitude 

Pc 

BETAS 

BETS 

Bifilar mass sine cyclic 
flapping amplitude 

X 

X 

X 

Hub X linear displacement 

Y 

Y 

Y 

Hub y linear displacement 

Z 

Z 

Z 

Hub z linear displacement 

Qi 

THTX 

THTX 

Hub roll displacement 

02 

THTY 

THTY 

Hub pitch displacement 

03 

THTZ 

THTZ 

Hub yaw displacement 
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Substructure BMl 


7.2.3 


Beam Model 1, BMl, is a uniform elastic beam. The 
substructure has 12 coordinates consisting of 6 linear 
and angular displacements at each of the two ends of the 
beam. The beam may be connected to other substructures 
at its ends. The unifo:rm mass of the beam, m , is 
approximated by two concentrated masses equal to m^2 at 
its ends . The beam is assumed to have zero mass moment 
of inertia associated with bending and torsion rota- 



tional 

displacements. Structural damping is assumed to 


be zero 

. Figure 12 illustrates BMl properties . 




BEAM MODEL BMl 




INPUT DATA 


Location 

Fortran 



No. 

Name 

Input Item Description 

Units 

1 

NCNl 

Connection mode number at end 1 

nd 

2 

NCN2 

Connection mode number at end 2 

nd 

3 

E 

Young ’ s modulus 

psi 

4 

A 

Cross sectional area 

ft2 

5 

lYA 

Second area moment about y-axis 

ft^ 

6 

IZA 

Second area moment about z-axis 

ft^ 

7 

GJ 

Torsional stiffness slug-ft/sec 

8 

MO 

Beam mass 

slugs 

9 

IXX 

Second mass moment about x-axis 

slug-ft2 



There is no provision in this pro- 
gram version for representation of 
rotary inertias in the BMl equations 
of motion. 


10 

lYY 

Second mass moment about y-axis slug-ft^ 

11 

IZZ 

Second mass moment about z-axis slug-ft^ 



There is no provision in this program 
version for representation of c.g. 
offset terms in the BMl equations of 
motion. 


12 

UCG 

Offset of beam c.g. from elastic 




axis in the x-direction 

ft 

13 

VCG 

Offset of c.g. in the y-direction 

ft 

14 

WCG 

Offset of c.g. in the z-direction 

ft 

15 

LENGTH 

Beam length 

ft 

16 

THETAl 

Euler pitch angle at node 1 ( yl ) 

deg 

17 

PHIl 

Euler roll angle at node 1 ( xl) 

deg 

18 

XSIl 

Euler yaw angle at node 1 ( zl ) 

deg 
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Input values of Euler angles at 
node 2 must be equal to the values 
at node 1 . 


19 THETA2 Euler pitch angle at node 2 ( y2 ) deg 

20 PHI2 Euler roll angle at node 2 ( x2 ) deg 

21 XSI2 Euler yaw angle at node 2 ( z2) deg 


Coordinate 
Symbol Code 

BEAM MODEL BMl 
COORDIEATE LABELS 

Description 

Output 

Ui 

XI 

XI 

Node 

1 

X displacement 

Vi 

Y1 

Y1 . 

Node 

1 

y displacement 

Wi 

Z1 

Z1 

Node 

1 

z displacement 

(0i)i 

THTXl 

THXl 

Node 

1 

roll displacement 

(02 ) 1 

THTYl 

THYl 

Node 

1 

pitch displacement 

(03 ) 1 

THTZl 

THZl 

Node 

1 

yaw displacement 

U2 

X2 

X2 

Node 

2 

X displacement 

V2 

Y2 

Y2 

Node 

2 

y displacement 

W2 

Z2 

Z2 

Node 

2 

z displacement 

(01 )2 

THTX2 

THX2 

Node 

2 

roll displacement 

( 02)2 

THTY2 

THY2 

Node 

2 

pitch displacement 

( 03)2 

THTZ2 

THZ2 

Node 

2 

yaw displacement 


7.2.4 Constraint Substructure CNl 


Constraint substructure CNl may be used to eliminate or 
constrain coordinates at a connection node . There are 
six coordinates at each connection node (see Figure 13). 
The coordinates consist of three linear displacements 
and three angular displacements . The user may constrain 
any five of the six coordinates at a connection node. 
At least one coordinate among the six must be retained. 

Constraint substructure CNl may not be used in conjunc- 
tion with modal substructure MSI at the same connection 
node. CNl may be used to eliminate physical coordinates 
at a node. Retained physical coordinates cannot then be 
replaced by modal coordinates. Simultaneous elimination 
of physical connection node coordinates and replacement 
of retained coordinates by modal coordinates is accom- 
plished by using zero mode shapes in an MSI substructure. 

Non-zero Euler angles (locations 8 to 10) are used to 
constrain local coordinates; zero values constrain 
global coordinates. 
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CONSTRAINT SUBSTRUCTURE CNl 
INPUT DATA 


Location 
No. . 

Fortran 

Name 

Input Item Description 

Units 

1 

2-7 

NCN 

CONSTR 

Connection node mumber 

Constraint control switch 

If = 0, then degree of freedom (d.o.f 

is constrained at connection node 

If = 1, then d.o.f. is not restrained 

Loc 2 controls x d.o.f. at node 

Loc 3 controls y d.o.f. at node 

Loc 4 controls z d.o.f. at node 

Loc 5 controls 0^ d.o.f. at node 

Loc 6 controls 0^ d.o.f. at node 

Loc 7 controls 0^ d.o.f. at node 

z 

nd 

nd 

.) 

8 

9 . 

10 

THETA 

PHI 

XSI 

Euler pitch angle (0 ) at node 

Euler roll angle at node 

Euler yaw angle (0_7 at node) 

2 

deg 

deg 

deg 


Constraint substructure CNl coordinate output labels are 
the same as the coordinate labels obtained from the 
component constrained by CNl. If constraint coordinate 
CNl is used to constrain 5 BMl connection coordinates at 
a connection node, and the remaining coordinate is w at 
node 1, the output label of this coordinate will be . 

7.2.5 

Eigensolution for Systems Without Damping EGl 



This eigensolution mode yields the frequencies and mode 
shapes of structures without damping. Matrices must be 
constant w.r.t. time. The effects of changes in 

physical properties on mode shapes and frequencies may 
be determined. These effects may be used to aid the 
interpretation of the effects of design variables on 
normal modes and vibratory response . Section 8 . 3 

describes the line printer output displayed by EGl. 



COMPONENT EGl INPUT DATA 


Location 

No. 

Fortran 

Name 

Input Item Description 

Units 

1 

I DEBUG 

Set = 0 for no debug printout - 
default value 

Set = 1 display M and K matrices 
in substructure and system equations 
of motion 

nd 
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7.2.6 Eiqensolution for Systems with Damping EG2 

The eiqensolution mode for systems with damping is used 
to provide information on the stability of coupled 
systems . These systems may include rotor representa- 
tions which embody dynamic and aerodynamic matrices 
derived from rotor aeroelastic analysis E927. Stability 
information may be used to supplement the steady state 
vibratory response data for a system to provide more 
data on system attributes . Section 8 . 3 describes the 
line printer output displayed by EG2. 

COMPONENT EG2 INPUT DATA 

Location Fortran 

No. Name Input Item Description Units 


1 IDEBUG Set = 0 for no debug printout - nd 

default value 

Set = 1 display M, C, and K matrices 
in substructure and system equations 
of motion. 


7.2.7 Substructure FAl 

The function of substructure FAl is to simulate a 
vibration absorber which may be attached to any other 
substructure in the fixed (non-rotating) system. 
Substructure FAl has seven coordinates . These are 
absorber mass displacement and six displacements of the 
base consisting of three linear displacements and three 
angular displacements. The mass of the absorber is free 
to move perpendicular to the base. The mass of the 
absorber is not free to move parallel to the base. 
Viscous damping is assumed. Figure 14 illustrates 
substructure FAl properties. Initial mass displacement 
and velocity may be specified for transient response 
analysis. 
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FIXED ABSORBER FAl INPUT DATA 


Location 

No. 

Fortran 

Name 

Input Item Description 

Units 

1 

NCN 

Connection node number 

nd 

2 

MO 

Absorber mass 

slugs 

3 

ZETA 

Absorber damping ratio 

nd 

4 

KO 

Absorber stiffness 

Ib/ft 

5 

IXX 

Second mass moment about x-axis 

slug-ft2 

6 

lYY 

Second mass moment about y-axis 

slug-ft2 

7 

IZZ 

Second mass moment about z-axis 

slug- ft2 

8 

THETA 

Euler pitch angle (0 ) at node 

deg 

9 

PHI 

Euler roll angle (8..^ at node 

deg 

10 

XSI 

Euler yaw angle (0 ; at node 

deg 

11 

DELTA 

Initial displacement of absorber 




mass 

ft 

12 

DDELTA 

Initial velocity of absorber mass 

ft/sec 


FIXED ABSORBER FAl COORDINATE LABELS 


Coordinate 

Description 


Symbol 

Code 

Output 



u 

DELTA 

X 

DELT 

X 

Absorber mass displacement 
Base X linear displacement 

V 

Y 

Y 

Base 

y linear displacement 

w 

Z 

Z 

Base 

z linear displacement 

01 

THTI 

THTX 

Base 

roll displacement 

02 

THT2 

THTY 

Base 

pitch displacement 

03 

THT3 

THTZ 

Base 

yaw displacement 


7.2.8 Steady State Forced Response FRl 

The harmonic steady state forced response mode of 
calculation FRl yields the harmonic coefficients of 
response to a steady harmonic excitation. The system 
equation reduces to an algebraic equation which is 
solved to obtain the harmonic coefficients of response. 
Speed of calculation in the FRl mode is high and the 
method is suitable for parametric design studies with 
component PVl. Section 8.3 describes the line-printer 
output displayed by FRl . 
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Location 

No. 

Fortran 

Name 

COMPONENT FRl INPUT DATA 
Input Item Description 

Units 

1 

OMEGA 

Forcing frequency 

hz 

2 

IDEBUG 

Debug control switch 
Set = 0 for no debug printout - 
default value 

Set = 1 to display M, C, and K 
matrices and force vectors for sub- 
structures and assembly, and solution 
vector for assembly at solution step. 

nd 

3 

ICONVG 

If = 0, results are displayed as 
displacements (ft) 

If = 1, results are displayed as 
accelerations (g's) 

nd 


7.2.9 General Control Mode, GEN 

The general control mode, GEN, is designed to provide a 
basis for the application of general controls to a case. 
In the present version of the Base Program, GEN provides 
input which may be used to suppress line printer output. 
When line printer output is suppressed, only the input 
echo, and coordinate solution list for one iteration 
(parameter variation or time step) are displayed. 
Solution data are stored on the plot file. The coordi- 
nate list is then used to identify, variables to be 
plotted. Suppression of line printer output minimizes 
hard copy line printer output. A further control now 
present in GEN permits a numerical value to be specified 
for an independent variable on the first axis of ab- 
scissas. The second axis of abscissas may contain the 
value of the mass of a vibration absorber varied in a 
design study. The axis of ordinates would contain 
vibratory response of the system. If the variable 
occurring on the second axis of abscissas is not speci- 
fied, the plot package displays curves corresponding to 
points on the second axis equidistantly . If the GEN 
string is not loaded, the GEN mode is not invoked and 
its input is not required. 
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COMPONENT GEN INPUT DATA 


Location 

No. 

1 


2 

3-9 

10 


Fortran 

Name Input Item Description Units 


CNTLl Final results control switch nd 

If = 0, suppress full line printer 
output-default value 
If = 1, get line printer output 

When output is suppressed a coordinate list of 
coordinate numbers and names is displayed for 
the first solution step. Coordinate numbers 
are tabulated against element number and 
coordinate name. The values of coordinate 
responses are stored for plotting at each 
solution step and are printed only at the 
first step. When coordinates are substructure 
connection node displacements, these dis- 
placements are referred to substructure local 
axes. Unconnected coordinates - e.g. vibra- 
tion absorber mass displacement - and modal 
coordinates are in-variant to axis rotations. 
The first solution step may yield the response 
vector obtained for the first value of a 
parameter in a parameter variation mode or for 
the first time instant in a time history 
solution. When the plot package is invoked, 
the user uses the first step coordinate list 
to identify and key in the coordinate numbers 
of the coordinates that are to be plotted. 

CNTL2 Input narrative control swtich 

If = 0, suppress printout - default value 

If = 1, display input narrative 

Open locations for future use 

XINDEP Independent variable for 3-dimensional plots. 
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7.2.10 


A three-dimensional plot generated by the plot 
program is based on two axes of abscissas and 
one axis of ordinates. The axis of ordinates 
contains the values of the responses of a 
coordinate of the system (referred to sub- 
structure local axes in the case of connection 
node displacements ) . A curve of variation of 
response with values along the first axis of 
abscissas is placed on the second axis at a 
position corresponding to the value in loca- 
tion 10. If a blank string is loaded in lo- 
cation 10 the plot package defaults to dis- 
playing curves equidistantly on the second 
axis of abscissas. 

Substructure GFl 

Generalized Force 1, GFl, permits the specification of a 
harmonic force and moment excitation acting at a point 
on a structure. The GFl substructure has six coordi- 
nates, comprising three linear and three angular dis- 
placements at a connection node. Forces and moments are 
resolved to a non-rotating axis system. The user 
specifies the cosine and sine coefficients of these 
forces and moments to define a harmonic excitation of 
known magnitude. The mass, damping, and stiffness 
matrices for substructure GFl are null matrices. 
Substructure GFl may be used to represent an excitation 
for forced response (FRl) and time history (THl) solu- 
tion modes . Substructure GFl may be included in eigen- 
solution (EGl and EG2) input streams without affecting 
the final eigensolution. Figure 15 illustrates GFl 
properties . Substructure GFl may be used to define a 
point at which vibrations are to be minimized by appli- 
cation of higher harmonic control inputs . The HHC mode 
is invoked by setting IHRESP =1. In this mode, a 
weighted sum of the six responses at the connection node 
to which GFl is attached is minimized. Weighting 
factors are placed in GFl loader locations following the 
HHC flag IHRESP. The performance index to be minimized 
and the basis of the minimization method is described in 
Reference 1. When the HHC mode is invoked, it is neces- 
sary to employ substructure RE3 to read HHC matrices 
transmitted to the Base Program (SIMVIB). Rotor aero- 
elastic analysis G400 must be run in a stand alone mode 
to place rotor impedance and HHC data on files. Cases 
14 to 16 in Table 13 illustrate the use of GFl without 
and with HHC. 
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GENERALIZED FORCE GFl INPUT DATA 


Location 

Fortran 



No. 

Name 

Input Item Description 

Units 

1 

NCN 

Connection node number 

nd 

2 

THETA 

Euler pitch angle (0 ) 

deg . 

3 

PHI 

Euler roll angle (0 ^ 

deg 

4 

XSI 

Euler yaw angle (0 ^ 

deg 

5 

FXCOS 

Coefficient of cosfne term in harmonic 
force excitation in x direction 

lb 

6 

FXSIN 

Coefficient of sine term in harmonic 
force excitation in x direction 

lb 

7 

FYCOS 

Coefficient of cosine term in harmonic 
force excitation in y direction 

lb 

8 

FYSIN 

Coefficient of sine term in harmonic 
force excitation in y direction 

lb 

9 

FZCOS 

Coefficient of cosine term in harmonic 
force excitation in z direction 

lb 

io 

FZSIN 

Coefficient of sine term in harmonic 
force excitation in z direction 

lb 

11 

FTICOS 

Coefficient of cosine term in harmonic . 
moment excitation in x direction 

ft-lb 

12 

FTISIN 

Coefficient of sine term in harmonic 
moment excitation in x direction 

ft-lb 

13 

FT2COS 

Coefficient of cosine term in harmonic 
moment excitation in y direction 

ft-lb 

14 

FT2SIN 

Coefficient of sine term in harmonic 
moment excitation in y direction 

ft-lb 

15 

FT3COS 

Coefficient of cosine term in harmonic 
moment excitation in z direction 

ft-lb 

16 

FT3SIN 

Coefficient of sine term in harmonic 
moment excitation in z direction 

ft-lb 

17 

IHRESP 

Control switch for HHC option 
= 0 do not activate option 
= 1 activate option to minimize a 

nd 



weighted sum of responses at the 
HHC connection node by application 
of HHC control inputs 


18 

WZX 

Weighting factor for X response in 
HHC mode 

nd 

19 

WZY 

Weighting factor for Y response in 
HHC mode 

nd 

20 

WZZ 

Weighting factor for Z response in 
HHC mode 

nd 

21 

WZWX 

Weighting factor for 0^ response in 
HHC mode 

nd 

22 

WZWY 

Weighting factor for 0 y response in 
HHC mode 

nd 

23 

WZWZ 

Weighting factor for 0„ response in 
HHC mode ^ 

nd 
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GENERALIZED FORCE GFl COORDINATE LABELS 


Coordinate 
Symbol Code • 

Output 


Description 

u 


X 

X 

Node 

X linear displacement 

V 


Y 

Y 

Node 

y linear displacement 

w 


z 

Z 

Node 

z linear displacement 

0 1 


THTl 

THTX 

Node 

roll displacement 

02 


THT2 

THTY 

Node 

pitch displacement 

03 


THT3 

THTZ 

Node 

yaw displacement 

7.2 

.11 

Vibration 

Isolator ISl 




Vibration isolator ISl is a nodal isolator which may be 
placed between an excitation source, like the rotor/ 
transmission system, and a fuselage to minimize vibra- 
tions at the fuselage. In its essential form, the 
isolator employs a pivoted bar weight to balance spring 
and inertia forces to create a node or point of zero 
vibration motion at the fuselage attachment. When 
damping is present, vibrations are minimized at the 
fuselage attachment at the excitation frequency (N per 
rev for an N bladed rotor). Nodal isolator ISl employs 
18 coordinates . Figure 16a illustrates a physical 
isolator and Figure 16b is a schematic corresponding to 
the physical isolator. The schematic shows two anti- 
resonant bar weights. The lower bar weight is a math- 
ematical bar weight permitting additional design flexi- 
bility. It has no counterpart in Figure 16a. There are 
two attachment (connection) nodes. Each connection node 
is allowed to have three translational and three rota- 
tional displacements,. The upper connection node A would 
be. attached to a rotor/transmission system. The lower 
connection node B would be attached to a fuselage; this 
is the point where vibrations are to be minimized. Each 
antiresonant bar weight has three translational degrees 
of freedom. Antiresonant bar flexibilities and bar 
dampings may be introduced by modifying the effective 
mass of the antiresonant bar weight with factors employ- 
ing the bar natural frequencies and dampings (locations 
50 to 55). 
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COMPONENT I SI INPUT DATA 


Location 

No. 

Fortran 

Name 

Input Item Description 

Units 

1 

NCNl 

Connection node number 1 

nd 

2 

NCN2 

Connection node number 2 

nd 

3 

AL 

Distance between node 1 & pivot 2 

in 

4 

BL 

Distance between pivot 2 & Y-Z 
isolator 

in 

5 

CL 

Distance between pivot 2 & node 2 

in 

6 

DL 

Distance between pivot 2 & X isolator 

in 

7 

W1 

Anti-resonant bar weight of Y-Z 
isolator (mass 1) 

lb 

8 

XI lY 

Y-Z isolator mass moment about y-axis 

in-lb-sec2 

9 

XI12 

Y-Z isolator mass moment about z-axis 

in-lb-sec2 

10 

W2 

Anti-resonant bar weight of X 
isolator (mass 2) 

lb 

11 

XI2Y 

X isolator mass moment about y-axis 

in-lb-sec2 

12 

TKTX 

Stiffness between pivot 1 & node 1 
in X direction 

Ib/in 

13 

TKTY 

Stiffness between pivot 1 & node 1 
in y direction 

Ib/in 

14 

TKTZ 

Stiffness between pivot 1 & node 1 
in 2 direction 

Ib/in 

15 

TKEX 

Stiffness between pivot 1 & node 2 
in X direction 

Ib/in 

16 

TKEY 

Stiffness between pivot 1 & node 2 
in y direction 

Ib/in 

17 

TKEZ 

Stiffness between pivot 1 & node 2 
in z direction 

Ib/in 

18 

TKAX 

Stiffness between pivot 2 & node 2 
in X direction 

Ib/in 

19 

TKAY 

Stiffness between pivot 2 & node 2 
in y direction 

Ib/in 

20 

TKAZ 

Stiffness between pivot 2 & node 2 
in z direction 

Ib/in 

21 

RKTHB 

Rotational stiffness of node 2 about 
y-axis 

in-lb/rad 

22 

RKTHT 

Rotational stiffness of pivot 2 about 
y-axis 

in-lb/rad 

23 

RKPHT 

Rotational stiffness of pivot 2 about 
z-axis 

in-lb/rad 

24 

TCTX 

Damping ratio of damper in parallel 
with spring TKTX 

nd 

25 

TCTY 

Damping ratio of damper in parallel 
with spring TKTY 

nd 

26 

TCTZ 

Damping ratio of damper in parallel 
with spring TKTZ 

nd 
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27 

TCEX 

Damping ratio of damper in parallel 
with spring TKEX 

nd 

28 

TCEY 

Damping ratio of damper in parallel 
with spring TKEY 

nd 

29 

TCEZ 

Damping ratio of damper in parallel 
with spring TKEZ 

nd 

30 

TCAX 

Damping ratio of damper in parallel 
with spring TKAX 

nd 

31 

TCAY 

Damping ratio of damper in parallel 
with spring TKAY 

nd 

32 

TCAZ 

Damping ratio of damper in parallel 
with spring TKAZ 

nd 

33 

RCTHB 

Damping ratio of spring in parallel 
with RKTHB 

nd 

34 

RCTHT 

Damping ratio of spring in parallel 
with RKTHT 

nd 

35 

RCPHT 

Damping ratio of spring in parallel 
with RKPHT 

NOTE: Damping is calculated from 

TCTX1=2 * TCTX * TKTX/WTX, lb sec/in 
If WTX=0, then TCTX1=0. 

Locations 36 to 47 permit the user 
to input values of frequencies re~ 
guired in the formulas for dampings 
contributed by damper elements. 

nd 

36 

WTX 

Frequency required to define value 
of damper in parallel with TKTX 

hz 

37 

WTY 

Frequency required to define value 
of damper in parallel with TKTY 

hz 

38 

WTZ 

Frequency required to define value 
of damper in parallel with TKTZ 

hz 

39 

WEX 

Frequency required to define value 
of damper in parallel with TKEX 

hz 

40 

WEY 

Frequency required to define value 
of damper in parallel with TKEY 

hz 

41 

WEZ 

Frequency required to define value 
of damper in parallel with TKEZ 

hz 

42 

WAX 

Frequency required to define value 
of damper in parallel with TKAX 

hz 

43 

WAY 

Frequency required to define value 
of damper in parallel with TKAY 

hz 

44 

WAZ 

Frequency required to define value 
of damper in parallel with TKAZ 

hz 

45 

WTHB 

Frequency required to define value 
of damper in parallel with RKTHB 

hz 

46 

WTHT 

Frequency required to define value 
of damper in parallel with RKTHT 

hz 
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47 WPHT Frequency required to define value hz 

of damper in parallel with RKPHT 

48 I3D2D Isolator type flag. Set = 2 for 

2-D or Y-Z isolator type. Set = 3 
for 3-D or Y-2 and X isolator type. 

When the flag is set to 2 , the 
terms contributed by the second 
(lower) anti-resonant bar weight 
are removed from the isolator 
equations of motion. When the 
flag is set to 3, both anti- 
resonant bar weights are repre- 
sented. 

49 IFLEX Isolator bar flexibility effects nd 

option. If = 0, no flexibility 
effects. If = 1, include flex- 
ibility effects. 

Note; isolator bar flexibility 
effects are applied to the isola- 
tor weights and mass moments of 
inertia using the amplification 
factor calculated from the fre- 
quency ratio ( forcing/natural ) and 
critical damping values input in 
the following locations . 

50 RXYMl Y-Z isolator bar frequency ratio in nd 

X-Y plane 

51 RXZMl Y-Z isolator bar frequency ratio in nd 

X-Z plane 

52 RXZM2 X isolator bar frequency ratio in nd 

X-Z plane 

53 DXYMl Y-Z isolator bar critical damping in nd 

X-Y plane 

54 DXZMl Y-Z isolator bar critical damping nd 

in X-Z plane 

55 DXZM2 X isolator bar critical damping nd 

in X-Z plane 

56 THETAl Euler pitch angle at connection deg 

node 1 (0yl) 

57 PHIl Euler roll angle at connection node deg 

1 (0x1) 

58 XSIl Euler yaw angle at connection node 1 deg 

(0zl) 

59 THETA2 Euler pitch angle at connection node deg 

2 (0y2). 

The Euler angles at connection node 
2 must be equal to the Euler angles 
at connection node 1 . 


27 



* t. 


60 

PHI 2 

Euler roll angle at connection 
node 2 (0x2) 

deg 

61 

XSI2 

Euler yaw angle at connection node 
2 (0z2) 

deg 

62 

XT 

Initial X displacement of anti- 
resonant bar weight 1 

in. 

63 

DXT 

Initial X velocity of anti- 
resonant bar weight 1 

in/sec 

64 

YT 

Initial Y displacement of anti- 
resonant bar weight 1 

in 

65 

DYT 

Initial Y velocity of anti- 
resonant bar weight 1 

in/sec 

66 

ZT 

Initial Z displacement of anti- 
resonant bar weight 1 

in 

67 

DZT 

Initial Z velocity of anti-resonant 
bar weight 1 

in/sec 

68 

XB 

Initial X displacement of anti- 
resonant bar weight 2 

in 

69 

DXB 

Initial X velocity of anti-resonant 
bar weight 2 

in/sec 

70 

YB 

Initial Y displacement of anti- 
resonant bar weight 2 

in 

71 

DYB 

Initial Y velocity of anti-resonant 
bar weight 2 

in/sec 

72 

ZB 

Initial Z displacement of anti- 
resonant bar weight 2 

in 

73 

DZB 

Initial Z velocity of anti-resonant 
bar weight 2 

in/sec 
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7.2.12 Modal Substructure MSI 

Modal Substructure MSI represents a dynamical system 
employing normal modes of a substructure. Each sub- 
structure of type MSI is limited to the use of one 
normal mode. The user must supply input defining the 
modal generalized mass, damping ratio, and frequency of 
the normal mode. These properties may derive from 
analysis or from test data. Processing in the progra;m 
forms damping and stiffness of the mode from the input 
data. Each substructure of type MSI is allowed to have 
up to five nodes. At each node, structure orientation 
angles and elements of the modal matrix for the normal 
mode are defined. At each such node, the elements of 
the modal matrix are three linear and three angular 
displacements for a unit value of the normal mode 
amplitude. The directions of the displacements at a 
node must be the same as the directions in Figure 13. 
The specification of the same node nvunber among two 
substructures of any type leads to the connecting of the 
two substructures at the node. More than one modal sub- 
structure of type MSI is connected to the same node when 
it is desired to represent the effects on a node of 
several normal modes. Nodes on MSI substructures may be 
used to define fuselage modes at up to five points. . If 
a physical substructure is not connected to a node at 
which responses are desired, these responses may be 
displayed by connecting a GFl substructure to the node 
of interest, with zero forcing coefficients. The GFl 
substructure contributes no matrices or forces to the 
system. It introduces physical displacements at the 
node which become members of the set of dependent system 
coordinates which are used to monitor the node response. 
Initial modal amplitudes and velocities may be specified 
to analyze transient system response. 

COMPONENT MSI INPUT DATA 

Location Fortran Input Item Description Units 

No . Name 


1 

ZETA 

Modal 

Damping ratio 

2 

MO 

Modal 

Generalized mass 

3 

OMEGA 

Modal 

Frequency 


nd 

lb-sec2/in 

hz 


4 NNODE Number of connection nodes allow- nd 

ing the substructure to be con- 
nected to other substructures 
(maximum of 5) 
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5-9 NODE Connection node numbers nd 

A representation of a normal mode 
structure with M (greater than one) 

. normal modes is achieved by con- 
necting M substructures of type MSI 
to the same connection node. The 
value of M can be any number which 
does not cause the system to be 
solved to exceed the working stor- 
age allocation in the Base Program 
(see section 9.7 for control of 
working storage size). 

10-39 GAMMA Mode shapes at connection nodes. in/in & 

rad/in 

For each node specified in locations 
5-9, define elements of the modal 
matrix at the node contributed by 
the normal mode. The 6 modal matrix 
elements are: U (longitudinal), V 

(lateral), W (vertical), THETAX 
(roll), THETAY (pitch), THETAZ 
(yaw). 

Ex. loc. 5-9 contain 2, 5, 6, 8, 9. 

Then, loc 10-15 has mode shape for node 2, 
loc 16-21 has mode shape for node 5, 

loc 22-27 has mode shape for node 6, 

loc 28-33 has mode shape for node 8, 

loc 34-39 has mode shape for node 9. 

The monitoring of the response at a 
point on a normal mode structure 
which is not connected to another 
physical structure may be achieved 
by connecting a GFl substructure to 
the point of interest. The forces 
in the GFl substructure input are 
equated to zero. The GFl substruc- 
ture has null mass damping, and 
stiffness elements and six dis- 
placements at the connection node. 

The GFl substructure does not affect 
the system response. The values of 
the displacements at the connection 
node are calculated from the assem- 
bled system solution coordinates and 
are displayed as line printer output 
or are saved for plotting. The user 
may constrain one or more of these 
displacements from being displayed 
by nulling the mode shape elements 
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at the node in the senses which are 
not of interest-e.g. , the rotational 
elements of the mode shape at the 
node would be equated to zero . 


40-54 EULER Euler angles input. 

For each node specified in loc 5-9, 
input in sequence THETA (pitch), PHI 
(roll), XSI (yaw). Ex. for modes 
specified above, then 

loc 40-42 has Euler angles for mode 2, 
loc 43-45 has Euler angles for mode 5, 
loc 46-48 has Euler angles for mode 6, 
loc 49-51 has Euler angles for mode 8, 
loc 52-54 has Euler angles for mode 9. 


deg 


55 


Not presently used 


56 

Q 

Initial 

modal 

amplitude 

in/in 

57 

DQ 

Initial 

modal 

velocity 

-1 

sec 


7.2.13 Parameter Variation Mode (PVl) 

The parameter variation mode, PVl, is used to identify a 
parameter which is to be changed in design studies and 
the range of values to be assumed by the parameter. 
This operational mode provides a convenient and general 
mode of selection of parameters to be varied among the 
substructures available in the Base Program and flagged 
in the input. Up to 10 parameters from any of the sub- 
structures may be varied simultaneously. The user 
specifies in the input to PVl the element number and 
loader input location number of a component to identify 
•the parameter to be varied. The three character string 
PVl triggers the input for PVl. If the PVl string is 
not loaded, the parameter variation mode is not invoked, 
and its input is not required. 

COMPONENT PVl INPUT DATA 

Location Fortran Input Item Description Units 

No . Name 


1 

2 


FIRSTV Starting value of parameter to be 
changed in design study 


( as appro- 
priate) 


FINALV Final value of parameter to be 
changed in design study 


(as appro- 
priate) 
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3 

NPTS 

Number of points of parameter 
variation 

nd 

4-23 

NEL 

Element number and corresponding 

nd 


location of parameter to be varied 
(up to 10 pairs can be loaded). 
Ex., conduct parametric studies of 
variations in fixed absorber weight 
(FAl) and forcing frequency in the 


range of 5 
steps of 5. 

to 

20 (slugs or Hz) in 

If FAl is 

element No. .3 and FRl is 

element No 

. 6, 

then PVl input is : 

Loc 1 

5. 

(Starting Value) 

Loc 2 

20. 

(End Value) 

Loc 3 

4. 

(No. of Cases) 

Loc 4 

3. 

(FAl Element No. ) 

Loc 5 

2. 

(Location of Weight) 

Loc 6 

6. 

(FRl Element No. ) 

Loc 7 

1. 

(Location of. Frequency) 


7.2.14 Rotor Aeroelastic Substructure RE2 

Rotor aeroelastic substructure RE2 is used to represent 
a multi-blade aeroelastic rotor in hover. The function 
of this component is to read rotor matrices transmitted 
to a file by rotor aeroelastic analysis E927 and to 
prepare the rotor data for assembly with other compon- 
ents. Standard application of RE2 is in conjunction 
with eigensolution mode EG2 . 

The use of component RE2 is preceded by stand-alone 
application of E927. Rotor aeroelastic analysis E927 is 
an external program in the program package yielding the 
mass, damping, and stiffness matrices of a multi-blade 
elastic rotor in hover. The analysis is based on the 
assumption of linear perturbations about a given rotor 
trim state. Blade coordinates comprise flapping, 
leading, pitching, coupled flatwise/edgewise bending 
modes, and uncoupled torsion modes. Five hub displace- 
ments are represented. These are shown in Figure 13, 
except that the hub is restrained from yawing. Blade 
coordinates are transformed to rotor modes. Rotor modes 
consist of collective and first cyclic sine and cosine 
functions of rotor speed. Use of rotor coordinates in 
the rotor equations, and the assumption of equal blades, 
equally disposed on the rotor disc, enables a constant 
coefficient equation set to be derived. Rotor coord- 
inate selection is made in the input specification to 
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E927. Rotor aeroelastic analysis E927 transmits to a 
file the number of coordinates, flags defining coord- 
inates selected, and corresponding mass, damping, and 
stiffness matrices. Matrices embody dynamic and aero- 
dynamic contributions- The three character string RE2 
in the Base Program input stream triggers reading of 
data and Base Program processing required to represent 
the aeroelastic rotor. File data are read by the Base 
Program. Numerical input data following the RE2 flag 
consist of hub connection node number and hub orienta- 
tion angles . Processing by the Base Program includes 
the addition of a dummy yaw coordinate. Zero mass, 
damping, and stiffness matrices are created for this 
. coordinate- This yields a six coordinate hub consistent 
with coordinate representations for other substructures . 
Six hub coordinates enable the rotor component to be 
coupled to other components without changing Base 
Program assembly logic. A valid physical model of rotor 
yaw motion is not achieved. Null rotor matrices are 
associated with this displacement. The user must 
eliminate the hub yaw coordinate from the system model . 
This may be accomplished with constraint component CNl . 
A valid model would occur also if fuselage normal modes 
in component MSI have no hub yaw displacement. In the 
present program version, only one rotor of type RE2 may 
be requested. 

AEROELASTIC ROTOR RE2 INPUT DATA 
Location Fortran 

No. Name Input Item Description Units 


1 

I EXEC 

E927 program execution control 
= 0. Read matrices from unit 10. 
= 1. Execute E927 and read 
matrices from unit 10 . This 
option is not available on CDC 
system. 

nd 

2 

NCN 

Rotor hub connection node 
number. 

nd 

3 

THETA 

Hub connection node Euler 
pitch angle (0y) 

deg 

4 

PHI 

Hub connection node Euler 
roll angle (0x) 

deg 

5 

XSI 

Hub connection node Euler 
yaw angle (0z) 

deg 
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Coordinate 
Output Label 

QFT1-QFT4 

THTE 

THTE 

BETA 

GAMA 

QTIS, QTIC 
QT2S, QT2C 
QT3S, QT3C 
QT4S, QT4C 
THRS, THRC 

THES, THEC 

BETS, BETC 

GAMS, GAMC 

X, Y, 2, THTX 
THTY, THTZ 

NOTES ; a) 

b) 

c) 

d) 


AEROELASTIC ROTOR RE2 COORDINATE LABELS 

Description 


Blade bending collective modal 
amplitude ^ ^ ^ 

Blade pitching^^^ 

Blade torsion^^^ 

/ c ) 

Flapping collective mode' ' 

Leading collective mode^*^^ 

Blade bending cyclic sine and cosine 
modal amplitudes^ 

Blade pitching sine and cosine cyclic 
amplitudes 

Blade torsion sine and cosine cyclic 
amplitudes 

Blade flapping sine and cosine cyclic 
/ 0 \ 

amplitudes' ' 

Blade leading sine and cosine 
amplitudes ^ ^ 

Hub displacements. See Figure 13 for 
positive directions 

Positive flatwise bending up and edgewise bending 
forward. 

Positive control nose down. 

Positive up. 

Positive forward. 
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7.2.15. Rotor Aeroelastic Substructure RE3 


Rotor aeroelastic substructure RE3 is used to represent 
a multi-blade aeroelastic rotor in forward flight or 
hover. The function of this component is to read the 
rotor hub excitation force and moment vector and the 
rotor impedance matrix transmitted to a file by rotor 
aeroelastic analysis G400. Component RE3 triggers the 
calculations required to assemble the G400 rotor force 
vector and rotor impedance with impedance elements 
contributed by other dynamical components. The use of 
RE3 is preceded by stand-alone application of G400. The 
coordinates represented in the RE3 substructure are six 
hub displacements shown in Figure 13. The RE3 rotor 
impedance substructure is consistent only with a steady 
state forced response solution mode, FRl, The frequency 
of force excitation specified in the input to FRl must 
be the same as that specified in the input to G400 to 
derive the rotor impedance matrix. This frequency will 
be an integer multiple of rotor speed and blade number 
(wf = nb ) . Component RE3 may be flagged to calculate 
Higher Harmonic Control (HHC) angles to minimize vibra- 
tions. Optimal control angles, blade stresses resulting 
from hub displacements, and HHC control inputs may be 
calculated. This is achieved by reading from a file the 
matrices transmitted from G400, deriving the effects on 
hub loads and blade stresses of HHC angle perturbations 
and hub displacements. When the HHC mode is invoked, it 
is necessary to employ a GFl substructure to define the 
node at which a weighted sum of response is to be 
minimized. Weight factors applied to the 6 nodal 
displacements are defined in GFl input. Weight factors 
applied to the HHC control angles to limit control 
excursions are defined in RE3 input. The basis of the 
HHC minimization method is described in Reference 1. 

Component RE3 may be used to introduce the impedance and 
load excitation of a substructure which is not a rotor. 
A sample application would be the representation of an 
empennage. Substructure RE3 has an option to permit a 
user to load in unit 12 a load excitation vector and 
impedance matrix. The data may derive from any suitable 
source. Flag IRDFIL (see below) would be set equal to 
one. When a rotor is not represented, the use of RE3 is 
not preceded by stand-alone application of G400. The 
HHC and stress flags (see below) must be turned off. 
Only one substructure of type RE3 may be requested. 
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AEROELASTIC ROTOR RE3 INPUT DATA 
Location Fortran 

No. Name Input Item Description Units 

1 IRDFIL Read option for rotor impedance nd 

matrix and force vector 

=0, read input from save file 11 

= I, read input from user file 12 

2 ISTFOR Rotor hub forces option nd 

= 0, leave force vector unchanged 
= 1, null force vector in the code. 

3 IDPFOR Display hub force option nd 

=0, display hub displacements 
= 1, display hub forces and moments 

Six hub translational and rotational 
displacements which are solutions to 
the assembled system equations are 
displayed if IDPFOR = 0. If this 
flag is 1.0, the interface reaction 
forces and moments between the hub 
and connected substructures - e.g., 
fuselage - consistent ’ with the 
motion of the assembled system are 
displayed (see Reference 1, Assembly 
of Rotor Impedance Matrix) 

There is an override instruction in 
the code which sets IDPFOR to zero 
if the user selects a display of 
system response in fractions of the 
acceleration of gravity (g's) 
(location 3 in FRl is equal to 1). 

In other words, the selection of a 
display of results in g's in FRl 
supersedes IDPFOR = 1 selection in 
RE3 and causes hub displacements to 
be displayed, even if hub forces and 
moments are requested. 

4 NON Rotor hub connection node nd 

5 THETA Hub connection node Euler deg 

pitch angel (0y) 

6 PHI Hub connection node Euler deg 

roll angle (0x) 
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7 

XSI 

Hub connection node Euler 
yaw angle (0z) 

deg 

8 

IHHC 

HHC control switches 
=0, do not activate option - 
default value 
= 1 , activate option 

nd 



This option is set to zero in- 
ternally if IRDFIL = 1 (loc. 1). 


9 

ITHHC 

Control input option for HHC 
applications 

= 0, compute optimal theta by 
quadratic minimization - no 
user input required - default 
value 

= 1, read HHC input angles as 
loaded in locations 11 through 
16 below 

nd 

10 

ISTRSS 

Calculate blade stresses option 

nd 



=0, do not calculate blade 
stresses. G400 location 933 
is zero. This option is set 
to zero internally if IRDFIL = 
1 (loc. 1). Default value is 
zero. 

= 1, calculate blade stresses. 
G400 location 933 is greater 
than zero. 




In the following, N refers to the 
ber of blades in locations 11 to 

nxim- 

19. 

11 

THNMIC 

HHC input angle - (N-l)th 
cosine component 

deg 

12 

THNMIS 

HHC input angle - (N-l)th 
sine component 

deg 

13 

THNC 

HHC input angle - Nth 
cosine component 

deg 

14 

THNS 

HHC input angle - Nth 
sine component 

deg 

15 

THNPIC 

HHC input angle - (N+l)th 
cosine component 

deg 
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16 


THNPIS 

HHC input angle - (N+l)th 
sine component 

deg 

17 


WTHNMl 

Weighting factor for (N-l)th 
HHC angle 

nd 

18 


WTHN 

Weighting factor for Nth 
HHC angle 

nd 

19 


WTHNPl 

Weighting factor for (N+l)th 
HHC angle 

nd 

7.2.16. 

Time 

History 

Solution Mode (THl) 



The time history solution mode THl yields the response 
of a system as a function of time. Transient responses, 
steady responses, and the stability of systems may be 
determined from the time history behavior of a system. 
The solution algorithm yielding the time history re- 
sponse is the Newmark Beta method described in Reference 
1. Section 8.3 describes the line printer output 
descriptors displayed by THl. 



COMPONENT THl INPUT DATA 


Location 

No. 


Fortran 

Name 

Input Item Description 

Units 

1 


BETN 

Newmark Beta factor. Range 
permitted is 0. to 0.25 

nd 

2 


BELT 

Time increment 

sec 

3 


TMAX 

Time limit for integration of 
the equations of motion. 

sec 

4 


OMEGA 

Forcing frequency 

hz 

5 


I DEBUG 

Debug printout switch 
= 0, for no printout - 
default value 
= 1 to display matrices 
for substructures and 
assembled system 

nd 

6 


IRSTRT 

Restart option 
= 0 for no restart case - 
default value 
= 1 for restart case. 

nd 




NOTE: Restart data are written 

to and read from file 8. 
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Section 8 


BASE PROGRAM OUTPUT DESCRIPTION 

Typical outputs from the Base Program are shown in Figures 17 
through 19. The outputs are presented in the same order as 
outputs obtained from program execution. 

Figure 17(a) is printed out as shown at the beginning of each run. 
It is the same for all runs; it is coded in subroutine INPUT 
(segment 3). Figures 17(b) and 17(c) are card image listings of 
the run titles and input data for the first case. If multiple 
cases are run, the new input data cards are printed at the be- 
ginning of each case. This sample run consists of one case only. 
It has 4 title cards with W (HZ) as the abscissa label to be used 
for plotting. Figure 17(c) shows that this case uses 5 com- 
ponents: FAl, GFl, THl, MSI, and GEN. 

8 . 1 Component Output 

Figures 17(d) through 17 (w) provide a narrative description of the 
input locations and the loader input values for each Base Program 
component. Component input narrative descriptions are listed in 
alphabetical order of the three character component names, start- 
ing with component BFl and ending with component THl. Figures 
17(d) through 17 (w) may be used to aid the preparation of com- 
ponent inputs. 

8 .2 Debugging Output 

After the component inputs are printed, a variety of debugging 
printouts are available to show the assembly of the substructure 
matrices, and intermediate calculations obtained for the different 
methods of solution. These printouts are obtained if the debug 
switch is set to 1.0 in any one of the components EGl, EG2, FRl, 
and THl. For all solution modes, siabroutine DSUB is called to 
identify the name of the subroutine providing the debug printout. 
Sxibroutine DMTRX is called to provide the debug output. Only 
non-zero quantities are displayed. A summary of the debugging 
printouts available for the four solution methods is presented in 
Table 5. Sample outputs are illustrated in Figures 18(a) through 
18(e). Figure 18(a) shows part of the assembly of the system 
mass, damping, and stiffness matrices as performed in subroutine 
ASMMCK before the real eigensolution EG2 is executed. Matrices 
MO, CO, and KO are svibstructure matrices referred to substructure 
local coordinates. Matrices M, C, and K are matrices for the 
assembled system obtained after mapping to independent coordin- 
ates. Figure 18(b) provides a sample output of the real eigen- 
solution EGl results obtained from subroutine S02EG1. In Figure 
18(c) results from the real eigensolution EG2 are shown. These 
are obtained in subroutine S02EG2 . Typical results obtained for 
the forced response solution FRl are illustrated in Figure 18(d) 
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in subroutine S02FR1. Finally, time history solution results are 
shown in Figure 18(e) in subroutine S02TH1. 

8.3 Solution Methods Output 

Figures 19(a) to 19(d) are samples of standard outputs obtained 
with the four solution methods EGl, EG2, FRl, and THl . For all 
solution modes, the display consists of seven fields of data 
containing alphanumeric name strings and numeric data keyed to the 
alphanumeric strings. These fields are listed by number below and 
contain the following data: 

1) Item counter for variables listed in each output row. 

2) The word ELEMENT. 

3) Element number of a substructure, listed in ascending 
order . 

4) Substructure three-character name or a blank string, or 
strings FREQ, DAMP, or TIME. 

5) Character name of a variable. Up to four characters may 
be used. If field 4 (above) is blank, the variable 
listed in each row is not a coordinate. The variable 
may be an eigenvalue or time. If field 4 contains a 
substructure name, the variable in the row is a coordi- 
nate . 

6) Solution mode description. 

a) For EGl, the descriptor is EIGV. This string is 
followed by the eigenvalue number, where the 
eigenvalue number is an integer ranging from one to 
the niunber of independent coordinates for the 
assembled system. (See Figure 19(a)). 

b) For EG2, the descriptor is EIGV when the string 
FREQ or DAMP occurs in field 5. The string FREQ is 
the descriptor for frequency of the eigenvalue in 
hz and the string DAMP is the descriptor for 
damping of the eigenvalue expressed as a fraction 
of critical damping (see equations (40) and (41) in 
Reference 1 ) . When the string in field 5 is the 
name of a coordinate, the string in field 6 is 
either REAL or IMAG, denoting the real and imagin- 
ary values of an element in the eigenvector. 
These strings are followed by the eigenvalue number 
(integer ranging from one to number of independent 
coordinates for the assembled system). The output 
displays are repeated until results for all eigen- 
values are displayed. (See Figure 19(b)). 
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,c) For FRl, the descriptors are AMPLITUD and PHASE 
denoting amplitude and phase of a coordinate 
response. If the display conversion factor ICONVG 
in FRl is equal to 0, amplitude measures are feet 
for translations, radians for rotations, and 
non-dimensional for modal amplitudes. Phase 
responses are expressed in degree measure. If a 
hub interface force and moment display is requested 
in component RE3 ( IDPFOR = 1 ) , forces are in lb 
units, and moments are in ft-lb units. Phase 
responses are in degree measure. If the display 
conversion factor ICONVG in FRl is equal to 1, 
displacement responses (translational, rotational, 
and modal) are converted to accelerations, divided 
by the acceleration of gravity, and displayed (see 
equations (28) and (29), Reference 1) as shown in 
Figure 19(c). 

The notation used in the calculations of the 
coordinate amplitude and phase angle is : 

X = AMPLITUD* cos(Wpt-PHASE) 

d) For THl, the descriptors are SEC and DISPMENT 
denoting time and displacement. Time is expressed 
in seconds and displacements are feet for transla- 
tions, radians for rotations, and non-dimensional 
for modal amplitudes. (See Figure 19(d)). 

7) The numerical value of the variable. 

The appropriate units are as follows : 


a) 

eigenvalue 

- 

frequency is in hz; 



- 

damping is non-dimensional 

b) 

eigenvector 

- 

translation is non-dimensional 



- 

rotation is 1/rad 

c) 

translation 

- 

ft, non-dimensional, or g*s 

d) 

rotation 

- 

rad, or (rad/sec2)/g 

e) 

hub reactions 

- 

forces are in lbs 



- 

moments are in ft-lbs 

f) 

phase 

- 

degrees 

g) 

time 

- 

sec 


After the solution results are printed, the computer program 
working storage requirement for the run is displayed in terms of 
the maximum nvimber of working storage words needed for the run as 
shown in Figure 19(e). The significance of this value is dis- 
cussed in the section 9.7, "Base Program Size". The solution 
results are sent to unit 2 to save data for the Tektronix plot 
package. The data are formatted to facilitate checking of the 
Base Program. Figure 19(f) was obtained by transmitting the plot 
file data to the line printer. These data are not displayed in 
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the course of a normal run performed to obtain plot results. The 
output data in Figure 19(f) are discussed below to aid the inter- 
pretation of the program by users. 

A sample output is presented in Figure 19(f) for a time history 
solution up to 0.05 seconds. The format is typical of all data to 
be plotted- The title cards, up to 4, are shown first. The first 
three title cards comprise the plot heading while the fourth title 
card defines the abscissa, which, for this example, is TIME (SEC). 
The next line provides 3 pieces of infojnnation; the total number 
of steps (or cases), (6.1), the total number of system coordi- 
nates, (15.1), and the parameter value for 3-dimensional plots 
(0.1200). Next, the variables are listed in order of increasing 
element number; in this case, element 1 is component FAl, element 
2 is component GFl, and element 4 is MSI. Element 3 is THl which 
is represented by the first variable, which is time in seconds. 
The data to be plotted follow the list of variables. The first 
card contains the step (or case) number (1.1 through 6.1) and the 
abscissa value (this is the parametric value if PVl component is 
activated or time if the time history solution is exercised). In 
this example, it is seen that the variable time has values from 
zero to 0.05 in steps of 0.01 second. The numerical values of the 
final system coordinates are then displayed. For this example, it 
is seen that the following plots could be generated versus time: 


1. 

variable 

2, 

FAl 

BELT 

2 . 

variable 

3, 

FAl 

X 

3 . 

variable 

9, 

GFl 

X 

4. 

variable 

15 

, MSI MODE 


(fixed absorber mass translation), 
(fixed absorber base translation), 
(generalized force connection node 
displacement). 

(normal mode amplitude). 


A plot of these four variables versus time is illustrated in 
Figure 19(g). Up to 5 curves can be shown on the same plot. 

When the time history solution is used, the final results are 
written to unit 8 in subroutine SV3TH1 for restart conditions. A 
sample output is presented in Figure 19(h) for time increments of 
0.01 up to 0.05 seconds. As shown, the displacement, velocity, 
acceleration, and force vectors are provided at each time step of 
the run. To restart a time history solution, only the data at an 
initial time are needed as input. These data are acquired from 
the last time step of a prior run. The results from Figures 19(f) 
and 19(h) originate from the same computer run. At a time of 0.05 
second, the displacements are the same, as indicated by the 
underlined quantities. 
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Section 9 


OVERALL PROGRAM STRUCTURE 

This section provides information on the overall program structure 
to. enable users and programmers to understand the program logic 
and the overlay segmentation structure. Flow diagrams for the 
different overlay segments are presented. Brief descriptions are 
provided for all the Fortran subroutines. The routines which 
require modification for CDC computer operation and the IMSL 
routines which must be replaced by the contractor are identified. 
The COMMON Blocks and input/output data files employed by the Base 
Program are described. Also, the determination of the computer 
storage requirements and of the program size are discussed as well 
as the method employed to change the computer program size. 
Finally, revisions to the computer program subroutines necessary 
to add a new substructure or to modify an existing substructure 
are presented in detail . 

9 . 1 Segmentation Structure 

The Base Program SIMVIB was overlaid to minimize the CDC computer 
storage requirements. The overlaid program size comprising a 
23 . 6 K 10 (decimal) word instruction bank and a 8.4Kio (decimal) 
data bank is 32Kio (decimal) words. This is equivalent to 100 Ks 
(octal) words. This limit complies with the requirement specified 
by the NASA Technical Monitor for interactive operation on the 
Langley Research Center CDC computer. The program size can be 
increased by specifying the size of the lA vector in the MAIN 
routine to enable the program to be applied to large problems (see 
section 9.7). 

The 18 overlay segments are shown in Figure 20. A brief descrip- 
tion of the function performed by each segment is provided in 
Table 5. The Fortran subroutines needed in the execution of the 
overlay segments are listed in Table 6. This table also includes 
a list of COMMON blocks. All COMMON blocks are stored in the main 
segment (number 1). IMSL routines which must be replaced by the 
contractor with equivalent routines are marked with an asterisk 
(*) while the routines which must be modified or replaced for CDC 
operation are indicated by a plus sign (+). 

9.2 Flow Diagrams 

Computer logic flow diagrams are presented in Figures 21 through 
41 for the high level subroutines. The flow chart in Figure 21 
illustrating the main program logic flow shows the parametric 
variation and the multiple cases loops. 
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9.3 Fortran Subroutines Descriptions 


To aid the understanding of the Base Program, an alphabetical list 
of Fortran subroutines is provided in Table 7; the appropriate 
segment number containing the routine is also provided in this 
table. The functions of the subroutines are given in Table 8. In 
addition, source code in the Base Program listing supplied to the 
NASA Structures Laboratory is commented to facilitate an under- 
standing of the program. An asterisk (*) denotes IMSL subroutines 
which have to be replaced by equivalent IMSL subroutines available 
at the Langley Research Center CDC facility. A plus sign ( + ) 
denotes system dependent subroutines. These routines contain CDC 
dedicated lines of code which are commented in column 1 . These 
comments must be removed and IBM dependent code corresponding to 
these lines in their immediate vicinity must be commented out. 


9.4 Fortran Subroutines Revisions for CDC Operation 

The revisions to the IBM Fortran subroutines which must be per- 
formed in order for the Base Program to execute on the CDC com- 
puter system fall into the three categories listed below: 

1. 3 routines must be replaced: ACONDP, ACON, and CORE. 

2 8 routines must be revised: FFLOAD, GVPLOT, INPUT, 

LBLEGl, LBLEG2, LOAD, MAINSV, and NAMES. 

The revisions are discussed in Table 8 and are also 
explained in the program Fortran code. 

3 . 18 IMSL routines must be replaced with equivalent IMSL 

routines by the contractor: EIGZF, EQZQF, EQZTF, EQZVF, 

LEQT2F, LINV2F, LINV3F, LUDATF, LUELMF, LUREFF, UERTST, 
UGETIO, VHSH2C, VHSH2R, VHSH3R, VXADD, VXMUL, and VXSTO. 

9.5 COMMON Blocks 

The Base Program utilizes a total of 18 COMMON blocks in its 
Fortran code. The Fortran subroutines in which the different 
COMMON blocks appear are presented in Table 9. Both the Fortran 
subroutines and the COMMON blocks are listed alphabetically for 
easy reference. The word size of a specific COMMON block is the 
same in all of the subroutines which reference it. 

9.6 Input/Output Data Transmission 

Input/output data transmission is accomplished in the Base Program 
through the use of nine files which are described in Table 10. 
The data files are assigned Fortran names in subroutine BLOCK, 
which is the first subroutine called in the main program. A 
description of the type of data transfer and the subroutines which 
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use a specific file are provided in the table. The input/output 
files are used to read input data (1, 5, 8, 10, 11 and 12), write 
output results (6), store solution results (2), and punch results 
on cards (7). The punch file is not presently used in the Base 
Program but can be readily activated. 

9 . 7 Control of Base Program Size 

The user of the Coupled Rotor/Fuselage Vibration Analysis may 
employ a convenient means of controlling the size of the Base 
Program. Size control may be achieved by changing two statements 
in the main routine. These statements as presently configured on 
the CDC system are: 



INTERACTIVE 

MODE 


BATCH MODE 

1. 

DIMENSION lA 

(8400) 

1. 

DIMENSION lA (40000) 

2. 

MAXA = 8400 


2. 

MAXA = 40000 


The first statement describes the size of the working storage 
vector (defined below) in terms of decimal words. The allocation 
of the working storage can be varied by changing the above state- 
ments. For example, 30000 decimal words may be allocated by 
modifying these statements to 

1. DIMENSION lA (30000) 

2. MAXA = 30000 

The program is recompiled and a new executable version created and 
stored for use. The new version will employ 30000 words of 
working storage. A user desiring to control the program size does 
not need to know any programming details beyond these statements . 

As a guide to programmers wishing to make other modifications, an 
introduction is given here of the space management capability 
contained in the Base Program. The working storage vector is a 
region of storage which is separate from regions used to store 
COMMON block data. It is a region of storage where arrays 
(literal strings, vectors of integers and real variables, ma- 
trices) may be placed or removed during program execution. Data 
are placed in the working vector lA or removed from lA by invoking 
in the program code two space management Fortran function subrou- 
tines. Function subroutine ISETSP allocates storage in lA for any 
array whose literal name, dimensions, and character (real, com- 
plex) are placed in the argument list of ISETSP. Function sub- 
routine IRMVSP may be used to remove the catalog entries (name, 
dimensions) applicable to the array referred to in ISETSP. 
Subroutine IRMVSP may be thought of as removing the space allo- 
cated in I A to that array. These procedures may be used to 
allocate or remove storage spaces allocated to arrays during the 
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execution of the program. In other words, they have the effect of 
controlling dynamically the amount of storage occupied in lA. The 
time history of variation of storage in lA with CPU time is shown 
schematically in Figure 42 for a hypothetical problem. The 
largest amount of storage taken up by arrays occupying lA is 
designated MAXSI2. The current storage in lA at any instant of 
time is ISIZE. The magnitude of ISIZE is always less than or 
equal to MAXSIZ. The program calculates ISIZE during execution 
and equates MAXSIZ to ISIZE if ISIZE exceeds a prior value of 
MAXSIZ. The program compares the value of MAXSIZ with MAXA. The 
parameter MAXA is equal to the size of lA. If MAXSIZ is less than 
MAXA, processing continues to completion. The program displays 
the value of MAXSIZ at the termination of the calculations for 
each case (see Figure 19(e)). If at any time in the calculations 
MAXSIZ is calculated to exceed MAXA, the processing terminates 
with a message informing the user to allocate more storage. 

If two or more cases are run successively, it is found that the 
value of MAXSIZ is larger for the second case as compared to the 
first case since computer storage allocated for the input com- 
ponents cannot be freed at the end of the calculations. As shown 
in the sample output in Figure 19(e), MAXSIZ is 1887 words after 
the first case and 2161 words after the second case - an increase 
of 274 words. Subsequent cases will have the same value of MAXSIZ 
as for the second case. 

Convenient size control provides a user with latitude to treat a 
wide range of problems. The version of the program used inter- 
actively on the CDC system employs 8400 decimal loads. A larger 
version of the program employs 40000 decimal words. This version 
can be used to solve large problems in a batch mode. For 
MAXA = 0, the program size is composed of storage required for 
instructions and common blocks, and is equal to 23600 decimal 
words on the CDC system. Thus, the version of the program with 
MAXA = 8400 would have a size equal to 23600 + 8400 = 32000 
decimal words. The larger program would have a size of 23600 + 
40000 = 63600 decimal words. 

Programmers must exercise care in the employment of function 
subroutines ISETSP and IRMVSP. A safe procedure for invoking 
these subroutines is to follow the applications in the program. 

A third space management routine is lASKSP which must be used by 
the programmer to return information on arrays which have already 
been allocated. Function subroutine lASKSP returns a pointer 
defining the position in the lA vector of an array named in its 
argument list. 

Programmers intending to add a new substructure to the code are 
not required to employ or understand the space management subrou- 
tines ISETSP, IRMVSP, and lASKSP. Section 9.8 describes the 
procedure for introducing code for a new svibstructure except a 
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substructure adding a forcing function such as GFl, or substruc- 
tures similar to RE2 or RE3 . The programmer must define in 
subroutine NSPACE the size (NDLOAD) of the input loader vector for 
the new substructure, and the size of a component dedicated work 
vector (NWORK). The value of NWORK is equal to 0 for all sub- 
structures (see Table 3) except RE2, RE3 , and THl . Substructures 
RE2 and RE3 involve the LINK subroutine (the RE3 substructure also 
includes logic to process impedance data) and the procedure for 
adding substructures of these types is different from the proce- 
dure for other substructures and is not described. The procedure 
for the addition of solution modes and operational modes capa- 
bilities is also not described. Modifications made to the program 
in these areas can be developed by following the code. 

The lA vector is equivalenced to a vector of real variables A. 
The programmer may recover integer variables by referring to the 
lA vector and real variables by referring to the A vector. The 
size of a word in the lA vector is the same as the size of a word 
in the A vector. Each subroutine called by the main routine 
MAINSV has an argument list containing lA and A as its first two 
members. The equivalence of lA and A is achieved by naming the 
first two arguments lA in the calls made to the subroutines in 
MAINSV. 

9.8 Addition of New Substructures and Modification of Existing 

Substructures 

The addition of a hypothetical substructure to the Base Program, 
referred to as BF4 for this example, requires two main steps: 

1. provide 7 new subroutines, as shown in Table 11. 

2. revise 18 existing subroutines, as shown in Table 12. 

The new subroutines to be developed for the BF4 component are 
listed in Table 11. They would be similar to those illustrated in 
Figures 43(a) through 43(g) for component BFl . The BFl routines 
are provided as a guideline only and therefore must be revised 
appropriately for each new substructure being added. 

The revisions to existing Base Program subroutines required to 
incorporate a new substructure are summarized in Table 12 . All 
the subroutines contain labelled COMMON CNAMES which is modified 
to include the new element BF4. In this illustration, assiutie that 
BF4 is the 21st element to be added to the Base Program. 

The procedure described above for adding code for a new svibstruc- 
ture applies to any substructure except substructures introducing 
a forcing function, like GFl or substructures introducing rotor 
effects, like RE2 and RE3 . The procedure described does not apply 
to solution or operational modes . 
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Modifications made to existing substructures involve the same set 
of routines listed in Tables 11 and 12. Many modifications would 
be limited to a subset of the routines listed in these two tables. 

Section 10 

TEST CASES DESCRIPTION AND RESULTS 


The Base. Program has been exercised with all program components. 
Check cases have been run to check out significant combinations of 
substructures and solution modes. A list of 16 test cases is 
provided in Table 13 . • The test cases include application of 

inplane and vertical bifilars, fixed system absorber, several 
isolator models (Cases 7 to 11 - see Figure 44), rotor contribu- 
tions from two different rotor models, and elastic beam represen- 
tation. The four methods of solution, which are eigensolutions 
with and without system damping, forced response, and time his- 
tory, are invoked and checked. Parametric variations are con- 
ducted for several design conditions. 

Restrictions affecting program applications are; 

1. to obtain non-zero responses, the GFl substructure must 
be used with the time history (THl) solution mode. 

2. when the time history (THl) solution is used, the 

parametric variation option (PVl) is not applicable. 

3. the only method of solution consistent with the rotor 
impedance representation (RE3) is the forced response 
solution (FRl). 

4. the GFl component must be connected to modal substruc- 
ture MSI at a point where it is desired to monitor 

responses and which is not connected to another struc- 
ture . 

5- only one substructure of type RE2 or RE3 may be request- 
ed at one time. However, the input may include both of 
these types. 

The program working storage requirement for each case is shown 
under the column heading "Size of Work Vector" in Table 13 . Most 

cases require less than 8000 words . Cases 11 and 12 have the 

largest storage requirement nearing 15000 words due to the utili- 
zation of four isolator units for case 11 and the complex eigen- 
solution application for case 12. Thus, all but cases 11 and 12 
can be run in the interactive mode. 


48 




From Table 13, the IBM central processing unit (CPU) time needed 
to execute a typical run is under 3 seconds . A time history 
solution of one second real time (case 5) requires 4.5 seconds of 
CPU time. This run also takes the longest to execute. 

The input and output for five representative cases (5, 7, 11, 12, 
and 13) are provided in Reference 11, Appendices A through E, 
respectively. They show results for different dynamic systems and 
methods of solution. 

Plots were developed for these cases except for case 12 using the 
Tektronix Plot Package described in Reference 7 . The plots are 
presented in Figures 45 through 47 for cases 5, 7, 11, and 13 
respectively. A summary of the plots generated and a brief 
description of each plot is presented in Table 14. Results for 
case 12 could be plotted if additional runs are made and all the 
results stored on unit 2. For example, the rotor blade collective 
pitch could be varied when executing the E927 program. The 
resulting matrices are then read from unit 10 by the Base Program. 
Finally, the final system eigenvalues and eigenvectors could be 
plotted against blade collective pitch. 

Up to 5 variables from one case may be plotted at the same time. 
Also, one variable from as many as 5 cases may be displayed on the 
same plot. 

Section 11 

IBM JOB CONTROL LANGUAGE (JCL) 

The IBM Job Control Language (JCL) needed to execute the Base 
Program or the IBM 370/168 computer system is shown in Figure 49. 
This setup reads basic input from unit 1, E927 matrices from unit 
10 if component RE2 is invoked, and G400 impedance matrix and 
force vector from unit 11 (or 12 if input by the user) if com- 
ponent RE3 is invoked. Time history restart conditions are either 
read from or written to unit 8 . Output results are sent to unit 
6. If plots are desired, the output results are stored on unit 2. 
Different options are presented in the JCL setup for writing of 
the output results to files on unit 2 for plotting or unit 8 for 
time history restart conditions. It should be noted that any JCL 
card starting with 3 or more asterisks is a comment card. Any 
card starting with // is a control card. The counter on the left 
hand side of Figure 49 shows that there are 13 control cards in 
the standard JCL setup required to execute the Base Program on the 
IBM system. 
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Section 12 


CPC PROGRAM INSTALLATION AND OPERATION 

This section provides the Job Control Language (JCL) for the 
installation and operation of the SIMVIB program on the NASA- 
Langley CPC NOS computer facility. The installation and operation 
procedures are based on protocol and system library facilities 
associated with the NASA-Langley CPC NOS computer facility. 

The JCL for installation of the SIMVIB program is presented in 
Figure 50 for the batch mode version { lA vector dimensioned 40000) 
and in Figure 51 for the interactive mode version (lA vector 
dimensioned 8400) . 

The JCL for execution of the SIMVIB program batch and interactive 
versions is shown in Figures 52 and 53 respectively. 

The JCL instructions apply to remote job entry operation (such as 
a telephone linkup). Revisions needed to run with card images are 
indicated in Figures 50 through 53. 
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Figure 1. Coupled Rotor/Airframe Vibration 
Analysis Data Flows. 











Coupled Rotor/Airframe Vibration Analysis Input/Output File Functions. 
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(a ) GENERAL INPUT CARD TORMAT 
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(b) SAMPLE INPUT CAROS FDR COMPONENT BFI 


Figure 4. Base Program Input Card Format and Examples. 


START 


^ * 














' 



BASE PROS:^AM - IHFUT SAMPLE X - X CASE 




ooox 

1 

FIXED 

SYSTEM ABS033ER - X 

MC3E SHAPE 




0002' 

1 ... 

FREC'JEKCY PARAMETRIC VASIATICM 




0003 

1 

FREQUE 

NCY (KZ) 






0004" 









0005 

1 

FAX 

X 




FIXED ABSORBER X 

*/ 

0006' 

3 . 

. X / 

X. 



/» 

COSNECTICM NODE NWBER 

*/ 

0007 

3 

2 / 

xo.o 



/» 

MASS (SLUGS) 


0008" 

3 

3 / 

0.X2 



/* 

DAMPING RATIO (N'O) 

»/ 

0009 

3 

4 / 

25000. 



/* 

SPRING STIFF.NESS (LB/FT) 

*/ 

0010- 

3 

5 / 

,1.0 

X.O..;,-, , 

.1.0 . 

/* 

INERTIAS (SLU3-FT*«2) 

*/ 

00X1 

3 

8 / 

0.0 

0.0 

0.0 

/* 

EULER ANGLES (DEG) 

»/ 

00X2:' 

1 

-XX / 

0.0 

0.0 

0.0 


INITIAL VALUES 

»/ 

00X3 









00X4' 

1 

GFX. 

2 



/« 

GENERALIZED FORCE I 

*/ 

00X5 

3 

X / 

X. 



/» 

COiNECTION NODE NUMBER 

*/ 

00X6' 

3 

2 / 

0.0 

0.0 


/* 

EULER ANGLES (DEG) 

*/ 

00X7 

3 

5 / 

1500. 

400.0 


/» 

X FORCE (LB) - COS , SINE 

*/ 

00X3' 

3. . 

. 7 / 

0.0 

0.0 



/» 

Y FORCE . - COS , SINE 

*/ 

00X9 

3 

9 / 

0.0 

0.0 


/* 

Z FORCE - COS , SINE 

»/ 

0020 

2 

XX / 

0.0 

0.0 


/* 

THETAX FORCE - COS , SINE 

*/ 

002X 

2 

X3 / 

0.0 

0.0 


/* 

THETA2 FORCE - COS , SINE 

«/ 

0022' 

1 .... 

_-X5 /. 

0,0 

0.0 



THETA3 FORCE - COS , SINE . 

*/ 

0023 









0024' 

1 

FRX 

3 



/* 

FORCED RESPONSE X 

»/ 

0025 

3 

X / 

X. 59X549 



/* 

FORCING FREQUENCY (HZ) 

*/ 

0026 

1 

-.2 /... 

X. 



./» 

GET. DEBUG PRINTOUT 

*/ 

0027 









0023' 

1 

MSX 

4 



/» 

MODAL STRUCTURE TYPE X 

■*/ 

0029 

3 

X / 

0.0 



/* 

DAMPING RATIO (ND) 

»/ 

0030' 

. 3 

2 / 

. X.O „ 



/» 

GENER. MASS ( LB-SEC»»2/IN) 

*/ 

003X 

3 

3 / 

0.0 



/* 

MODAL FRE0UE):CY (HZ) 

*/ 

0032 

3 

4 / 

X. 



/* 

NUMBER OF NODES 

*/ 

0033 

3 

5 / 

X. 



/* 

NODE NUr-SERS 

*/ 

0034 

2 

XO / 

...X.O 

0. 

3. 

/* 

MODE SHAPE 

*/ 

0035 

2 

X3 / 

0. 

0. 

0. 


MODE SHAPE 

*/ 

0036 

2 

40 / 

0. 

0. 

0. 

/'* 

EULER ANGLES (DEG) 

*/ 

0037 

1 

-56 / 

0. 

0. 

0. 

/* 

INITIAL VALUES 


0033 









0039 

X 

PVX 

5 



/» 

PARAMETRIC VARIATION I 

*/ 

0040 

3 

X / 

X. 59X549 



/* 

STARTING VALUE (HZ) 


004X 

3 

2 / 

9.549297 



/» 

FINAL VALUE (HZ) 

*/ 

0042 

3 

3 / 

6. 



/* 

NUMBER OF POINTS 

V/ 

0043 

3 

4 / 

3. 



/* 

ELEMENT NUMBER 

*/ 

0044 

1 

- 5 / 

X. 



/» 

LOADER LOCATION 

*/ 

0045 









0046 

X 

GEN 

6 




general element 

»/ 

0047 

3 

I / 

0. 



/* 

=0 NO PRINTOUT OF RESULTS 

»/ 

0040 

X 

-XO / 

0.X2 



/■* 

VALUE FCR 3-D PLOTS - DAM? 


0049 









0050 

X 

STOP 







0C5X 


Figure 5. Input Data Sample 1. 


56 



START 











* 









1 

BASE PROGRAM 

- INPUT SAMPLE 2 - 

3 CASES 



0001" 

1 

FIXED SYSTEM 

ABSORBER - 2 

MODE 

SHAPES 




0002 

1 

TIME HI 

STC.RY 

SDLUTIC.M . .. 






0003- 

1 

TIME tSEC) 







0004 










OOOS" 

1 

FAX 

X 





FIXED ABSORBER X 

*/ 

0006 

3 . 

. X /. 

X. , 




/» 

CCKTIECTICM NODE NUMBER 

*/ 

0007' 

3 

2 / 

xo.o 




/» 

MASS (SLUGS) . 

»/ 

OOOS 

3 

3 / 

0.X2 





DAMPING RATIO (NO) 

»/ 

0009- 

3 

4 / 

25COO. 




/* 

SPRING STIFFNESS (L3/FT) 

*/ 

0010 

3 .. 

. 5 / 

X.O.- 

x.o 

.. X.O 



INERTIAS (SLUG-FT*»2) 

»/ 

OOXl' 

3 

3 / 

0.0 

0.0 

0.0 


/* 

EULER ANGLES (DEG) 

*/ 

0012 

X 

-XX / 

0.0 

0.0 

0.0 



INITIAL VALUES 

»/ 

C0X3’ 

00X4 

, .1 

GFX 

2 





/» 

GENERALIZED FORCE 1 

•/ 

00X5' 

3 

X / 

X. 




/* 

CONNECTION NODE NUMBER 

*/ 

0016 

3 

2 / 

0.0 

0.0 



/* 

EULER ANGLES (DEG) 

*/ 

00X7" 

3 

.. 5 / 

X500. 

400.0 



/» 

X FORCE (LS) - COS , SINE 


0016 

.... _ 3 .. 

7 / 

0.0 

0.0 



./» 

Y FORCE - COS , SINE 


00X9" 

3 

9 / 

0.0 

0.0 



/* 

Z FORCE - CCS , SINE 

»/ 

0020 

Z 

XX / 

0.0 

0.0 



/* 

THETAX FORCE - CCS , SINE 

*/ 

0021' 

2 

X3 / 

0.0 

0.0 



/» 

THETA2 FORCE - COS , SlNE 

*/ 

0022 



--X5 

0.0 ... 

0.0 



/» 

THETA3 FORCE - CCS , SINE 

*/ 

0023' 

0024 

1 

THI 

3 




/* 

TIME HISTORY 1 

*/ 

0025' 

3 

I / 

0.25 




/* 

NENMARK BETA FACTOR 

*/ 

0026 

. . . 3 .. 

..... 2 .,/.... 

o.ox ... 




./» 

TIME I.NCREMENT (SEC) ..J 

*/ 

0027' 

3 

3 / 

.05 




/» 

MAX TIME (SEC) 

*/ 

0028 

3 

4 / 

6.0 




/* 

FORCING FREQUENCY (HZ) 


0029 

3 

5 / 

0.0 




/* 

DEBUG SELECTOR 

»/ 

0030 

...1 

- 6/ 

0.0 




/* 

RESTART FUG 

*/ 

0031 

0032 

I 

MSX 

4 




/* 

MODAL STRUCTURE TYPE 1 

*/ 

0033 

3 

X / 

0.0 




/» 

DAMPING RATIO (ND) 

*/ 

0034 

. 3 

2 / 

.X-0 .... 




/* 

GENER. MASS ( LS-SEC4*2/IN) 

*/ 

C03S 

3 

3 / 

0.0 




/* 

MODAL FREQUENCY (HZ) 

*/ 

0036 

3 

4 / 

X. 




/* 

NUMBER OF NODES 

*/ 

0037 

3 

5 / 

X. 




/* 

NODE NUMBERS 

*/ 

0033 

2 

10 / 

1.0 

0. 

0. 


/» 

MODE SHAPE ... 

*/ 

0039 

2 

X3 / 

0. 

0. 

0. 


/* 

MODE SHAPE 

»/ 

0040 

2 

40 / 

0. 

0. 

0. 


/* 

EULER ANGLES (DEG) 

*/ 

0041 

1 

-5S / 

0. 

0. 

0. 


/» 

INITIAL VALUES 


0042 










0043 

I 

MSX 

5 

COPY 

4 


/» 

MODAL STRUCTURE TYPE X 

*/ 

0044 

1 

-XO / 

0. 

x.o 

0. 


/» 

MODE SHAPE 

*/ 

0045 










0046 

1 

GEM 

6 




/* GENERAL ELEMENT 

*/ 

0047 

3 

X / 

0. 




/* 

= 0 f'O PRDfrOUT OF RESULTS 

*/ 

0048 

1 

-10 / 

0.12 




/* 

VALUE FOR 3-D PLOTS - DAMP 

*/ 

0049 










0050 

1 

EK'D 

X 

CASE 

X 


/* 

EK'D OF INPUT FOR CASE X 

■*/ 

0051 










0052 

X 

FAX 

1 




/* 

FIXED ABSORBER X 

*/ 

0033 

X 

- 2 / 

15.0 





MASS (SLUGS) 

»/ 

0054_ 










0055 

X 

END 

2 

CASE 

2 


/* 

END OF INPUT FOR CASE 2 

■*/ 

0056. 

0057 

1 

FAX 

i. 




/* 

Fixep AasORSEA 1 


0C53... 

0051 

■■ '1 

li. 

p 

b 

: ;■ 

- r 

r- ■ ■■ 


HASi _■ 


..DQfcO_ 

' f 




i 




OObl 


STOP 

• 

. 




, ' • 1 

5 

OA&Z. 

X 





* *■ 


. . 1 , 




Figure 6. Input Data Sample 2 



DEBUG PRINT SUBROUTINE ASItICK 


HO I !• 
HO ( I, 
HO 'I e; 
no c 3. 
no ( <1. 
no ( 5, 

HO- ( • s; 


no 

no r 6. 
no < ?• 
HO C n 
no ( fl. 

( a. 

( 9 , 

rr» 
no < 10. 
no ( 10 
no i 11 
HO •• I 11 r 
no ( 12, 
no ( 12. 
no ( 13. 
no ”( 13,“ 


no 

no 

H0“ 


( i<» 

( 19. 
I 15, 
(161 
f 16, 
( 16. 
no I 17. 
HO “r"I7V 


l|a 1 
lSls-3, 
41= 6. 
3)= 1. 
21 = 6 . 
71S-2. 

141= a: 
ai=- 2 . 

141= 3. 
71= 1. 
l4)«-9. 

ai= 1 . 

141«-2. 

7la-3. 

161=“ 3; 
101 = 1 . 
171S-3. 
111 = 1 , 
171=-i: 
121= 1. 
171= 5. 
61=-4. 
'131= 3; 
71=-4, 

121 = a. 

21= 3. 
'6l=-l. 
101=-1. 
161= 3. 
61= 3. 
"121= 2 i 


040490*02 

769500*01 

314640-02 * 

512540403 

314640*02 

096630-05 

004910*03^ 

096830*05 

696620*02 

746330*02 

104380-03* 

746330*02 

323380*01 

279930*01 

243510403* 

5125^*0403 

243510403 

516990403 

6S916D402 ~ 

516990403 

124120-01 

552190*03 

169640401' 

552190-03 

604940401 

657650*03 

5imo-o3 ■ 

556970402 

535790403 

241030*01 

562060*01“' 


HO 

1 1, 

2 |s-2. 096030*05 

no 

( 1 

no 

( 2, 

n=-2. 096630*05 

no 

1 2, 

HO 

1 ~i, 

“isli i'.S910tD-02 

ho ' 

( *3, 

HO 

I 3, 

413-2.045230*01 

no 

1 3, 

no 

( 4, 

313*2.045230*01 

no 

C 4< 

HO 

1 5, 

91= 5.236290-02 

no 

( S, 

HO 

• ( -5; 

"14)5-i;S347£B-'or 

TIO " 

r-£; 

no 

( 6. 

101= 5.236290-02 

no 

( 6, 

no 

1 ii 

16)3-3.023950-03 

no 

( 6, 

no 

1 7. 

9)3-3.279930-01 

HO 

( 7. 

Ho 

t 7rWl=-i.482650-»l’“ 

H,”f~7 

HO 

1 «■ 

10)3-3.279930-01 

no 

( a, 

MO 

1 «■ 

161« 6.051640-02 

no 

( a, 

no 

C 9. 

9)3 1.512540403 

HO 

( 9, 

HO 

r 9 , 

17)s-3. 113950402 

-Mr" 

TIB, 

HO 

( lOi 

12)3-2.045230-01 

no 

( 10, 

HO 

1 111 

5)3-6.007410-03 

no 

( 11, 

Ho 

1 111 

13)3-1.720990402 

no 

( 111 

no 

■*( 12, 

6)3-6:007410-03 

— MO • 

"f 12 

no 

1 12. 

131= 1.772650401 

no 

1 12i 

no 

1 13, 

5)3-1.949310-02 

no 

i 13. 

no 

r 13, 

10)3-5.276900400 

HO 

1 13, 

HO 

1 13, 

1/1= 2.1iU/6U401 

' MO ■ 

'ri4; 

no 

C 14. 

61S-1. 161690*01 

no 

< 14, 

MO 

1 14. 

141= 3.169640401 

no 

f 14, 

MO 

1 15, 

31= 6.614670401 

no 

C 15, 

MO 

1 “lO’; 

" 7)3-3.141030-01“ ■ ■ 

— HO" 

r 16; 

HO 

( 16, 

11)3*2.562060*01 

no 

f 16, 

no 

C 17, 

513-1.511960-03 

no 

( 17, 

no 

( 17, 

9)3-1.556970402 

no 

1 17, 

MO 

) 17, 

'13 S.11074DI0I 

— MO” 

TI7; 


. 31= 5 

. 21 = 1 
. ir= s 

. 151= 3 

> 41= 1 

. lll=-6 

ri7r»-5 
. 1213-6 
. 171= 1 
. Ill* 6 

. 121* 6 

> 171= 6 
I lll=-2 

; 

. 131=*1 
. 71= 6 

I 141= 1 

r'eii '6 
. 141= 1 
I 61= 4 
. llls-8 
r“5Is“4 
I 91=-5 

I 161=*2 

> 151= 3 

; 81S-3 

. 121=-9 
. 61= 6 
, 101=*1 
"171= “3 


.236290*02 
.746330*02 
:23a290-02 ' 
.445670402 
.516990403 
.007410*03 
623950-03^ 
607410*03 
336760*01 
.314640*02 
.631640-02" 
.314640*02 
,482050*01 
.045230*01 
:23d290-0r' 
.055360401 
.314840*02 
.772650401 
: 314840-02 ■ 
.720990402 
.002460*03 
.604940401 
:002460*03"“ 
.276900400 
.110760401 
,169840401 
.625820-02'' 
,295790401 
683620-02 
.621760403 
; 535790 J 03' 


CO ( 1, 11= 2.871030*01 CO I 

C0“ r l,“15I=-2. 058530400— — CO “I 
CO 12, 41= 4.058040400 CO 

CO 1 3j 31= 5.025660404 CO 

CO 1 4, 2)=*1. 206700400 CO 

CO ~ I '• 5;- 5 2 ;671030*01 CO 

CO I 5. 91= 2.337830401 CO 

CO ( 5, 131=-1. 637380400 CO 

CO ( 6. 51= S. 622330-01 CO 

CO- r li 91- 2.630540400 ~C0 *< 

CO ( 6, 131= 2.128160-01 CO I 

CO 1 7, 51= 1.167550-01 CO ( 

CO I 7, 91=-1. 013460401 CO ( 

C0 "-| - 7, *131=-4;23C670-01 C0“ (“ 

CO 1 8, 51=-1. 134120-03 CO I 

CO ( 6. 91=-1. 772330401 CO ( 

CO ( 8, 131= 1.543620-01 CO ( 

CO- ( 9, - 5)= 2.236930401 CO 


1. 2)3 1.351140-01 

" 2, "ns 1.167550-01 “ 

2. 15)3-3.502750400 

3. 41=*4. 165200403 

4. 3)3*6.720900403 

5. '"61s*5;622330*0I“ 
5, 10)3*2.^030540400 

5, 141= 2.126160-01 

6. 61= 2.071030-01 

or 10)3 2.337830401 - 

6, 141= 1.637360400 

7. 61= 1.134120-03 

7. 101= 1.772330401 

“7i 141= 1.543620*01-- 
6, 61= 1.167550*01 

6, lOls-l. 013460401 
6. 141= 4.236670*01 


C -9 ,^ 6)s-2:830540400 — ' 


CO ( 
CO ■"( 
CO ( 
CO I 
CO c 
"XO "f" 
CO ( 
CO ( 
CO I 
“CO I" 
CO ( 
CO I 
CO f 
* CO--t “ 
CO I 
CO f 
CO ( 


CO 

( 9. 

91= 5.025660404 

CO 1 

9, 10)3*6.173100404 

CO 

( 9, 

CO 

1 9i 

131= 1.586630403 
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Figure 7. E927 Aeroelastic Rotor Matrices. 
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Figure 7. Concluded. 
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Figure 8. Sample User Input of Rotor Impedance 
Matrix and Force Vector. 
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“ .12600*^ ' 

"^5062'7 

“^2697".’’ 

YAH 

HOMT 

(*, OMEGA DIR.) 

SIN 

-17.10 

.15370*06 

-.69610*05 

.35510*05 

9127. 

.19600*06 

-.67960*06 






VXPS 

WPS 

VZPS 

HXPS 

HYPS 

UZPS 

LOHG^. 

SHR 



_COS^ 



_5162. 

^280. 

_1.*31655*0.5_ 

^2262^05. -7672. 

LAT. 

SHR 

(*, STRBO) 

cos 


5665. 

.39560*05 

-1685. 

.20760*05 

-.16620*05 

6162. 

VERT. 

SHR 

(♦t UP) 

COS 


.10960*05 

9765. 

.23950*05 

-1066. 

-3633. 

-.21660*05 

ROLL 

HOMT 

(i, STRBO UP) 

cos 


-.16760*05 

6603. 

6639. 

-3193. 

.16150*05 

-7566. 

PIJCH 

_HOMT 

(♦» NOSE UP) 

__cos 


-3666^ 

_1566^ 

_^1M6. 

_:^.«9OT*05 

-6655_. ^ 

J2680*05 

YAH 

HOMT 

(«. OMEGA OIR.I 

COS 


-.16560*06 

.63050*05 

-.25330*05 

-.35730*05 

-.26710*06 

.75770*06 

LONG. 

SHR 

f«. AFT) 

SIN 


-.66600*05 

.65750*05 

1763. 

-1656. 

.30500*05 

.19160*05 

UT. 

SHR 

(»• STRBO) 

SIN 


-.66560*05 

-.75670*05 

-2601. 

7650. 

-.10260*05 

.52110*05 

• VERT. 

SHR 

UP)_ 

SIN 


^2001D*0£ 

_62M.^ 

-.2nio_*p^ 

_.1161D*05_ 


_J508._ _ 

ROLL 

HOMT 

STRBO UP) 

SIN 


-5066. 

-.20690*05 

2762. 

.26030*05 

-2737. 

-.39260*05 

PITCH 

HOMT 

i*, HOSE UP) 

SIH 


5396. 

313.0 

6602. 

-656.5 

.39290*05 

.36200*05 

YAH 

HOMT 

!♦» OMEGA DIR.l 

1 SIN 


-.79600*05 

-.11220*06 

5996. 

.16000*06 

.17350*06 

.12150*06 


Figure 9. Sample Output, of Rotor Impedance Matrix 
and Force Vector* 
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i 



N MASSES. ONLY ONE MASS SHOWN 

y =^y COSv^ + y SIN \j/ )(2/N) 
c s 


Figure 10. Horizontal Linear Bifilar BFl. 
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a 

t 



N MASSES. ONLY ONE MASS SHOWN 

/S =/Sq/N + (/3c COS«//+/3g SINv//)(2/N) 


Figure 11. Vertical Linear Bifilar BF2 . 
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NODE I 



Figure 12. Uniform Beam Segment BMl 



(EXTERNAL APPLIED FORCE) 



Figure 13 . Displacements and Forces Acting 

on a Connection Node of a Substructure. 
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OD 


TRANSMISSION 






NODE 2 



•MM 


0B IS ROTATION COORDINATE 


FOR 2ND ANTI- RESONANT BAR WEIGHT 
(W, NOT SHOWN) 


NODE 2 


SIDE VIEW 


TOP VIEW 


Two Degree-of-Freedom Vibration Isolator ISl. 


Figure 16 









« M 

M SIHVIB PROGRAM n 

M 

« „. . INPUT AND . IMTERMAL. CALCULATIOMS — M. 
» EMPLOY STANDARD UNITS * 

)( (FOOT* P0U}40» SECOt^D) ^ 




(a) Program Name and Units 


TITLE I - SIHVIB PROGRAM TEST RUN 
TITLE 2 - IHPLANE BIFILARS WITH 5 MOOES 
TITLE 3 - FORCED RESPONSE SOLUTION 
TITLE 4 - H (HZ) 


(b) Run Title Cards 


Figure 17. Base Program Output Samples 



..... I...N.-P..U.T D_E..CJ< C...A..SJ) 

KMKItKtOtMMMIHHtKtHtMtHHtMitItIHOHOliOiNliKKMKIOfItMMtttiMiHIKMKiHHtNKMKMNKIHHIlHHHHHIitMKMMNItKMMMWNIHtiHOiNNMMKi 


_PAL 


1 


- J. 



.QQfi60QJ) 

1 

/ 

i.odobtoo 

0.0 

0.0 

/* CONNECTION NODE NUMBER 

»/ 

0007000 

a 

/ 

l.OOODtOl 

0.0 

0.0 

/* MASS 

*/ 

0006000 

3 

/ 

1.2000-01 

0.0 

0.0 

/* DAMPING RATIO 

*/ 

0009000 


/ 

2.5000+04 

0.0 

0.0 

/* SPRING STIFFNESS 

*/ OOlOOOO 

5 

/ 

1.0000+00 

1.0000+00 1.0000+00 

/» MOMENTS OF INERTIA 

*/ 

6011000 

8 

/ 

0.0 

0.0 

0.0 

/» EULER ANGLES 

*/ 

0012000 

-H 

/ 

1.0000-02 

2.000O- 

02 0.0 

/* INITIAL VALUES 

*/ 

0013000 

"GFi 

— 

2 



_ . 



— *— 

0015000 ~ 

1 

/ 

1.0000+00 

0.0 

0.0 

/» CONNECTION NODE NUMBER 

*»/ 

0016000 

2 

/ 

0.0 

0.0 

0.0 

/» EULER ANGLES 

«/ 

0017000 

5 

/ 

1.5000+03 

4.0000+02 0.0 

/» X FORCE - COS t SINE 

*/ 0016000 

7 

/ 

0.0 

0.0 

0.0 

/* r FORCE - COS , SINE 

*/ 

0019000 

9 

/ 

0.0 

0.0 

0.0 

/» Z FORCE - COS , SINE 

*/ 

0020000 

11 

/ 

0.0 

0.0 

0.0 

/* THETAl FORCE - COS t SINE 

M/ 

0021000 

13 

/ 

0.0 

0.0 

0.0 

/» THETA2 FORCE - COS , SINE 

*/ 

0022000 

-15 

/ 

0.0 

0.0 

0.0 

/» THETA3 FORCE - COS > SINE 


0023000 

THl 


3 


0 



0025000 

1 

/ 

2.5000-01 

0.0 

0.0 

/* NEUMARK BETA FACTOR 

*/ 

0026000 

2 

/ 

1.0000-02 

0.0 

0.0 

/* TIME INCREMENT SEC 

*/ 

0027000 

3 

/ 

1.0000+00 

0.0 

0.0 

/» MAX TIME SEC 


0026000 

4 

/ 

6.0000+00 

0.0 

0.0 

/* FORCING FREQUENCY HZ 

»/ 

0029000 

5 

/ 

0.0 

0.0 

0.0 

/« DEBUG SELECTOR 

»/ 

0030000 

-6 

/ 

0.0 

0.0 

0.0 

/» RESTART FLAG 

»/ 

0031000 

nsi 


4 


0 

/« MODAL STRUCTURE TYPEl 

*/ 

0033000 

1 

/ 

0.0 

0.0 

0.0 

/» DAMPING RATIO 


0034000 

2 

/ 

1.0000+00 

0.0 

0.0 

/* GENER. MASS LB-SECKM2/IN 

«/ 

0035000 

3 / 

0.0 

0.0 

0.0 

/» MODAL FREQUENCY HZ 

*/ 

0036000 

4 

/ 

1.0000+00 

0.0 

0.0 

/* NUMBER OF NODES 

*/ 

0037000 

5 

/ 

1.0000+00 

0.0 

0.0 

/» NODE NUMBERS 

*/ 

0036000 

10 

/ 

1.0000+00 0.0 

0.0 

/« MODE SHAPE 

*/ 

0039000 

13 

/ 

0.0 

0.0 

0.0 

/« MODE SHAPE 

»/ 

0040000 

40 

/ 

0.0 

0.0 

0.0 

/» EULER ANGLES 

*/ 

0041000 

-56 

/ 

2.0000-02 

0.0 

0.0 

/» INITIAL VALUES 

*/ 

0042000 

GEN 


5 


0 

/* GENERAL ELEMENT 

*/ 

0044000 

1 

/ 

0.0000+00 

0.0 

0.0 

/» =0 ND PRINTOUT OF RESULTS 

*/ 

0045000 

-10 

/ 

1.2000-01 

0.0 

0.0 

/» VALUE FOR 3-D PLOTS- DAMP 

*/ 

0046000 


STOP 


(c) Input Data Listing 


Figure 17 


Continued 



COHPOHEHTSBIFILAR HORIZONTAL LINEAR BIFILAR »»« 

1 NCNl CONNECTION NODE NUMBER END 1 (NDl 

2 N TOTAL NUMBER OF BIFIURS (NO) 

3 MB BIFILAR MASS . (SLUGS) 

4 RR DISTANCE FROM CENTER OF BIFILAR TRACKING HOLE TO CENTER OF ROTATION (INCH) 

5 R _ .EQUIVALENT PENDULUM ARM (INCH) - 

6 ZETGAM INPLANE BIFIUR DAMPING RATIO C/CCRIT (NTD) 

7 OMEGA BIFIUR ROTATION SPEED (RAD/SEC) 

8 THETA EULER PITCH ANGLE AT NODE - ROTATE SECOND ABOUT THE Y-AXIS (DEGREES) 

? PHI,, EULER ROLL ANGLE AT NODE . - ROTATE THIRD ABOUT THE X-AXIS (DEGREES) 

10 XSI EULER YAH ANGLE AT NODE - ROTATE FIRST ABOUT THE Z-AXIS (DEGREES) 

11 GAMMAC COSINE COEFFICIENT OF INITIAL PENDULUM CYCLIC DISPUCEMENT (DEGREES) 

12 DGAMHC COSINE COEFFICIENT OF INITIAL PENDULUM CYCLIC VELOCITY (DEGREES/SEC) 

. 13,_. GAMMAS SINE COEFFICIENT OF INITIAL PENDULUM CYCCLIC DISPUCEMENT (DEGREES) 

14 DGAMMS SINE COEFFICIENT OF INITIAL PENDULUM CYCLIC VELOCITY (DEGREES/SEC) 


(d) Component BFl Data 
Figure 17. Continued. 


ELEMENT: 2 

1 

4 

9.31680D-01 

1.84200D401 

1. 960800+00 

2.500000-03 

2.70177D+01 

0.0 

0.0 - 

0.0 

0.0 

0.0 

0.0 

0.0 



«~J 

ro 


C0HP0HEKT:BIFILAR »»»*»»**» VERTICAL LINEAR BIFILAR *»* 


1 

2 

,...3,. 

4 

..-5,. 

6 

7 . 

8 
. ? 
10 
11 
12 
13.. 
14 
15. 
16 


NCNl . CONNECTION NODE NUMBER END 1 (ND) 

N TOTAL NUMBER OF BIFILARS (NO) 

MB BIFIUR MASS (SLUGS) 

RR DISTANCE FROM CENTER OF BIFIUR TRACKING HOLE TO CENTER OF ROTATION (INCH) 

R .EQUIVALENT PENDULUM ARM (INCH) 

ZETBET VERTICAL BIFILAR DAMPING RATIO C/CCRIT (NTD) 

OMEGA BIFILAR ROTATION SPEED (RAO/SEC) 

THETA EULER PITCH ANGLE AT NODE - ROTATE SECOND ABOUT THE Y-AXIS (DEGREES) 

PHI EULER ROLL ANGLE AT. NODE ...-.ROTATE THIRD ABOUT THE X-AXIS (DEGREES) - 

XSI EULER YAM ANGLE AT NODE - ROTATE FIRST ABOUT THE Z-AXIS (DEGREES) 

BETAO BIFILAR MASS COLLECTIVE INITIAL FUPPING ANGLE (DEGREES) 

DBETAO BIFILAR HASS COLLECTIVE INITIAL FUPPING VELOCITY (DEG/SEC) 

BETAC BIFIUR MASS COSINE COEFFICIENT OF INITIAL FUPPING ANGLE (DEGREES) 

DBETAC BIFILAR MASS COSINE COEFFICIENT OF INITIAL FUPPING VELOCITY (DEG/SEC) 

.BETAS.,.._ BIFILAR. MASS SINE COEFFICIENT OF INITIAL FUPPING ANGLE (DEGREES) 

DBETAS BIFILAR MASS SINE COEFFICIENT OF INITIAL FUPPING VELOCITY (DEG/SEC) 


(e) Component BF2 Data 


Figure 17. Continued. 


element: 

1 

4 

9.31660D-01 
1.84200D+01 
1.22800DtOO 
2.50000D-03 
• 2.70I77D*01 
0.0 

.. 0.0 ..^... 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 . — 

0.0 



C0MP0NEMT:BEAM1 


BEAM ELEHEMT TYPE 1 ««>«»««««*» 


ELEMENT 


1 NCNl . CONNECTION NODE NUMBER END 1 (NDl 1 

Z NCN2 COmECTION NODE NUMBER END 2 (ND) 2 

3 E . BEAM EUSTIC MODULUS (PSIJ - 3.00000D+07 

4 A BEAM CROSS-SECTIONAL AREA (FTJHf a 1 1.66700D-0I 

5 lYA BEAM SECOND AREA MOMENT OF INERTIA ABOUT THE Y-AXIS (FT»»^) S.03760D-06 

6 IZA BEAM SECOND AREA MOMENT OF INERTIA ABOUT THE Z-AXIS IFTftK'i) 3.21500D-05 

7 6J . BEAM TORSIONAL STIFFNESS (SLUG-FT/HAD ) 2.00000D+08 

8 MO HASS OF BEAM (SLUGS) 8,OOOOOD-02 

...S . IXX. BEAM SECOND MASS MOMENT OF INERTIA ABOUT THE X-AXIS (SLUG-FTiwE) 0.0 - 

10 lYY BEAM SECOND MASS MOMENT OF INERTIA ABOUT THE Y-AXIS (SLUG-FT«»2) 0.0 

11. IZZ ... BEAM SECOND MASS MOMENT OF. INERTIA ABOUT THE Z-AXIS (SLUG-FT»«»2) 0.0 

12 UCG OFFSET OF BEAM CENTER OF GRAVITY FROM EUSTIC AXIS IN X-DIRECTION (FT) 0.0 

,.13.. VCCi OFFSET OF BEAM. CENTER OF. GRAVITY FROM EUSTIC AXIS IN Y-DIRECTION (FT) 0.0 

1^ NCG OFFSET OF BEAM CENTER OF GRAVITY FROM EUSTIC AXIS IN Z-DIRECTION (FT) 0.0 

,15 lENGTH.. length . OF BEAM (FEET) - - — <i.00000D-01 

16 THETAl EULER PITCH ANGLE AT END 1 - ROTATE SECOND ABOUT THE Y-AXIS (DEGREES) 0.0 

17 PHIl EULER ROLL ANGLE ATEND 1 - ROTATE THIRD ABOUT THE X-AXIS (DEGREES) 0.0 

18 XSIl EULER YAM ANGLE AT END 1 - ROTATE FIRST ABOUT THE Z-AXIS (DEGREES) 0.0 

.19 ..theta? EUEER.PITCH ANGLE AT END 2...- ROTATE SECOND ABOUT THE Y-AXIS (DEGREES)- - 0.0 

20 PHI2 EULER ROLL ANGLE AT END 2 - ROTATE THIRD ABOUT THE X-AXIS (DEGREES) 0.0 

21_ XSI2 EULER .YA)( .ANGLE AT , END 2. - ROTATE, FIRST. ABOUT THE. Z-AXIS (DEGREES) 0.0 

(f) Component BMl Data 
Figure 17. Continued. 





COMPONENT:CONSTRI 


CONSTRAINT COMPONENT 1 »)hhh*)hhhh» 


ELEMENT 


1 NCN CONNECTION NODE NUMBER (ND) I 

Z CONSTR DEGREE OF FREEDOM SELECTOR 

= 0 DEGREE OF FREEDOM CONSTRAINED 

= 1. DEGREE OF FREEDOM . FREE 

X Y Z THTX THTY THTZ 
0 0 0 0 0 1 

8 THETA... .EULER PITCH ANGLE CDEGREES) - ROTATE SECOND ABOUT THE Y-AXIS 0.0 

9 PHI EULER ROLL ANGLE (DEGREES) - ROTATE THIRD ABOUT THE X-AXIS 0.0 

10 XSI EULER YAM ANGLE (DEGREES) - ROTATE FIRST ABOtn" THE Z-AXI5 0.0 

(g) Component CNl Data ‘ 

COMPONENT :EIGEN1 »»»»»»«»«» REAL EIGENSOLUTION ELEMENT: 

1 IDEBUG DEBUG SELECTOR 1 

= 0 ==> NO DEBUG PRINTOUT 

= 1 ==> TRACE MATRIX ASSEMBLY AND SOLUTION 

(h) Component EGl Data 

component: EIGEN2 »hhhhmhm(x» COMPLEX EIGENSOLUTION tiKitxitxiUMm ELEMENT: 


1 IDEBUG DEBUG SELECTOR 1 

..... __ ._ .y. ofBOG-pRlNTOUr" " ' 

= I ==> TRACE MATRIX ASSEMBLY AND SOLUTION 

(i) Component EG2 Data 


Figure 17. Continued 


COMPONENT JFIXABSl i m)nn nnnnn( FIXED SYSTEM ABSORBE R T YPE 1 ELEMENT? 


1. 

2 

3 

.. 

MO 

ZETA 

_,.CQNNECTJPN .NODE, HUMBERT mi. .. 

ABSORBER MASS (SLUGS) 

ABSORBER DAMPING RATIO (NO) 

l.OOOOOD+01 

1.20000D-01 


KO 

ABSORBER STIFFNESS (LB/FT) 

2.50000D+0<i 

.5 

.ixx_ 

ABSP.RPER. S_E.CQNP .MASS,. MOMENT .OF JNERTIA. ABOUT. THE X. AXIS.J SLUGrFT.««2.) . 

,.._1.000DOD*00. 

6 

lYY 

ABSORBER SECOND MASS MOMENT OF INERTIA ABOUT THE Y AXIS (SLUG-FT^^E ) 

l.OOOOOD+00 

7 

IZZ 

ABSORBER SECOND MASS MOMENT OF INERTIA ABOUT THE Z AXIS (SLUG-FT««21 

X.OOOQOD.tilQ. 

6 

THETA 

EULER PITCH ANGLE (DEGREES) - ROTATE SECOND ABOUT THE Y-AXIS 

0.0 

9 , 

_ PHI 

EULER ROLL ANGLE (DEGREES) - ROTATE THIRD ABOUT THE X-AXIS 

... 0.0 _ 

10 

XSI 

EULER YAM ANGLE (DEGREES) - ROTATE FIRST ABOUT THE Z-AXIS 

0.0 

11 

DELTA 

ABSORBER MASS INITIAL DISPLACEMENT (FT) 

l.Q0.0ODD.r02_ 

12 

DOELTA 

ABSORBER MASS INITIAL VELOCITY (FT/SEC) 

2.000000-02 


(j) Component FAl Data 


component: FORCER l mihcnkjhow* FORCED RESPONSE SOLUTION TYPE 1 ELEMENT: 

1 OMEGA FORCING FREQUENCY (HERTZ) 1. 952000+01 


2 IDEBUG DEBUG SELECTOR 0 

s 0 ==> NO DEBUG PRINTOUT 

, = 1 ==> TRACE MATRIX ASSEMBLY AND SOLUTION 

3 ICONVG OUTPUT DISPUY SELECTOR 0 

= 0 ==> OISPLACEMENTS (FEET) 

:= 1 ==> ACCELERATIONS (G'S) . 

(k) Component FRl Data 


Figure 17. Continued 



<Tl 


COMPPKEHTIPENIHPUT - ...t****!**"*! SEHERAL ..INPUT FOR . PROGRAM CONTROL ELEMENT; 7. 


..1 „,JCHTL1 PRINT SELECTOR . FOR FINAL. RESULTS 0 

so ==> SUPPRESS LINE PRINTER OUTPUT 

= I ==> FULL LINE PRINTER OUTPUT 

2_,....ICMr.L2._ERINT_.SELECTOR. FOR CQMPONENT.INPUTS — - - - - 1 , 

= 0 ==> SUPPRESS LINE PRINTER OUTPUT 

s 1 ss> full LINE PRINTER OUTPUT 

_3-9 =-sr..-_..0PEH..L0CATI0N3..F0tt FUTURE USE— ..... — ..._ - ...._ _ — — - - 

10 XINDEP INDEPENDENT VARIABLE FOR 3-D PLOTS 0.0 


(1) Component Gen Data 


Figure 17. Continued. 



CpMPPNENTjGENFCRl 


«»»»»**»»«.» GENERALIZED FORCE TYPE I t USED WITH FORCERI) 


element; 


.1 JMCN CONNECTIOH .NODE ,.NUMBER-tND)„. ..... i 

2 THETA EULER PITCH ANGLE (DEGREES) - ROTATE SECOND ABDUT THE Y-AXIS 0.0 

_EULER.RQLL_ANGLE..I DEGREES) _r..JJDTA.TE_THIRD ABOUT....TOE.X.-AXIS - — 0.0 — - 

4 XSI EULER YAH ANGLE (DE^EES) - ROTATE FIRST ABOUT THE Z-AXIS 0.0 

„5 FXCOS ... .COSINE...CQMPONENT. OF. X. DIRECTION. FORCE .. ILB)_. _ ... - _ 1.70000D+03 

6 FXSIN SINE COMPONENT OF X DIRECTION FORCE (LB) 0.0 

2_...F.YCOS _COSINE.^COMPONEMT_OF...Y. DIRECTION FORCE. (LB) — 0.0- 

8 FYSIN SINE COMPONENT OF Y DIRECTION FORCE (LB) l.lOOOOD+03 

_ ...9 .EZCOS„.. COSINE .COMEONENT-.OE...Z. DIRECTION. FORCE ILB). o.O — 

10 FZSIN SINE COMPONENT OF Z DIRECTION FORCE (U) 0.0 

JH....EI1C0S-..C0SINE.. COMPONENT. Ot.THETA...l. MOMENT— — .tlB).. — — 0.0 - 

12 FTISIN SINE COMPONENT OF THETA 1 MOMENT (IN-LB) 0.0 

13__FT2CQS .JCOSINE. COMPONENT. OF. THETA 2 MOMENT.. ._ (INrLB) ..... — 0.0- 

14 FT2SIN SINE COMPONENT OF THETA 2 MOMENT (IN- LB) 0.0 

.: i5_...EI3C0S— COSINE. CDMPONENI...DE THETA .3 MOMENT —tINrLB) - — o.O 

16 FT3S1N SINE COMPONENT OF THETA 3 MOMENT (IN- LB) 0.0 

_ 17 JfHRESP.. HHC..EU6 ..„ _ - - 0-. 

= 0 HHC NOT ACTIVE 

= 1 HHC ACTIVE 


.IS 

HZX 

WEIGHT. FOR X . RESPONSE 

19 

HZY 

HEIGHT FOR Y RESPONSE 

,20_ 

_... 

HEIGHT FOR Z RESPONSE 

21 

UZHX 

HEIGHT FOR THETAl RESPONSE 

22 

HZUY 

HEIGHT FOR THETA2 RESPONSE 

23 

UZUZ 

HEIGHT FOR THETA3 RESPONSE 


(m) Component GFl Data 
Figure 17. Continued. 


0.0 


0.0 


.. 0.0 

0.0 


._p,0 

0.0 



COnPOHCKnlSOUTEl ISOLATOB typ e 1 «»»iiiiin««im UEmHU 1. 


_ _l_ 

-•SHL. 

, _F(mcjlQH .HPP£- -L 

il« L _ 

1 

z 

NCN2 

CONttECrZOH NODE NUrtBER 2 

IND) 

5 

3 

AL 

DISTANCE BETWEEN NODE 1 A PIVOT 2 

IIH) 

1.250000*00 

4 

BL 

DISTANCE BETWEEN PIVOT 2 A Y-Z ZSQUTOR 

IIN) 

l.OOOOODtOI 

i X!-_ 

OISTAMCe BtTWEH PIVOT t t HOOT t 

I«U_ 


6 

OL 

DISTANCE BETWEEN PIVOT 2 A X ISQUTCR 

(IN) 

0.0 

- 2. 

Ml 

WEIGHT OF Y'Z ISOLATOR 

UBS» 

4.^«OOCO*«1 

a 

XIIY 

INUnA OF T-Z ZSOUTOa ABOUT THE T-AXIS 

fXH-LB-SEC«*2) 

0.0 

^ 

_xn^ 

INERTIA Of Y-Z ISOLATOR ABOUT THE Z-AXXS 

IZN-LB-SEC>i«2) 

0.0 

10 

M2 

WEIGHT OF X ZSOUTOR 

(LBS) 

0.0 

11 

XI2Y 

INERTIA OF X ISQUTOR ABOUT THE Y-AXIS 

IIN-LB-SECMH2) 

0.0 

It 

TKTX 

X STIFFNESS BETWEEN NODE 1 A PIVOT X 

ILfiS/IN) 

1.00000D«01 

_n 

JTCTY 

Y STIFFNESS BETWEEN NODE 1 A PIVOT 1 

(LfiS/IN) 



14 

TKTZ 

Z SnFFHESS BETWEEN NOOE 1 A PIVOT 1 

(LBS/JN) 

1.00000D«09 

IS 

TK2X 

X STIFFNESS BETWEEN PIVOT 1 A NOOE 2 

1 LBS/IN) 

1. 000000*05 

u 

TKEY 

T STIFFNESS BETWEEN PIVOT I A NOOE A 

<IBS/INI 

1.000000*05 


_mz_ 

Z STIFFNESS BETWEEN PIVOT 1 A NODE 2 

•LBSnNJ 

1.000000*05 

la 

TKAX 

X STIFfNESS BETWEEN NOOE 2 A PIVOT 2 

(LBS/IH) 

1.000000*09 

19 

7KAY 

Y STIFFNESS BETWEEN NOOE 2 A PIVOT 2 

ILBS/IN) 

1.000000*09 

to 

TKAZ 

Z STIFFNESS BETWEEN NODE 2 A PIVOT 2 

(LBS/IN) 

1.500000*04 


RKTO 

ROTATIONAL STIFFNESS OF NODE 2 ABOUT THE Y-AXIS 

C IN- LB/RAO) 

0.0 

tt 

RKTHT 

ROTATIONAL STIFFNESS OF PIVOT 2 ABOUT THE Y-AXIS 

( IN-LB/RAO 1 

0.0 

23 

RKPHT 

ROTATIONAL STIFFNESS OF PIVOT 2 ABOUT THE Z-AXIS 

fIN-LD/RAO) 

0.0 

24 

TCTX 

X OANPINS BETWEEN HXE 1 A PIVOT 1 

(HO) 

5.000000-02 


.lETL 

Y OANPING BETWEEN NODE I A PIVOT I 

(ro) 

5.000000-02 

26 

TCTZ 

Z OANPING BETWEEN NXE 1 A PIVOT 1 

()f)) 

5.000000-02 

27 

TCEX 

X OANPING between PIVOT 1 A NODE 2 

no) 

5.000000-02 

ta 

TC2Y 

Y OANPING BETWEEN PIVOT 1 A NOOE 2 

(HD) 

5.000000-02 


(n) Component ISl 


Figure 17. Continued 



29 

TCEZ 

Z OAnPIHG BETWEEN RIVOT 1 1 NODE 2 

(NDI 

S.BeSSSBrO* 

30 

TCAX 

X OAKPING BETUEEN NXE 2 ft PIVOT £ 

CNDI 

5.000000-02 


.T£*t_ 

_y .EAtjpm esmifttiiQSf 

J1P.L _ 


32 

TCAZ 

Z DAHPING BETWEEN NOOE 2 ft PIVOJ 2 

(NOt 

5.00000D-02 

33 

RCTHB 

ROTATIONAL OAKPING OF NODE 2 ABOUT THE Y-AXIS 

(NOl 

52QQQQ0DtO2_ 

34 

RCTHT 

ROTATIONAL DAMPING OF PIVOT 2 ABOUT TH Y-AXIS 

(ID) 

5.000000-02 

_ a. iPQu. 

RomioKAi o*npiHs Of PIVOT 2 about hl . t-axib 

inu _ 

— UQMOJtPi. 

36 

UTX 

X FR20UENCT BETMEEH NOOE 1 2 PIVOT 1 

IHZ) 

2.4S0O0D.01 

37 

KTY 

Y FREQUENCY BETWEEN NODE 1 ft PIVOT 1 

(HZI 

2.4eoooPFOi 

3ft 

KTZ 

Z FREQUENCr BETUEEN NODE 1 2 PIVOT 1 

fHZ) 

2.48000Dt01 

— ??_ 


X FREQUENCY BETUEEN PIVOT 1 ft NOOE 2 

JHZJ 


40 

U£Y 

Y FREQUENCY BETWEEN PIVOT 1 ft NOOE 2 

CHZl 

2.460000*01 

41 

UEZ 

Z FREQUENCY BETUEEN PIVOT 1 ft NOOE 2 

fHZ> 

2.460000*01 

42 

MAX 

X FREQUENCY BETUEEN NOOE 2 ft PIVOT 2 

(HZ) 

2.400000*01 


i**I_ - 

Y FREQUENCY BETUEEN NOOE 2 ft PIVOT 2 

JHZ) _ 

_2j^ft000^0i_ . 

44 

HAZ 

Z FREQUENCY BETWEEN NODE 2 ft PIVOT 2 

(HZ) 

2.420000*01 

45 



(HZ) 

0.0 

46 

WtHT 

ROTATIONAL FREQUENCY OF PIVOT 2 ABOUT TH Y-AXIS 

(HZ) 

0.0 


w»*nnr 

ROTATIONAL FREQUENCY OF PIVOT £ ABOUT TH Z-AXIS 

iHZJ 

0.0 

4ft 

13020 

CCKIROL SUITCH 2 FOR 2-0 DAVI 2 »> I FOR S-D DAVI ISOLATOR 

(ND) 

2 

49 

IFLEX 

ACCOUNT FOR FLEXIBILITY OF ISOUTOR BARS • NO ft ss>x . YES 

(ND) 

0 

50 

RXYHl 

FREQUENCY RATIO OF Y-Z ISOLATOR IN X-Y PLANE 

(ND) 

0.0 

_®L 

£XM_ 

FREQUENCY RATIO OF Y-Z ISOLATOR IN X-Z PLANE 

(ND) 

_0^ 

52 

RXZM2 

FREQUENCY RATIO OF X ISOUTOR IN X-Z PLANE 

(»)) 

0.0 

53 

DXYfll 

CRITICAL DAMPING OF Y-Z ISOUTOR IN X-Y PLANE 

(ND) 

0.0 

54 

oxzni 

CRITICAL DAMPING OF Y-Z ISOUTOR IN X-Z PLAtlE 

(ND) 

0.0 

_®L. 

oxai2_ 

CRITICAL DAMPING OF X ISOLATOR IN X-Z PLANE 

(W) 

0.0 

56 

THETAl 

EULER PITCH ANGLE AT END 1 - ROTATE SECOND ABOUT THE Y-AXIS 

(DEG) 

0.0 

57 

PHIl 

EULER ROLL ANGLE AT END 1 - ROTATE THIRD ABOUT THE X-AXIS 

(OE6) 

0.0 

5ft 

XSll 

EULER YAH ANGLE AT EFO 1 - ROTATE. FIRST ABOUT THE Z-AXZS 

(DEG) 

0.0 


(n) Component ISl Data 




Figure 17. Continued 




59 

THETA2 

EULER 

PITCH ANGLE AT END 2 - ROTATE SECOND ABOUT THE Y-AXIS 

(DEG) 

0.0 

60 

PHI2 

EULER 

ROLL ANGLE AT END 2 - ROTATE THIRD ABOUT THE X-AXIS 

(DEG) 

0.0 

_61_ 

XSI2 

-EULER 

JAM ANGLJ.AI.EW 2_-_ROJAlTI-FlRSI.ABpULTHS Z:AXIS„,., _ 

(DEG.) 

_fl.O 

62 

XT 

PIVOT 

1 

INITIAL X DISPUCEMEHT 

( INCH ) 

0.0 

63 

DXT 

PIVOT 

1 

INITIAL X VELOCITY 

(IN/SEC) 

0.0 

64 

YT 

PIVOT 

1 

INITIAL Y DISPLACEMENT 

(INCH) 

0.0 

_45. 

_DYT 

piyoL 

1_ 

INITIAL Y VELOCITY 

JJN/.§E.£} _ 

_ILJ . 

66 

ZT 

PIVOT 

1 

INITIAL Z DISPUCEMEKT 

(INCH) 

0.0 

67 

DZT 

PIVOT 

1 

INITIAL Z VELOCITY 

(IN/SEC) 

0.0 

68 

XB 

PIVOT 

2 

INITIAL X DISPUCEMENT 

(INCH) 

0.0 

_69_ 

DJ<B 

fiyOT_ 

Z 

INITIAL X VELOCITY 

JIN/SECJ 

_p._^ 

70 

YB 

PIVOT 

2 

INITIAL Y DISPUCEMENT 

(INCH) 

0.0 

71 

DYB 

PIVOT 

2 

INITIAL Y VELOCITY 

(IN/SEC) 

0.0 

72 

ZB 

PIVOT 

2 

INITIAL Z DISPLACEMENT 

(INCH) 

0.0 

.7L 

^Z^ _ 

piyoL 

2_ 

INITIAL Z VELOCITY 

iJN^EC) 

_o.q , 


(n) Component ISl 
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COMPONEMT:MODSTR1 


MODAL STRUCTURE TYPE 1 ****»)hhh(k 


element: 


.1 2ETA DAMPING RATIO (ND) . - 

2 MO GENERALIZED HASS (LB-SEC»»2/IN) 

. 3 OMEGA... MODE FREQUENCY (HERTZ) ... . 

4 NNODE NUMBER OF NODES DESCRIBED BY THIS MODE 

5 . NODE CONNECTION NODE NUMBERS OF NODES DESCRIBED BY THIS MODE 

1 0 0 0 0 

10 GAMMA MODE SHAPE. ENTER U> V, U. THETAX, THETAY. THETAZ FOR EACH NODE: 


3.40000D-02 

1.80000D+00 

1.61000D+01 

1 


55 

56 R 

57 DQ 


EMPTY LOCATION 

INITIAL. .NODAL AMPLITUDE (IN/INI.. 

INITIAL MODAL VELOCITY (1/SEC) 


(o) Component MSI Data. 
Figure 17. Continued. 




U 

. V 

U THETAX 

THETA Y 

T1 

10 - 15 

NODE 1 

1.0000D400 

1. 00000-02 

-2.40000-01 0.0 

0.0 

0.0 

16 - 21 

NODE 2 

0.0 

0.0 

0.0 0.0 

0.0 

0.0 

22 - 27 

NODE 3 

0.0 

0.0 

0.0 0.0 

0.0 

0.0 

28 - 33 

. NODE 4 

0.0 „. 

0.0 

. .. 0.0 0.0 

0.0 

0.0 

34 - 39 

NODE 5 

0.0 

0.0 

0.0 0.0 

0.0 

0.0 

40 EULER EULER ANGLES AT CO)(NECTION NODES. 

enter: 




THETA - 

EULER PITCH ANGLE. ROTATE SECOND ABOUT THE Y-AXIS (DEGREES) 



PHI 

EULER ROLL 

ANGLE. ROTATE 

third about 

THE X-AXIS (DEGREES) 



XSI - 

EULER YAM 

ANGLE. ROTATE 

FIRST ABO(n’ 

THE Z-AXIS (DEGREES) 





THETA 

phi 

- XSI 



40 - 42 

NODE 1... 

0.0 

..0.0 . 

0.0 



43 - 45 

NODE 2 

0.0 

0.0 

0.0 



46 - 48 

NODE 3 

0.0 

0.0 

0.0 



49 - 51 

NODE 4 

0.0 

0.0 

0.0 



52 -. 54, 

. NODE, 5 

0.0 

.0.0 

0.0 „.. 

.. 



0.0 

0.0 



00 

f\i 


component; PARMV l 


PARAMETRIC VARIATION TYPE 1 »»»»»»»)*»)» 


ELEMENT 


10 


1 FIRSTV STARTING VALUE FOR PARAMETRIC VARIATION 1.000000*01 

2 FINALV FINAL VALUE FOR PARAMETRIC VARIATION 3.000000*01 

3 , NPTS NUMBER OF POINTS IN PARAMETRIC VARIATION 10 

4 NEL GLOBAL ELEMENT NUMBER AND CORRESPONDING LOADER LOCATION FOR INDEPENDENT VARIABLE 
TO BE PARAMETRICALLY VARIED (UP TO 10 PAIRS) 

3 3 ... 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

- . ■ . ^ ^ 

(p) Component PVl Data 


COMPONEN T ; ROTORE L2 oimiminnm EU STIC ROTOR TYPE 2 witiHHtinnnnt .EiEMENT? Jl 

„ U.XEC _E927 JX^UTJONJCONTRO!. JF4.A5 0 

= b READ ROTOR MATRICES FROM INPUT FILE. DO NOT RUN E927. 

= 1 Rm E927 TO CALCUUTE ROTOR MATRICES -:NOT AV«IUSBLe *. :r ’ 

2 NCN ROTOR CONNECTION NODE NUMBER ^ .ft 

EULER PITCH ANGLE - ROTATE SECOND ABOUT THE Y-AXIS (DEGREES) 0.0 

EULEJ, mL_AN.G|.E_- R0TATE_JHIBD_ftE!fi5rtJ^ ISEJiSEig) _(bjp 

EULER YAM ANGLE - ROTATE FIRST ABOUT THE Z-AXI3 (DEGREES) 0.0 

(q) Component RE2 Data 


3 THETA 

± ^L_ 
5 XSI 


Figure 17. Continued 



CaaXMEHTiROTCREU xmuhh .. {USTIC ROTOR TTPI J «mh»wh ^ EUimTI 


I IROFIL. REIO FILE. FUR ... _.. ..... ... .... ... . . ... _.. .... ... ... 1 

> I READ ROTOR JIIPEOAHCE AW HUB FORCE VECTOR FROM SAVE FILE IT IMOS FROGRAIII 
■ 1 REAS ROTOR IHPEOAIICE AW HUB FORCE VECTOR FROM IMIT It CUSER THPUTI 
WTE I • INPUT DATA FORMAT IS tSlt.A. 

WTE t. r IF IROFXL>l> .XHHC.IUIC SI AW IMRSS | UK ESI ARC SET TO EERO UfrCRNAUTf - - - 

t ISTFOR SET ROTOR HUB FORCE FLAS 

• S LEAVE HUB FORCE VECTOR UlCHAHGEO 

..... ..... .1 .HULL THE HUB FORCE VECTOR IH THE CODE — ... 

1 IDPFOR DISPUT HUB FORCE FLAG 

• S DISPUT HUB DISPLACENEHTS AW SAVE OISPUCEHEHTS FOR PIOTTIW 

T...1...0ISPUT.HUB .INTERFACE . FORCES AW HOHEHTS AW SAVE FOR PLOT ROUTINE 

I ABOVE FUG 13 SET TO ZERO INTERHALLT FOR ACCELERATIONS OUTFUTI 


4 

NCM 

ROTOR HUB CGtllECTIOIJ NODE NUMBER 


1 

S 

Iweta'” 

EULER PITCH WCIE - ROTATE SECOIO ABOUT THE Y-MIS lOECREESI 

sTs ~ 


4 

PHZ 

EULER ROLL A»IGLE - ROTATE THIRD ABOUT THE X-AXZ9 IDE6REESI 

I.S 




EULER' TAU Aiisu - ROTATE FIRST ABOUT THE Z-AXIS IDEGREESI 

I.S 


0 

IHHC 

HHC FUG 


0 


« HHC.WT ACTIVE— 

> I HHC ACTIVE 

9 ITNHC READ THETA FOR HHC FUG S 

!...« . CONFUTE. OFTIIUL.THETA. ST . QUAORATIC.niNimZATIOH 

> I READ HHC INPUT AISLES - LOCATIONS II TO U BELCH 

IS ISTRSS DISPUT BUOE STRESSES FLAG - SET INTERHALLT. FROH GAOS INPUT LOCATION 991 S 

». »_ 00 WT. DISPLAT STRESSES t G4IIS LOCATION 991. IS ZERO 

■ 1 DISPUT STRESSES - G400 UKATIOH 911 IS GREATER THAN ZERO 

11 THItllC HHC INPUT AWU - IH-1I1H COSINE CONPONEHT IDE6REESI. H ■ MISER OF BLADES I.S 

IE iHStis~ HHC liiwir «siE~rN-i iiii sine coiiKiitOT^ ToinEEsi.' h ouser of budes Tb 

IS THC HHC INPUT ANGLE - I H ITH COSINE CONPONEHT IDEGREESI. H ■ MISER OF BLADES S.S 

“ T4~THtts“ ■»(C"IHPW'AWLE'*-T iTITH sine MWOHEHt «DEGRtE8l.“ H • MISER «F BLADES T.t 

IS THNPIC HHC IiaVT ANGLE - INlllTH COSINE CONPONEHT IDEGREESI. H ■ MISER OF BLADES S.O 

ii~TOiiPis “HHC iNfuf AWLi ^SVioBER OF Sub iTi 

17 UTHim HEIGHT FOR IH-IITH CONPalEHT OF HHC INPUT AIKLES. H • MISER OF BLADES 9.D 

~fs 'iilHil ■" HEIGHT "FOR-iTI ml CONPOliEHf OF^ h“» miser of' BLAOEB t.S 

19 UTHNPl HEIGHT FOR IHlllTH CONPONEHT OF HHC INIVT ANGUS. H « MISER OP BLADES S.l 


(r) Component RE3 Data Input 



CD 

Ca> 


Figure 17 


Continued 



HUB irtPEDANCE MATRIX FOR HHC PITCH ANGLES (UNITS = DEG) 


- ASS 


LONG. 

SHR 

(♦» 

AFT) 

COS 

-6.122 

UT. 

SHR 

(♦, 

STRBO) 

COS 

5.983 

VERT. 

SHR 

(♦, 

UP) 

COS 

-.9781 

ROLL 

MOMT 


STRBO UP) 

COS 

-3.113 

PITCH 

MOMT 

(♦f 

NOSE UP) 

COS 

1.450 

YAH 

MOMT 

(♦f 

OMEGA DIR.) 

COS 

-.5519 

LONG. 

SHR 

(♦. 

AFT) 

SIN 

-6.731 

LAf.” 

SHR 

‘(♦," 

STRBO ) 


-4.019 

VERT. 

SHR 

(♦. 

UP) 

SIN 

4.041 

ROLL 

MOMT 

(♦f 

STRBO UP) 

SIN 

-1.292 

PITCH MOMT 

i*. 

NOSE UP) 

SIN 

-3.460 

YAW " 

homt" 

(♦. 

OMEGA OIrI) 

SIN 

-.2991 


A4S 

. ASS- 

- B3S 

B4S 

■B5S 

-2.269 

-10.49 

6.087 

.3472 

6.623 

-6.202 

-7.944 

2.679 

1.216 

-9.454 

-12.70 

.1284 

-4.017 

3.210 

-3.038 

.6585 

-1.804 

1.172 

-2.267 

1.968 

-.2336 

-2.126 

3.585 

.2067 

-1.958 

6.095 

-3.664 

.1123 

1.048 

.7380 

.1032 

-6.567 

-6.156 

-1.600 

-9.874 

-i.996 

8.^ 

5.674 

-6.713 

-8.297 

-2.930 

3.490 

-.7051 

-12.99 

.1188 

2.294 

-1.932 

-3.076 

.5805 

-1.643 

-.3038 

* 1.856 

1.376 

-.1772 

-2.101 

-1.254 

-.9614 

-.6519 

6.058 

-3.819 


(s) Hub Impedance Matrix for HHC Pitch Angles for RES Component 


Figure 17. Continued. 



i I I P 


»»*«hhhhhi SUmARY OF HHC RESULTS iijhhhhhuhhi 

HZX HEIGHT. FOR ..X. RESPONSE-- - 

HZY HEIGHT FOR Y RESPONSE 

HEIGHT...FOR , Z ..RESPONSE 

HEIGHT FOR THETAl RESPONSE 

HEIGHT...FOR TOETAZ .RESPONSE — - 

HEIGHT FOR THETAS RESPONSE 

HTHNMl HEIGHT-FOR-lHrl JTH COHPONEMI OF HHC INPUT- ANGLES -H = NUtCER OF BLADES -- 

HTHN HEIGHT FOR I N )TH CONPONENT OF HHC INPUT ANGLES. N = NUMBER OF BLADES 

HTHNPl . HEIGHT- FOR- (MtUTH COMPONENT-OF HHC INPUT, ANGLES. N = NUMBER OF BLADES 

THSTR(l) OPTIMAL ANGLE - CN-DTH COSINE COMPONENT (DEGREES). N = NUMBER OF BLADES 

THSTRCZ) OPTIMAL ANGLE ..T. I N )TH COSINE. COMPONENT- (DEGREES). N s NUMBER OF BLADES 

THSTR(3) OPTIMAL ANGLE - (N*1)TH COSINE COMPONENT (DEGREES). N s NUtSER OF BLADES 

THSTR(4)_ OPTIMAL. ANGLE.. -.tN=I)TH SINE COMPONENT (DEGREES).-. N s. NUMBER OF BLADES - 

THSTR(5) OPTIMAL ANGLE - ( N )TH SINE COMPONENT (DEGREES). N s NUMBER OF BLADES 

THSTRIfil OPTIMAL.ANGLE_-..INtl)TH SINE COMPONENT. (DEGREES) N .S NUMBER OF BLADES — 

QUADRATIC PERFORMANCE INDEX 

(t) Summary of HHC Results for RE3 Component. 


Figure 17. Continued. 


- 2.00000D*04 
Z.OOOOODt04 
1.50000D+05 
0.0 

0.0 - - 

0.0 

— X.OOPOOD-02 - 
l.OOOOOD-02 
-.... l.OOOOOD-02 
3.50&81D-01 

6.126160-02 

2.52027D-02 

_-4.28766D-01 

2.16364D-01 

^1.416750-01 

1.165750-06 



00 

CT> 


»*•»»«»«»»*» BLADE STRESS HARHGNICS FOR AEROEUSTIC ROTOR RES ihuhuhhhih 

HARMONIC =1 2 3 4 5 6 7 8 

FLAT STN 1 COS -0.26050«01 0.11960*01 -0.89490*00 0.47740-01 -0.13580*00 0.90990*00 -0.60890*00 -0.32310*00 

FUT STN 1 SIN 0.36300*01 -0.12230*00 0.92150-01 0.10090*01 -0.11500*01 -0.75730*00 -0.25790*00 0.15810*00 

FLAT STN 2 COS 0.90620*00 0.32000*01 -0.53440*00 0.28690*00 -0.49800*00 -0.38660-01 -0.11970*00 0.24620*00 

FLAT STN 2 SIN 0.10440*02 -0.74170*00 0.13930*01 -9.39870*00 0.88980*00 -0.45130*00 -0.16030*00 0.27800*00 

FUT STN 3 COS -0.79160*01 0.56880*01 0.19610*01 -0.14110*01 -0.38430*00 0.51950-01 -0.21620*00 -0.805^-01 

FUT STN 3 SIN 0.94690*01 0.22030*00 -0.95110*00 -0.77050*00 0.39520*00 0.83100*00 -0.13510*00 -0.12550*00 

UGE STN 1 COS -0.91560*01 0.14750*01 0.94410*00 -0.35980*01 -0.43500*01 0.11690*00 0.87280*00 0.56970-01 

EDGE STN 1 SIN 0.14230*02 -0.19780*01 -0.85540*00 0.43410*01 -0.25310*01 -0.18640*01 0.12250*01 -0.28510*00 

EDGE STN 2 COS -0.54310*01 0.33680*01 0.10460*01 -0.29700*01 -0.42910*01 0.61230*00 0.10200*01 -0.36760-01 

EDGE STN 2 SIN 0.93680*01 -0.24930*01 -0.18220*01 0.55290*01 -0.16630*01 -0.27680*01 0.14260*01 -0.26310*00 

cbs '-oTiiii d.64610*00 -0.11040*01 -0.126SO*01 0.43830*00 0.37910*00 -0.10880*00 

EDGE STN 3 SIN 0.17330*01 -0.84030*00 -0.97490*00 0.18330*01 -0.63760*00 -0.10730*01 0.45160*00 -0.96770-01 

Tons STN i cos ” b!79B8b*00 -0.96710*00 0.14910*01 -0.61120*00 0.13040*01 -0.44450-01 -0.61640-01 -0.61120-01 

TORS STN 1 SIN 0.30150*01 -0.61990*00 0.82950*00 0.76090*00 0.69880-01 -0.24770*00 0.52820-01 0.66310-01 

fciis" i”’c03 0.11670^01 -0.92590*00 0.io990*01 -0.66980*00 0.10400*01 -0.73660-01 -0.76550-01 -0.25940-01 

TORS STN 2 SIN 0.29620*01 -0.44620*00 0.79050*00 0.47050*00 0.27950*00 -0.17760*00 0.37320-02 0.78290-01 

TmS~STN 3 7;OS 0.24480*di -0.66840*00 0.34230*00 -0.26790*00 0.57650*00 -0.59660-01 -0.45770-01 0.97550-02 

TORS STN 3 SIN 0.29560*00 -0.27710*00 0;65640*00 -0.74650-01 0.39800*00 -0.94040-01 -0.22670-01 0.29920-01 

ol290OT*00 -0.347^ l).6327iD*00 -0.24860+00 "oisOWD+OO -0.89ho-02 -0.29190-01 -0.22110-01 

PUSH STN 1 SIN 0.91460*00 -0.22470*00 0.15410*00 0.40070*00 -0.10670*00-0.98060-01 0.23460-01 0.33990-01 


(u) Blade Stress Harmonics for RES Component. 


Figure 17. Continued 



HALF PEAK TO PEAK VALUES 


FUT 

STN 

1 

0.6387Dt01 

FLAT 

STN 

2 

O.llXZDtOZ 

FUT 

STH 

3 

0.1787Dt02 

EDGE 

STH 

1 

0.Z477D«02 

EDGE 

STH 

2 

0.2103Dt0Z 

EDGE 

STH 

3 

0.63070*01 

TORS 

STH 

1 

0.55060*01 

TORS 

STH 

2 

0.67300*01 

ims~ 

STH 

~3 ~ 

0.37260*01 

PUSH 

STH 

1 

0.20180*01 


(v) Half Peak-To-Peak Values of Stresses and Pushrod Load 
for RE3 Component 
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00 

00 


... COMPONENT :TIMEH1_S1 ; < Hmx»nnnnnt NEWMARK .IMIEGRATJON...MnHPfl_.MS*S ELEhfNT.? ._ 3 


... HIRMARKJ^ETA Fjft.CT.pPt.,PAJ«E.., PERMITTED... =„Q..OJP.fl _ 2.S0000D-01 

2 DELT TIME INCREMENT (SECONDS) l.OOOOOD-02 

.3-_IM*>L__.TLM.E..LIMnJL0R..IMTEG54JIp|^ l.OOOOODtOD 

4 OMEGA FORCING FREQUENCY (HERTZ) 6.00000D+00 

. ... ..5 .IDE.BU5,_. P.EBUP. SEJ.EQTQR. 0 

= 0 ==> NO DEBUG PRINTOUT 

= 1 ==> TRACE MATRIX ASSEMBLY AND SOLUTION 

6 IRS TRT RESTAR T F LAG q 

= d ==> NO RESTART 
= 1 ==> RESTART 


(w) Component THl Data 


Figure 17. Concluded 





DCBU6 PHIKT SUBROUTINE 

HO ( If 1|3 l.OOOOOD^Ol HO ( 3> 31s l.OOOOOO^Ol 

HO ( 6* 6I« l.OOOOOOfOO HO . I 7» 7)s 1.000000*00 . 

CO ( 1. l)s 1.200000*02 CO I I, 2)s-1.2Q000D*02 

KO I 1. 1)3 2.500000*04 KB < 1* 2)>-2.500000*04 

H f !• 1)> 1.000000*01 

C I If 1)3 1.200000*02 C C If 2)3.1,200000*02 

K ~ f~ if II* 2isOodoD*04 K'”i” u 2)»-2.50000D*04 

HO ( If ll3 1.200000*01 




HO f 

4f 

4)3 1.000000*01 

HO 

C 

Sf 

5)3 

1.000000*00 

CO ( 

2 . 

1)3-1.200000*02 

CO 

( 

2* 

2)3 

1.200000*02 

KO ( 

2f 

1)3-2.500000*04 

KO 

f 

2f 

2)3 

2.500000*04 


C 1 

1 2, 

1)3-1.200000*02 

c 

( 2, 

2)3 

1.200000*02 

K 1 

1 2, 

1)3-2.500000*04 

K 

( 2. 

2)3 

2.500000*04 


H C If 1)3 1.000000*01 H I 2f 2)* 1.200000*01 

C I 1, 1)3 l.lDOOOOtOl C I 1| 2)3-1. 200000*02 C ( 2t 1)3-1.200000*02 C I 2i 2» 1.200000*02 
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Figure 20- Base Program Segmentation Structure. 
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Figure 21. Flow Chart for the Main Program. 
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Figure 22. Flow Chart for Subroutine INPCDC. 
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Figure 24. Flow Chart for Subroutine PINPUT 
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Figure 25. Flow Chart for Subroutine LINK. 
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Figure 26. Flow Chart for Subroutine RECOVR. 
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Figure 27. Flow Chart for Subroutine TRNSFM 
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Figure 28. Flow Chart for Subroutine COORD. 
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Figure 29. Flow Chart for Subroutine LABEL. 
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Figure 30. Flow Chart for Subroutine COMPUT. 
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Figure 31. Flow Chart for Subroutine UPDATE 
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Figure 37. Flow Chart for Subroutine MATDl 













Figure 38. Flow Chart for Subroutine RESTOR. 
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Figure 39. Flow Chart for Subroutine RESFRl . 
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Figure 42. Variation of Working Storage with CPU Time. 
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l._,c_ (.V(101,XSI._ 1, 0019 

1 C (V(111,GAMMAC1,(V(121,DGAMMC1 0020 

1 C (V(131,GAMMAS1,(V(141,DGAMMS1 0021 

1C* VERTICAL * (V( 11,XNCN 1,(V( 21, XN 1,(VC 31, MB 1, 0022 

1_C _ ..( V(._41,RR... _1,(V( S1,R 1,(V( _61,ZETBET), 0023 

1 C (V( 71, OMEGA 1,(V( 81, THETA 1,(V( 91, PHI 1, 0024 

1 C (V(101,XSI 1, 0025 

1 C (V(111,BETA0 1,(V(121,0BETA01 0026 

1_.C ( V( 13 1 , BETAC ...1., C V( 14 1 .DBETAC 1 0027 

1 C (V(151, BETAS 1,(V(161,0BETAS1 0023 

1C 0029 

7 V(41 = V(41 / 12.0 0030 

.7 V(51. .V(5). /_12.0_.. .... .0031 

7 11=14 0032 

7 IF(N.GT.ll 11=16 0033 

7 DO 10 1=8,11 0034 

7 V(I1=V(I1/RC 0035 

4 10 CONTINUE 0036 

1 C 0037 

7 RETURN 0033. 

7 El®. _ .. 0039 

(A1 ROUTINE CINBFl 


Figure 43 . Fortran Subroutines for Component BFl . 
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START 

...COL ,-".r+-r:~l + -2 +— ~3 + 4—--+ -5 *——6 7 +- 

7 SUBROUTINE RECBF1(V>NG.XDEP>DXDEP>KN) 0001 

1 C Mm t wMKWWMi t MMMiimn t KmoMMmunmMiDcioiincicMWMmticK i iinHiinomKiiii i iicmtmiidtmcKKmotmc ooOE 

" 1 C» ^RoijfriNi NAME : RECBFl »0004 

1 C» *0005 

1C* FUNCTIONS *0006 

. _1_,.C* RECOVER BIFILAR DISPUCEMEMTS AND VELOCITIES... — .... - - — *0007 

1 C* FOR RESTART WITH TIME HISTORY SOLUTION *0006 

1 C* FOR INPLANE BIFILARS IF KN=1, *0009 

1 C* FOR VERTICAL BIFILARS IF KN=2 *0010 

1 C **************** * ***«***«M*** ** ***«*M * ***I H I******» * *****M*M » * ’ * *** *** * ** 0012 

1 C 0013 

7 DIMENSION V(1),XDEP(11,DXDEP(1J 0014 

_ ...1 ..C - .0015 

1 C INITIALIZE THE COORDINATE INDEX FOR THIS ELEMENT 0016 

1 C 0017 

7 NN = NG ♦ 1 0018 

I.....C, 0019 

1 C INITIALIZE THE BIFILAR DISPUCEMENTS AND VELOCITIES 0020 

1 C 0021 

7 V(ll) = XDEP(NN) 0022 

7_ V( 12 )....= .0XDEP(NNJ._ ...„ : 0023 

7 NN = NN + 1 0024 

7 V(13) = XDEP(NN) 0025 

7 V(14) = DXDEP(NN) 0026 

7 IF(KN.LT..2.)..GO...TO_10 0027 

7 NN = NN + 1 0028 

7 V(15) = XDEP(NN) 0029 

7 V(16) = OXDEP(NN) 0030 

_ ...1 .c ..... _ .; ..._ 0031 

1 C INCREASE THE COORDINATE IM3EX TO ACCOUNT FOR THE HUB COORDINATES 0032 

1C 0033 

4 10 CONTINUE 0034 

I NS,_=,.NN.+..6 0035 

1C 0036 

7 RETURN 0037 

7 END 0038 


(B) ROUTINE RECBFl 


Figure 43. Continued. 
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SUBROUTINE PINBFKV.IELEM.KPRINT.KN) 


SUBROUTINE NAME : PINBFl 


unction: 

..PRIMT...THE INPUT 

FOR INPUNE BIFILARS IF KNsI> 
FOR VERTICAL BIFILARS IF KN=2 


DIMENSION V(l) 


..IF,(KN.EQ.1).WRITE(KPHINT,1000). lELEM 
IF (KN.EQ.2) MRITECKPRINT,2000) lELEM 
NCNl = V(l) + .001 
MRITEtKPRINT.lOOX) NCNl 

„N = . V( 2 ), .001 . .. 

WRITE(KPRINT,1002) N 
WHITE CKPRINT, 1003) V(3) 
URITE(KPRINT>1004) V(4) 

..WRITE(KPHINT,1005). VtS),. . 

IF (KN.EQ.l) WHITE(KPRINT,1006) V(6) 

IF (KN.EQ.2) KRITE(KPRINT,2006) V(6) 
HRITEfKPRINT.lOO?) V(7) 

_WRITE( KPRINTilOOB ) . ..V( 8 ) 

WHITE(KPRINT,1009) V(9) 
WHITE(KPRINT,1010) V(IO) 


IF(KN.GT.l) SO. TO. 
WRITE(KPRINT,1011) 
WRITE(KPRINT,1012) 
WHITE(KPRINT,1013) 
.lffiITE(KPRINT,l014). 
GO TO 20 
CONTINUE 

WRITE(KPRINT,2011) 
WRITE (KPRINT, 2012) 
'WRITE(kPHlNt,2dl3) 
WRITE(KPRINT,2014) 
WHITE(KPRINT,2015) 
WRITE( KPRINT, 2016) 
"CONTIhWE 


10..... _ 

V(ll) 

V(12) 

V(13) 

_.y.J19).. 


FORMAT STATEMENTS 


) F0RMAT(55H1C0MP0NENT:BIFILAR « 

1 , 55HAL LINEAR BIFIUR *** 

2 , 110 ,//) 

J F0RMAT{55H1C0MP0NEMT:BIFILAR « 

1 ■■■TSSHAL linear BIFIUR 

2 , 110 ,//) 

I FORMATC55HO 1 NCNl CONNECTION NODE NUMBER 


END 1 (hO) 


Figure 43. Continued. 
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START 


6 1 

2 1002 FORHAT(55HO 

.. 6 „ 

2 1003 FORMAT(55HO 

6 1 

2 1004 FORMAT(55HO 

. 6 ., 1. .... ... ,35H 

2 1005 FORMAT(55HO 

6 1 

2 1006 FORMAT(55HO 

,..6_ _.l 

2 2006 FORHAT(55HO 

6 1 

2 1007 FOHMAT(55HO 

...6 1 

2 1008 FORMAT(55HO 

6 1 

2 1009 FORMAT(55HO 

.6 1 

2 1010 FORKAT(55HO 

6 1 

2 1011 FORMAT(55HO 

,. 6 _ ..... .1 _. .. 

2 1012 FORKAT(55HO 

6 1 

2 1013 F0RrUT(55H0 

„6 1 - 

2 


1014 FORMAT(55HO 
6 1 

2 2011 F0RKAT(55H0 

6 _ . 1 .. 

2 2012 FORt1AT(55HO 

6 1 

2 2013 FORMAT(55HO 

.6 _1.... 

2 2014 FORMAT! 55H0 

6 1 

2 2015 FORMAT! 55H0 

_ 1 _ 

2 2016 FORMAT! 55H0 

6 1 

1 C 

7 RETURM 


,40XiI15) 

2 M TOTAL NUMBER OF BIFILARS !ND) 

,40X.I1S) 

3 MB BIFILAR MASS ! SLUGS) 

,40X.1PE15.5) 

4 RR DISTANCE FROM CENTER OF BIFILAR TRACKING 

.HOLE TO. CENTER OF ROTATION ! INCH) ,5X,1PE15.5) - - 

5 R EQUIVALENT PENDULUM ARM !INCH) 

,40X,1PE15.5) 

6 ZETGAM INPLANE BIFILAR DAMPING RATIO C/CCRIT !N 

3HTD ) . T96 , 1PE15 . 5 ) - 

6 ZETBET VERTICAL BIFILAR DAMPING RATIO C/CCRIT !N 

,3HTD),T96,1PE15.5) 

7 OMEGA BIFILAR ROTATION SPEED ! RAO/SEC) 

.._ .... ... ...,40X.1PE15.'5) 

8 THETA EULER PITCH ANGLE AT NODE - ROTATE SECO 
,30HNO ABOUT THE Y-AXIS !DEGREE5) ,10X,1PE15.5 ) . 

9 PHI EULER ROLL ANGLE AT NODE - ROTATE THIRD 

,30H , AB0UT..THE X-AXIS I DEGREES ) ,10X .IPEIS.S ) 

10 XSI EULER YAW ANGLE AT NODE - ROTATE FIRST 

,30HABOUT THE Z-AXIS !DEGREES) ,10X,1PE15.5) 

11 GAMMAC COSINE COEFFICIENT OF INITIAL PENDULUM C 

,30HYCUC..DISPIACEMENT !DEGREES) , .,10X,1PE15.5) 

12 DGAMMC COSINE COEFFICIENT OF INITIAL PENDULUM C 

,30HYCLIC VELOCITY !DEGREES/SEC) ,10X,1PE15.5) 

13 GAMMAS SINE COEFFICIENT OF INITIAL PENDULUM CYC 

iSOHCLIC.DISPUCEMENT...! DEGREES). jlOX,lPE15.5) 

14- OGAMMS SINE COEFFICIENT OF INITIAL PENDULUM CYC 
,30HLIC VELOCITY !DEGREES/SEC) ,10X,1PE15.5) 

11 BETAO BIFILAR MASS COLLECTIVE INITIAL FLAPPING 

,.30H ANGLE...! DEGREES). ,10X,1PE15.5).. 

12 DBETAO BIFILAR MASS COLLECTIVE INITIAL PUPPING 

,30H VELOCITY (DEG/SEC) ,10X,lPE15-5) 

13 BETAC BIFIUR MASS COSINE COEFFICIENT OF INITI 

,,30HAL..fUPPING ANGLE. !DEGREES)_ ,10X,.-1PE15.5) 

14 DBETAC BIFILAR MASS COSINE COEFFICIENT OF INITI 
,3CHAL PUPPING VELOCITY ! DEG/SEC ),10X,1PE15. 5) 

15 BETAS BIFILAR MASS SINE COEFFICIENT OF INITI 

,.30HAL_FUPPING.ANGLEJDEGREES) „,10X,1PE15.5) .. 

ii5 DBETAS BIFILAR MASS SINE COEFFICIENT OF INITI 
,30HAL FUPPING VELOCITY !DEG/SEC),10X,1PE15.5) 


END 


0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 
0074' 

0075 

0076 

0077 

0078 

0079 

0080 
0081 
0082 

0083 

0084 
.0085 
0086 

0087 

0088 
.0089 

0090 

0091 

0092 
...0093 

0094 

0095 

0096 
...0097 

0098 

0099 

0100 


!C) ROUTINE PINBFl 


Figure 43. Continued. 


128 



START 

_„COL rr+--~l- + 2— -+—,—3 -* — — 4— -+ 5 + 6 -+— -7 +- 

7 SUBROUTINE CRDBF1(DATA>HLOCAL>NNODEL>NCNL.EULER>KN) 0001 

1_..C»* K0003 

1 c» SUBROUTINE NAME : CRDBFl »0004 

1 C» »0005 

1 C» F U N C T I 0 N : «0006 

INFORM THE MAIN PROGRAM OF THE REQUIRED NUMBER OF COORDINATES - *0007 

1 C» FOR INPLANE 8IFILARS IF KN=1, )»0008 

1 C» FOR VERTICAL BIFILARS IF KN=2 »0009 

1 c» *0010 

1_.„ C *>(illlKintm»m<<llMK)HH(MKl(IH(ICmCllKKm(MX)IMMMKMK)(XMl(ICM I H H HH<)HCKK)ll(MX)<im i ( ' )H()CIC)mim 0011 

1 C 0012 

7 REAL MB 0013 

7 DIMENSION DATA( 1 } ,V( 10 ) 0014 

...7. ...DIMENSION NCNU1),EULEH(3»1J 0015 

7 EQUIVALENCE(V( D.XNCNl ),(VC 2),XN »,(V( 3), MB ), 0016 

6 1 (V( 4),RR ),(VC 5),R ),(V( 6),ZETGAM), 0017 

6 2 (V( 7), OMEGA ),(V( 8), THETA ),(Vl 9), PHI ), 0018 

6 3 ( VCIO l.,XSt_. } :. 0019 

1C 0020 

1 C RETRIEVE INPUT DATA 0021 

1 C 0022 

D0„20,..l=ltl0 .... ..0023 

7 V(I) = DATA(I) 0024 

A 20 CONTINUE 0025 

1C 0026 

l__C_J.PECIFr...NUM8ER_OF...COORDINATES i_CONNECTION. NODE. NUMBERS, t. EULER ..ANGLES0027 

1 C 0028 

7 NLOCAL = KN 0029 

7 NNODEL =1 0030 

_. _ __ 

7 NCNL(l) = 0 0033 

7 NCNL(2) = 0 0034 

7 HCNH.3)_=..0 0035 

7 KNM5 = KN - 5 0036 

7 DO 10 I = KNM5.KN 0037 

7 NCNL(I) = NCNl 0038 

,_J*_ ._10...CP!Tri.!|JE_ ..... _ .....0039 

1C - . - 

7 EULERCl.l) = THETA 0041 

7 EULER(2>1) = PHI 0042 

7 EULER(3t;y_= >S 0043 

i c '0044 

7 RETURN 0045 

7 END 0046 


(D) ROUTINE CRDBFl 


Figure 43. Continued. 
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START 

7 SUBROUTINE RENBFKV.NS.XDEP.DXDEP.KN) 0001 

1 C XM i d i mm i no t Mi t mowimiotumiiHDnuwMMm i im i m x iMdniinDmmmimiim H CMMKimMiniini i iKMii 31 0002 

1 C» »0003 

I C» SUBROUTINE NAME : REN8FI *0004 

1 C3f *0005 

1C» function: »0006 

._.I C» . INITIALIZE. DISPUCEMEHT3 AM) VELOCITIES — *0007 

1C* FOR RESTART WITH TIME HISTORY SOLUTIONS *0008 

1 C* FOR INPLANE BIFILARS IF KN=1, *0009 

1 C» FOR VERTICAL BIFILARS IF KN=2 *0010 

..._...* 0011 . 

1C OOlS 

7 DIMENSION V( 1) ,XDEP(1),DXDEP( 1) 0014 

_ ...1 ..c: .0015 

1 C INITIALIZE THE COORDINATE INDEX FOR THIS ELEMENT OOlfr 

1 C 0017 

7 NN = NG ♦ 1 0018 

i “c iNifiALizE "™’e ifFiEiR'iSsp"^^ ~ 0020 

1C 0021 

7 XDEP(NN) = V(ll) 0022 

7_ ..0XDEP.(NN)=.y(12,) .... ..._ ..... _ .... __ ..._ 0023 

7 NN = NN ♦ 1 0024 

7 XDEP(NN) = V(13) 0025 

7 DXDEP(NN)= V(14) 0026 

_,.7 JF(KN.LT.,.2.)._OO...TO..JO ^ .0027 

7 NN = NN + 1 0028 

7 XDEP(NN) = V(15) 0029 

7 DXDEP(NN)= V(16) 0030 

...1 ,C ...... ..... _... _ ..... ._ .... _ ..... 0031 

1 C INCREASE THE COORDINATE INDEX TO ACCOUNT FOR THE HUB COORDINATES 0032 

1C 0033 

4 10 CONTINUE 0034 

;.....7 NO JE...NN J;_..6. 0035 

1 C 0036 

7 RETURN 0037 

7 END 0038 


(E) ROUTINE REN8F1 


Figure 43. Continued. 
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START 

-COL + 2— — +— r-3 + 4 T-+ 5 + 6 7— -4- 

7 SUBROUTINE NAMBFKALPHA.NLOCAU 0001 

1 Cx.. *0003 

1 CX SUBROUTINE NAME S NAMBFl X000<4 

1 Cx X0005 

ICx function: xoo 06 

.._1„...CX_ NAME. VARIABLES FOR INPLANE. BIFILARS - — X0007 

1 CX . XOOO8 

1 CXXX XX X XXXX X XXXX X X XXX- XXXXXXXX X XXXXXXXXXXXXXXXXXXKKXXKXXXKXXKXXXXXXXXXKXX 0009 

1 c 0010 

...7_ DIMENSION. ALPHACS.l) ,GAMMAC( 2 ) .GAMMASt 2 ) , 0011 

6 1 Xt2),Y(2),Z(2),THTX(2)tTHTY(2J,TMT2(2) 0012 

7 DATA GAMMAC/AHBFl ,4HGAMC/t 0013 

6 1 GAMMAS/4HBF1 ,4HGAMS/, 0014 

...6 _ 2... ... X. ..... /4HBF1 , ,4HX 0015 

6 3 Y /4HBF1 ,4HY /, 0016 

6 42 /4HBF1 ,4H2 /, 0017 

6 5 THTX /4HBF1 ,4HTHTX/, 0016 

_ 6 6 THTY...y4HBFl....i4HTNTY/, 0019 

6 7 THT2 /4HBF1 .4HTHT2/ 0020 

1 C 0021 

1 C RETURN THE COORDINATE NAMES 0022 

~"7~ NLOCAL ^8 ” ~~ ' 0024 

1C 0025 

7 00 10 1=1.2 0026 

I1...E...I...+ ..1 0027 

7 ALPHACIl, 1) = GAMMAC(I) 0028 

7 ALPHACIl. 2) = GAMMAS(I) 0029 

7 ALPHAdl, 3) = X(I) 0030 

_ 7 ALPHA(I1,_4)....=_Y(I) _. __ 0031 

7 ALPHAdl, 5) = 2d) 0032 

7 ALPHAdl, 6) = THTXd) 0033 

7 ALPHAdl, 7) = THTYd) 0034 

7 ALPHA.dlx..S)_r....THT2(.I) 0035 

4 10 CONTINUE 0036 

1 C 0037 

7 RETURN 0038 

_7 END ....0039 

(F) ROUTINE NAMBFl 


Figure 43. Continued 



START 

7 SUBROUTINE CPBFI(DATA,NC,M,C,K,F.IDIH) 0001 

1 CK wimmiii i iMKMXM i cKMMWi i MiiidciiiODtiiiDOcmcMmimnnoiKm H umi t mmMmomi i n H uniK ii mimmMit ooog 

1 Cx SUBROUTINE NAME : CPBFl XOOOA 

1 Cx X0005 

1 CX function; X0006 

1— CX .CALCULATE_THE MASS». STIFFNESS,. AND DAMPING MATRICES X0007 

1 CX AND THE FORCE VECTOR FOR INPLANE BIFILARS XOOOB 

1 CX K0009' 

I C WXKX X XMXKXXXXXXMXXXXMXXXXXXXX X KXMXX X X XX X X XMXXXXXXXMXKMXXXMXMXXXXXKX X XX - X OOlO' 

l_.C OOll! 

7 REAL K,M,M8 0012 

7 DIMENSION VdOl.DATAd) 0013 

7 DIMENSION M(IDIM,l),CtIDIM,l),K(IDIM,l),FdJ OOIA 

7 _ ...EQUIVALENCEtV(..l.),XNCNl_),(V(..2),XN J,(V( ..3),MB ), _ 0015 

6 1 CV( A),RR ),(V( 5),R ),(V( 6),2ETGAM), 0016 

6 2 (V( 7), OMEGA ),(V( 8),THETA ),(V( 9), PHI ), 0017 

6 3 (VdO).XSI ) 0018 

l—C 0019 

1 C RETRIEVE INPUT DATA 0020 

1 C 0021 

7 DO 20 1=1,10 0022 

7_ V(.I).=. OATAtI)_ .... ..._ _ 0023 

20 CONTINUE OOE^f 

1 C 0025 

1 C INITIALIZE MATRICES TO ZERO 0026 

l__C ^ ^ 0027 

7 DO 10 1=1,8 0028 

7 F(I) = 0.0 0029 

7 00 30 J=l,8 0030 

.7 Mdi,J)..= ,.0.0 .0031 

7 Cd,J) = 0.0 0032 

7 K(I,J) = 0.0 0033 

4 30 CONTINUE 0034 

4 ISL.CONHNUE ; 0035 

1C 0036 

1 C CALCUUTE THE MATRICES 0037 

1 C 0038 

7 NC_= 8 _ _ 0039 

7 "n" " = »< ■“ 0040 

7 OMEGAM = OMEGA x SC3RT(RR/R) 0041 

1C 0042 

7 Cd,l) = 2.0XZETGAMXOMEGAM _ 0043 

7 Cd,2) = 2.0xbMEGA “ 0044 

7 C(2,l) =-2.0XOMEGA 0045 

7 Ct2,2) = 2.0XZETGAMXOMEGAM 0046 

1_C _ _ „ 0047 

7 kt 1, i) = 0MEGAXX2X(RR/T?-I.di 0048 

7 Kt 1, 2) = 2.0XZETGAMXOMEGAMXOMEGA 0049 

7 Kt 2, 1) = -2.0XZETGAMXOMEGAMXOMEGA 0050 

7 Kt 2, 2) = OMEGAX«2X(RR/R-1.0J 0051 

1 t 0052 

7 Mtl.l) = 1.0 0053 

7 Mtl,4) = N/t2.0XR) 0054 

7 Mt2,2) = 1.0 0055 

7 Mt2,3) =-N/t2.0XR) 0056 

7 Mt 3,2 ) --MBXR 0057 

7 Mt3,3) = N«MB 0058 

7 Mt4,l) = MBXR 0059 

7 Mt4,4) = N«M8 0060 

' ~ ~7 ~ RETURN - - 0061. 

7 END _ 0063 

; t&) ROUTINE CPBFl 


Figure 43. Concluded 
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Figure 44. Vibration Isolator Models. 
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BASE PROGRAH CASE 5 
FIXED ABSORBER TINE HISTORV RESPONSE 
% NODE SHAPE - FX FORCE (COS B SIN) 


COORD UARIABLE tCEY 
a DELT DISP 1 
3 X DISP 2. 
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Figure 45. Base Program Results for Test Case 5. 
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BASE PROGRAn CASE 7 
1-P SINGLE ISOLATOR - OUTPUT IN G'S 
I TRANS, i a pus. NODES - F2 FORCE * 


COORD UARIABLE key 



FREQUENCY <HZ) 


Figure 46 . Base Program Results for Test Case 7 . 
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BASE PROQRAH CASE U 
1 a-D ISOLATORS - 3 FORCES | 3 nORENTS 
S TRANSt1ISS;ON Ik 6 FUSELAGE HOPES 

ISoMTigft ± 


COORD URRIABLE KEY 



as XBF ARPL. i 

a? VBF ARPL, i 

as ZBF ARPL. 3 

CASE 1 



FREQUERCY <HZ) 


(a) Isolator 1 Response 


Figure 47 . Base Program Results for Test Case 11 . 
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BASE PROGRAM CASE 11 
<4 3-D ISOLATORS - 3 FORCES & 3 MOMENTS 
6 TRANSMISSION % S FUSELAGE MOPES 



COORD UARIABLE 
39 ZBF AMPL. 

65 ZBF AMPL, 

8F AMPL. 
BF AMPL. 

CASE 1 


TSOLATDR 

1 . 

2 . 

3 
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(b) Vertical Response at Four Isolator Locations 


Figure 47. Continued. 
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e.e 


PROQRAh CASE li 

4 a-D ISOLATORS - a FORCES ( a nOHEHTS 
6 TRANSMISSION L R FUSELAOE NODES 



FREQUENCY <HZ) 


COORD VARIABLE H^SOlAToR 
as ZBF AMPL. d. 

65 ZBF AMPL. 

101 ZBF AMPL. 

137 ZBF AMPL. 


CASE 1 
gERO PeNsiTY 


z 

3 
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(c) Vertical Response at Four Isolator Positions 


Figure 47. Continued 
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BASE PROGRAM CASE U 
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ZBF 
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CASE 1 
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(d) Vertical Response at Four Isolator Positions 


Figure 47. Continued. 
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BASE PROGRAM CASE 11 
A a-D ISOLATORS - 3 FORCES i 3 nORENTS 
a TRANStllSSIOH L 6 FUSELAGE RODES 



COORP UARIABLE XSDlaTok 
137 ZBF AHPL. ± 

101 2BF ARPL. Z 

65 ZBF ARPL, 3 

as ZBF ARPL. A- 

CASE 1 
DEPTH REUERSED 


FREQUEHCV (HZ) 


6.00 tl0 


(e) Vertical Response at Four Isolator Positions 


Figure 47. Continued. 
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BASE PAOCRAn CASE 11 
4 e-D ISOLATORS > 9 FORCES k 3 nOREHTS 
6 TRANSRISSION k S FUSELAGE RODES 



COORD UARIABLE ISOL/ITok 


89 

?BF 

ARPL. 

1. 

S5 

ZBF 

ARPL. 

1 

101 

ZBF 

ARPL. 

3 

137 

ZBF 

ARPL. 

4- 


CASE 1 


FREQUENCV (HZ) 


e.d i.eo 8.00 3.00 <1.00 s.oo s.os xio 


(g) Vertical Response at Four Isolator Positions 


Figure 47. Continued, 
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BASE PROGRAM CASE U - TO 6<1 LBS 

^ e-D ISOLATORS - 3 FORCES B 3 MOMEHTS 
6 TRANSMISSION B 6 FUSELAGE MODES 

COORD UARIABLE ISol^ToR 
39 2BF AMPL. ± 
CASES 



FREQUENCY (HZ) 


(h) Vertical Response for Different Isolator Masses 


Figure 47. Concluded. 
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COORD 

UARIABLE 

KEY 

9 

RX 

ANPL. 

1 

11 

RV 

ANPL. 

Z 

13 

RZ 

ANPL. 

3 

CASE 

1 



±- j 1 


100 

80 



10 15 20 25 30 


3RD NODE FREO (Hit) 



(a) Force Results 


Figure 48. Base Program Results for Test Case 13. 



lSri.4CEM£-N.T J & S 


BASE PROCRAH CASE 13 
inPEDANCE - resuits in 
4 FUSELAGE NODES AND 6 ROTOR D.O.F, 


COORD uariable Ker 

9 X ANPL. 1 

It V ANPL. 2. 



3RD NODE FREQ iUi) 


(b) Displacement Results 


Figure 48. Concluded, 
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//ET<»73:;UI< JCO t<*0^5,SrAL.19»15),*CASC'Pr»ra - X73DcA', 

// CLASSED. Ti::ns3i.tt3CLcVi;L=a,o),RCCic:;=33oo:<» 

// MGGCUG3=T,K0TIFY=ET473 

»*» . R0T07CRAFT OYMAHIC nVSTSrrS Ar.'ALtSIG rZ2C'',\t\ 

n»» EXCCUTCS SIMVIO FRCClAfl WITH OVERLAY STPUCTU-E TCR CDC CC.*:rUTCR 

//STEPOl EXEC PCn=FIi:AL01 

.//STEPLIO .... DO 0$.MsET47J.SimO.LOAO,DISPsCHPi 

//FTOSFOOl 00 OUJCIY 

...READ INPUT DATA FPCtl WIIT 1 fCR NORMAL R I'M - 

//FTOlFCOl DO 0SHsETA73.SIMVID.DATA(5VIC3l),DlSPsCHR 

PEAO.C-400 INTUr.DATA FROM 11.CPRC3RAM) CR WIIT 12 lUGEO) 

• 

//FTllFOOl DO Dl'XMY <ss ALL CASES EXCEPT 1'4»15.X6 

• in'FTUFOOl 00 OSN5ETA73.SIir/IB.OATA(C'iOOIft?;$ ) ,OISP=i:i{R <== CA52 W 
*4»FT11F001 00 05H=ET<*73. SIMVIO. DATA(CAOOIi:XS)»DICP=Su7 <»* CASES 1SA16 - 

#«k- 

USE UJIIT 12 FOR ALL CASES EXCEPT FOR CASES 1A,15»1$ 

«w V 

//FTlCF0Ol_.DO.„.OSN5ET<»73.SILr/IB.DATAIC'VID13X).DISP=SHP. - 

k Xk-.- — — - 

««*» READ INPUT OATA FROM Will 10 FCR COUPLING WITH E-927 

kxk-- — — — 

//FTIOFOCI , DO __ DSMtETA73.5IMVIB.DATA(SVIB12X)»0ISP=SMR .... . 

«»K — 

«»*» fiCPMAl. OUTPUT IS WRITTEM OM FILE 6 



//FT06FC01 ,,D0 _... 5YS0UT=T 

» * k ........... ...... 

«** FL7XHE0 CLfTPl/r COES TO FILE 7 

kkk-..---- 

.//FT07F001, 00.. _DU{fMY_ ... _ ..... ... ... .... .. -. 

kXk 

kkk OUTPUT USED FOR PLOTTING IS WRITTEM TO FILE 2 

— — 

//FTC2F001 _ DO SYSOUTaT 

kkk— — — 

O'JTFUT USED FCR TINE HISTORY RESTART IS WRITTEN TO FILE 3 

kkk ~ 

//FTOSFCOl 00 SYSOUTsT <=== FOR ALL CASES EXCEPT CASE 6 

»*»i»FT0er001 00 OSWaET^73.SVlC04 X. data, OISP= n;CO, keep) <== CASE 6 

kk k------------ 

vXkkkUttitkkkNkkvH CPTICM3 FCR OUTPL/T RESULTS «k»«kMk«kki.'),kkkVi«k»K«kt«4*« 
#k««x»«« UNIT Z TO PLOT RESULTS IS USED AS AN EXAMPLE . k»kk»*kkkk>*kkkk« 

UH* OPTICM 1. WRITE OUTPUT TO PAPER 
kkkpTOEFOOl DO SYSOTTaT 

kkW>. 

'■'kkk * OPTICM 2. iOlITE OUTPUT TO A f.'CW FILE' ‘ 

»»»FT02F001 00 DSH=ETA73,PLOTSV01.{)ATA,CISP=(NEM,CATLG), 

«tin» L'MITaTSO, SPACEa C TPK . < 5 , 5 I) , 

OCD=(RECFliafB.U?eCLaa0.eLKSIZEs3120t 

»»» OPTION 3. i:stTE cumrr over am old file 

»UkFTO2FO0I CO 05 MsETh 73.PLOT3V01. OATA. OISPstOLO, KEEP) 

Kkk — - 

“«kk OPTION WRITE OUTPUT IM DACX OF A)l OLD FILE 
*«»FT02F001 00 DGMsET<i73.PLOTSVOl.OATA»DlSP=tMOO,KEEP) 

kkk— 

_k«k OPTICM 5, 00 NOT WRITE AMY OUTrUT . 

»*»»FT02F001 00 OUr-IY 

l»»ki»kNtfXkkk»«»«<(»k«kMkk»k»,«tN)*»«»«K«w«kk»»k«k»k»k»»«»kw»ik«k4(kkWk»kk»kkkWk» 

// <sss LAST JCL DIPUT CARO 


Figure 49. Base Program IBM JCL 



* /JOB 

^ /HOSEQ 

SVLINKB,T1400,CM300000. 

USER , . 

CHARGE, ,LRC. 

RFL<300000) 

GET.SVBIB. 

GET,SVB2. 

FTN( I=SVB1B,R=0,L=TAPE1 , E=NSVB1G0, A) 

FTN< I=SVB2, R=0,L=TAPE2,B=NSV2GO, A) 

REPLACE(NSVBIGO) 

REPLACE (NSV2GO) 

ATTACH( ALTMLIB/UN=LIBRARY) 

SEGLOAD(B=SVBABS) 

LDSET(LIB = ALTMLIB , PRESET = 2ER0,MAP = SBEX/LMAP ) 

LOAD(NSVBIGO) 

LOAD(NSV2GO) 

HOGD. 

REPLACE(SVBABS) 

DAYFILE, JCLOUT. 

REPLACE , JCLOUT. 

REPLACE(TAPE1=F0RTB1 ) 

REPLACE(TAPE2=FORTB2) 

REPLACE, LMAP . 

EXIT. 

DAYFILE, JCLOUT. 

REPLACE, JCLOUT. 

REPLACE(TAPE1=F0RTB1 ) 

REPLACE( TAPE2=FORTB2 ) 

REPLACE, LMAP . 

* /EOR 

SIMVIB GLOBAL ACOH VR , CBET A , CDAT A , CDE B UG , CF R 1 , CGEN , CGF 1 , Cl HOU T 
SINVIB GLOBAL Cl TER , CL AB EL , CMAX A , CN AMES , CP LOT , CR E2 , CR E3 
SIMVIB GLOBAL CR E3H 1 , CR E3H2 , CR STR T 

SUB2 TREE COMP UT- ( UP DATE , AS SMBL , S OL VE- ( S OLEGM , SOLFR T ) , R ESTOR , S AVE ) 
SUBl TREE INPCDC-( INPUT ,PIMPUT-(PRINT1 ,PRINT2n 
TREE S IMVIB-C SUBl , LINK, RE COVR ,TRNSFM, LABEL, SUB2) 

END SIMVIB 

* /EOF 

^ For card images operation; 

1. delete /JOB and /NOSEQ instructions (first 2 cards). 

2. replace /EOR with a 6/7/S multiple punch card. 

3. replace /EOF with a 6/7/S/9 multiple punch card. 


FIGURE 50. CDC JCL FOR INSTALLATION OF THE BASE PROGRAM BATCH MODE. 
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¥: /JOB 

* /NOSEQ 

SVL1NKI,T1^00,CM300000. 

USER,- . 

CHARGE, ,LRC. 

RFLC300000) 

GET,SVE1T. 

GET,SVB2. 

FTH<I=SVB1T,R=0,OPT=2,L=TAPE1 ,B=NSVT1G0, A) 

FTfK I=SVB2, R = 0,OPT = 2,L = TAPE2, B = NSVT2GO, A) 

REPLACE(NSVTIGO) 

REPLACE(NSVT2GO> 

ATTACH( ALTMLIB/UH=UIBRARY) 

SEGLOAD(B=SVOABS) 

UDSET(LIB=ALTMLIB,PRESET=ZERO,MAP=SBEX/LMAP ) 

LOAD(NSVTIGO) 

LOADcnSVT2GO) 

NOGO. 

REPLACE(SVOABS) 

DAYFILE, JCLOUT. 

REPLACE, JCLOUT. 

REPLACE(TAPE1=F0RTT1 ) 

REPUACE(TAPE2=F0RTT2) 

REPLACE, LMAP. 

EXIT. 

DAYFILE, JCLOUT. 

REPLACE, JCLOUT. 

REPLACE(TAPE1=F0RTT1 ) 

REPLACE(TAPE2=FORTT2) 

REPLACE, LMAP. 

* /EOR 

SIMVIB GLOBAL ACOH VR , CBET A, CDAT A , CDEBUG , CFR 1 , CGEN , CGF 1 , CIHOUT 
SIMVIB GLOBAL Cl TER , CL ABEL , CMAXA , CNAMES , CP LOT , CR E2 , CR E3 
SIMVIB GLOBAL CR E3H 1 , CR E3H2 , CR STR T 

SUB2 TREE COMP UT- ( UP DATE , ASSMBL , SOL VE- ( SOLEGH , SOLFR T ), RESTOR , SAVE > 
SUBl TREE INPCDC-(INPUT,PINPUT-(PRINT1,PRINT2) ) 

TREE SIMVIB- ( SUBl , LINK, RECOVR , T R NSFM , LABEL , SUB2 ) 

END SIMVIB 

* /EOF 

* For card images operation: 

1. delete /JOB and /NOSEQ instructions (-first 2 cards). 

2. replace /EOR with a 6/7/S multiple punch card. 

3. replace /EOF with a 6/7/S/9 multiple punch card. 


FIGURE 51. CDC JCL FOR INSTALLATION OF THE BASE PROGRAM INTERACTIVE MODE 



/JOB 

» /NOSEQ 

SIMVIB,T1^00,CMSOOOOO. 

USER , . 

CHARGE, ,LRC. 

RFL(300000) 

GET.SVBABS. 
GET,TAPE1=SVE12. 
GET,TAPE10 = E927h1AT. 


NOTE - REPLACE CARD ABOVE WITH FOLLOWING CARD(S>: 


GET,TAPES=SVE06X. 

< = = = 

FOR 

TIME 

HISTORY RESTART 

- SEE 

TABLE 

13, 

CASE 

6 

GET,TAPE11=SVB14X. 

11 

It 

II 

V 

FOR 

G400 

IMPEDANCES 

- SEE 

TABLE 

13, 

CASE 

14 

GET,TAPE11=SVB15X. 

< = = z 

FOR 

G400 

IMPEDANCES 

- SEE 

TABLE 

13, 

CASES 

isai6 

GET,TAPE12=SVB13X. 

< = = = 

FOR 

USER 

INPUT IMPEDANCES 

- SEE 

TABLE 

13, 

CASE 

13 


SVBAES. 

DAYFILE , JCLOUT . 
REPLACE , JCLOUT . 
REPLACE(OUTPUT=OUTR 
REPLACE(TAPE2=SVBPI 
EXIT. 

DAYFILE, JCLOUT. 
REPLACE, JCLOUT. 
REPLACE(OUTPUT=OUTR 
REPLACE(TAPE2=SVBP1 
DMD(200) 

/EOF 


UNI ) 
2 ) 


UNI) 

2 ) 


For card images operation: 


1 . 
2 . 


delete /JOB and /NOSEQ instructions (-first 
replace /EOF with a 6/7/S/9 multiple punch 


2 cards ) . 
car d . 


FIGURE 52. CDC JCL FOR EXECUTION OF 


THE BASE PROGRAM BATCH MODE. 
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FOLLOW NORMAL LOGON PROCEDURE 


GET.SVOABS. 

GET,TAPE1 = DATA1 . 

GET,TAPES=DATAS. 

GET,TAPE10=DATA10. 

GET,TAPE11 = DATA11 . 

GET,TAPE12=DATA12. 

SVOABS. 

NOTE 1; DATA FILES BELOW ARE USED AS NEEDED FOR EACH RUN. 

DATAl = INPUT DATA FILE. 

DATAS = TIME HISTORY RESTART DATA FILE. 

DATAIO = E927 MATRICES FILE. 

DATAll = GAOO ROTOR IMPEDANCES FILE. 

DATA12 = USER INPUT OF ROTOR IMPEDANCES FILE. 

. . . AFTER "?" APPEARS, KEY IN THE LETTER "T" AND CLEAR THE SCREEN. 

. . . OUTPUT RESULTS WILL BE PRINTED OUT IN BLOCKS OF 30 LINES. 

. . . A PAUSE ALLOWS THE USER TO MAKE A COPY IF DESIRED. 

. . . TO CONTINUE, CLEAR THE SCREEN AND KEY IN A (PERIOD). 

. . . AT THE END OF THE CALCULATIONS, KEY IN A TO END EXECUTION. 

. . . PLOT FILE RESULTS ARE ON TAPE2. 

. . . OUTPUT PRINTOUT IS ON TAPES. 

REPLACE, TAPE2=DATA2. <==z SEE NOTE 2 BELOW 

REPLACE, TAPE3=DATA3. <=== ” 

NOTE 2: 

DATA2 = PLOT DATA FILE. 

DATAS = OUTPUT DATA FILE. 

. . . TO OBTAIN A HARD COPY OF THE OUTPUT, REPLACE CARD ABOVE WITH: 

REWIND, TAPES. 

COPYSBF, TAPES, LIST. 

DELIVER. BLDG JOHN DOE 

ROUTE, LIST, DC=LP. 

. . . TO STACK UP PLOT FILES, REPLACE TAPE2 CARD ABOVE WITH: 

APPEND, DATA2,TAPE2, 

GET,DATA2. THIS PROCEDURE TAKES OUT THE 

PACK,DATA2. <=== -EOR- MARK AT THE END OF THE 

REPLACE, DATA2. PLOT DATA FILE ”DATA2”, 

. . . NOW THE PLOT DATA FOR 2 RUNS IS ON ONE FILE! 

FIGURE SZ. CDC JCL FOR EXECUTION OF THE BASE PROGRAM INTERACTIVE MODE. 



150 



151 


TABLE 1 


BASE PROGRAM SUBSTRUCTURES 


Component 

Flag 

Descriptive 

Name 

BFl 

BIFILAR 1 

BF2 

BIFILAR 2 

BMl 

BEAM MODEL 1 

CNl 

CONSTRAINT 1 

FAl 

FIXED ABSORBER 1 

GFl 

GENERALIZED FORCE 1 

ISl 

ISOLATOR 1 

MSI 

MODAL SUBSTRUCTURE 1 

RE2 

ROTOR ELASTIC 2 

RE3 

ROTOR ELASTIC 3 


component Description 


Inplane linear bifilar vibration absorber for reduction of 
inplane rotor vibrations 

Vertical linear bifilar vibration absorber for reduction of 
out-of-plane rotor vibrations 
Uniform elastic beam segment 

Component allowing elimination of various degrees of free- 
dom at a connection node 

Single mass, spring and damper vibration absorber for re- 
duction of aircraft vibrations at specific locations 
Force component allowing specification of harmonic force 
and moment input excitations 

Modal isolator employing two anti-resonant bar weights 

Normal mode representation of a dynamic structure 

Aeroelastic rotor model represented by mass, damping 
and stiffness matrices derived from program E927 (hover only) 
Aeroelastic rotor model represented by an impedance matrix 
derived from program G400 (hover and forward flight). 

Higher Harmonic Control (HHC) option available with RE3 . 



TABLE 2. BASE PROGRAM SOLUTION AND OPERATIONAL MODES 


Mode 

Flag 

Descriptive 

Name 

EGl 

EIGENSOLUTION 1 

EG2 

EIGENSOLUTION 2 

FRl 

FORCED RESPONSE 1 

THl 

TIME HISTORY 1 

GEN 

GENERAL 

PVl 

PARAMETER VARIATION 1 


Component Description 


Eigensolution mode for systems without damping (real) 

Eigensolution mode for systems with damping (complex) 

Forced response solution mode yielding harmonics of steady 
state response 

Time history solution mode using Newmark Beta Method 
(Reference 1) 

General control operational mode providing general control 
for a computer run 

Parametric variation operational mode facilitating variation 
of parameters for design studies 
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TABLE 3. SUMMARY OF ATTRIBUTES OF BASE PROGRAM COMPONENTS 


Number 

Name 

Number 
of Input 
Locations 
(NDLOAD) 

Number 
of Work 
Locations 
(NWORK) 

Number of 

Connection 

Nodes 

Number of 
Coordinates 

Output Coordinates Fortran Names 

1 

BFl 

14 

0 

1 

6 

CAMMAC , GAMMAS , X , Y , 2 , THTX . THTY . TUTZ 

2 

BF2 

16 

0 

1 

9 

BETAO , BETAC , BETAS , X , Y , 2 , THTX , THTY , THT2 

3 

BMl 

21 

0 

2 

12 

Xl,yi,21.THTXl.THTYl,THT21 (at end 1>, 

4 

CNl 

10 

0 

1 

6 

X2,y2,22,THTX2,THTY2,THTZ2 (at end 2) 
X , Y , 2 , THTX , THTY , THTZ 

5 

*EC1 

1 

0 

- 

- 

- 

6 

>*EG2 

1 

0 

- 

- 

- 

7 

FAl 

12 

0 

1 

7 

DELTA , X , Y , 2 , THTl , THT2 , THT3 

8 

*FR1 

3 

0 

- 

- 

- 

9 

+CEN 

10 

0 

- 

- 

- 

10 

GFl 

23 

0 

1 

6 

X , Y, 2 , THTl , THT2 , THT3 

11 

ISl 

73 

0 

2 

18 

XT, YT , 2T , XB , YB , 2B , XDT , YBT , 2BT, TXBT, TYBT , T2BT, 

12 

MSI 

57 

0 

1 to 5 

1 

XBF, YBF, 2BF, TXBF,TYBF, T2BF 
MODE 

13 

4pVl 

23 

0 

- 

- 

- 

14 

RE2 

5 

30+3*NC*NC 

1. 

1+NC 

£927 names (see section 7) - NC=29 max. 

15 

RE3 

19 

232 (ISTRSS=0) 

1 

6 

Displacements are X,Y,Z,WX,WY,WZ. 

16 

*TH1 

6 

3269 (ISTRSS=1) 
2 

- 

- 

Forces/moments are RX,RY,RZ,RMX,RMY,RMZ 

^Indicates solution mode. 
■^Indicates operational mode. 
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TABLE 4. 

SUMMARY OF BASE 

PROGRAM DEBUG OUTPUT OPTIONS 

Segment 

Subroutine 

Solution Methods 

Description of Debug Printout 

No. 

Name 

EGl 

EG2 FRl THl 



8 

RECOVR 




X 

Time history restart input data 

9 

TRNSFM 

X 

X 

X 

X 

Matrix coordinate transformations 

9 

RENTHl 




X 

Initial displacements and velocities 

9 

REITHI 




X 

As above 

12 

UPDMDl 




X 

Displacements and velocities at each time step 

13* 

ASMMCK 

X 

X 

X 

X 

Assembly of M,C, and K matrices 

15* 

S02EG1 

X 




Eigenvalues and eigenvectors - real 

15* 

S02EG2 


X 



Eigenvalues and eigenvectors - complex 

16* 

S02FR1 



X 


Forced response solution intermediate results 

16* 

S02TH1 




X 

Time history solution intermediate results 

16 

ASMZ2 



X 


G400 impedance matrix and force vector assembly 

16 

ASMH2 



X 


G400 H matrix relating hub reaction loads. to HHC 
inputs 

17 

RESFRl 



X 


HHC angle inputs used in the calculations of 
responses 

17 

HTRE3 



X 


Intermediate matrix P and optimal HHC angles 

17 

HXRE3 



X 


Dependent coordinate response - cosine and sine 
components 

17 

TDPRE3 



X 


T matrices in independent and dependent domains 
for HHC application 

17 

VDEP 



X 


Baseline vector in dependent domain without HHC 

17 

WZDGFl 



X 


Weighting factors for responses with HHC 
application 

17 

STRRE3 



X 


Baseline vector, HHC angles, stress vector and 
and one-half peak-to-peak stresses 

17 

RE2EG1 

X 




Displacements as referred to a local axis system 

17 

RE2FR1 



X 


As above 

17 

RE2TH1 




X 

As above 

18 

SAVTHl 




X 

Time history data saved for restart 

* Indicates that 

sample outputs are included 

in this report. 

X Indicates which 

solution 

method 

is 

applicable. 



TABLE 5. DESCRIPTION OF HIGH LEVEL ROUTINES 
Segment Leading 

Number Routine Segment Description 


1 

MAINSV 

Main program 

2 

INPCDC 

Controls subroutines reading, printing, and 
converting to internal program units of com- 
ponent input data 

3 

INPUT 

Reads the component input data 

4 

P INPUT 

Controls component dedicated routines reading 
component input data 

5 

PRINTl 

Prints out input for all components except ISl 
and RE3 

6 

PRINT2 

Prints out input for components ISl and RES 

7 

LINK 

Controls the reading of input data transmitted 
by external programs E927 and G400 

8 

RECOVR 

Recovers component initial displacements and 
velocities for a time history restart con- 
dition 

9 

TRNSFM 

Forms matrix transformations relating to 
system coordinates 

10 

LABEL 

Forms labels for output coordinates and solution 
methods 

11 

COMPUT 

Controls routines processing the data and solves 
the system equations 

12 

UPDATE 

Updates system parameters or value of time 

13 

ASSMBL 

Assembles system mass, damping, and stiffness 
matrices 

14 

SOLVE 

Controls choice of solution methods (eigen- 
solutions or forced response/time history) 

15 

SOLEGN 

Calculates eigensolutions (real or complex) 

16 

SOLFRT 

Calculates forced response or time history 
solution 

17 

RESTOR 

Restores system coordinate displacements from 
a global to a local axis system 

18 

SAVE 

Saves final results for printing, plotting, 
and for time history solution restart conditions 



TABLE 6. SEGMENTATION STRUCTURE AND COMMON BLOCKS 


Segment 

Fortran Subroutine 


COMMON Block 

Number 

Name 



Name 

1 

MAINSV 



ACONVR 


BLOCK 

DMTRX 


CBETA 


CINBFl 

DSUB 


CDATA 


CINBMl 

lASKSP 


CDEBUG 


CINFRl 

IRMVSP 


CFRl 


C INGEN 

ISETSP 


CGEN 


CINGFl 

MATMUL 


CGFl 


CINISl 

RECFBl 


CINOUT 


CINKLl 

RECCNl 


Cl TER 


CINMSl 

RECFAl 


CLABEL 


CINRE2 

RECISl 


CMAXA 


CINRE3 

RECKLl 


CNAMES 


CON INF 

RECMSl 


CPLOT 



RECRE2 


CRE2 





CRE3 





CRE3H1 





CRE3H2 ■ 





CRSTRT 

2 

INPCDC 



- 

3 

+ INPUT 





+FFLOAD 

+LOAD 


- ■ 


FREAD 

NSPACE 



4 

P INPUT 




5 

PRINTl 


. 



PINBFl 

PINKLl 


- 


PINBMl 

PINMSl 




PINCNl 

PINPVl 




PINEGl 

PINRE2 




PINFAl 

PINTHl 




PINFRl 





P INGEN 





PINGFl 




6 

PRINT2 





PINISl 

PINRE3 



7 

LINK 





LK2RE2 

PRTRE3 


- 


LK2RE3 

PRTSTR 



8 

RECOVR 





RECTHl 

- 


- 
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TABLE 6. CONTINUED. 


Segment FORTRAN Subroutine COMMON Block 

Number Name Name 


9 TRNSFM 


ASCND 

CRDRE3 

BET 

DEPCl 

BETO 

EULTRN 

BETl 

INDCl 

BET2 

INDC2 

COORD 

RENBFl 

CRDBFl 

RENCNl 

CRDBMl 

RENFAl 

CRDCNl 

RENISl 

CRDFAl 

RENKLl 

CRDGFl 

RENMSl 

CRDISl 

• RENRE2 

CRDKLl 

RENTHl 

CRDMSl 

REITHI 

CRDRE2 



10 


LABEL 



+ACON 

NAMBMl 

- 


+ACONDP 

NAMCNl 



+CORE 

+NAMES 



+LBLEG1 

NAMFAl 



+LBLEG2 

NAMGFl 



LBLFRl 

NAMISl 



LBLTHl 

NAMMSl 



NAMBFl 

NM2RE2 



NAMBF2 

NM2RE3 


11 


COMPUT 



*LINV3F 

*UERTST 



*LUDATF 

*UGETIO 



*LUELMF 



12 


UPDATE 



FFREAD 

UPDPVl 

- 


UPDMDl 

UPDTHl 


13 


ASSMBL 



ASMALL 

CPFAl 

- 


ASMMCK 

CPISl 



CP 

CPKLl 



CPBFl 

CPMSl 



CPBF2 

CP2RE2 



CPBMl 

CP2RE3 


14 


SOLVE 
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TABLE 6. CONCLUDED 


Segment FORTRAN Subroutine COMMON Block 

Number Name Name 


15 

*EIGZF 

*EQZQF 

*EQZTF 

*EQZYF 

SOLEGl 

SOLEG2 

SOLEGN 

S02EG1 

S02EG2 

*VHSH2C 

*VHSH2R 

*VHSH3R 

16 


SOLFRT 


ASMH2 

NMARK2 


ASMZ2 

SOLFRl 


FORGFl 

SOLTHl 


*LEQT2F 

S02FR1 


*LINV2F 

S02TH1 


LUREFF 

*VXADD 


MATDl 

*VXMUL 


NMARKl 

*VXSTO 

17 


RESTOR 


HBFRE3 

RE2EG1 


HJRE3 

RE2EG2 


HTRE3 

RE2FR1 


HXRE3 

RE2TH1 


PRTHHC 

STRRE3 


PSTRSS 

TDPRE3 


RESEGl 

VDEP 


RESEG2 

RESFRl 

WZDGFl 


RESTHl 


13 


SAVE 


+GVPLOT 

SAVTHl 


PRINT 

SV2EG1 


SAVEGl 

SV2EG2 


SAVEG2 

SV2FR1 


SAVFRl 

SV3TH1 

*IMSL routines to 

be replaced by the contractor. 

+Fortran routines 

to be revised or replaced for CDC system operation. 



TABLE 7 ALPHABETICAL LIST OF BASE PROGRAM FORTRAN ROUTINES 


No. 

Routine 

Name 

Seg. 

No. 

No. 

Routine 

Name 

Seg. 

No. 

No 

Routine 

Name 

Seg. 

No. 

1 

+AC0N 

10 

33 

CPFAl 

13 

65 

HXRE3 

17 

2 

+ACONDP 

10 

34 

CPISl 

13 

66 

lASKSP 

1 

3 

ASCND 

9 

35 

CPKLl 

13 

■ 67 

INDCl 

9 

4 

AS MALL 

13 

36 

CPMSl 

13 

68 

INDC2 

9 

5 

ASMH2 

16 

37 

CP2RE2 

13 

69 

INPCDC 

2 

6 

ASMMCK 

13 

38 

CP2RE3 

13 

70 

+ INPUT 

3 

7 

ASMZ2 

16 

39 

CRDBFl 

9 

71 

IRMVSP 

1 

a 

ASSHBL 

13 

40 

CRDBMl 

9 

72 

ISETSP 

1 

9 

BET 

9 

41 

CRDCNl 

9 

73 

LABEL 

10 

10 

BETO 

9 

42 

CRDFAl 

9 

74 

+LBLEG1 

10 

11 

BETl 

9 

43 

CRDGFl 

9 

75 

+LBLEG2 

10 

12 

BET2 

9 

44 

CRDISl 

9 

76 

LBLFRl 

10 

13 

BLOCK 

1 

45 

CRDKLl 

9 

77 

LBLTHl 

10 

14 

CINBFl 

1 

46 

CRDMSl 

9 • 

78 

*LEQT2F 

16 

15 

CINBMl 

1 

47 

CRDRE2 

9 

79 

LINK 

7 

16 

CINFAl 

1 

48 

CRDRE3 

9 

80 

*LINV2F 

16 

17 

CINFRl 

1 

49 

DEPCl 

9 

81 

*LINV3F 

• 11 

18 

C INGEN 

1 

50 

DMTRX 

1 

82 

LK2RE2 

7 

19 

CINGFl 

1 

51 

DSUB 

1 

83 

LK2RE3 

7 

20 

CINISl 

1 

52 

*EIGZF 

15 

84 

+LOAD 

3 

21 

CINKLl 

1 

53 

*EQZQF 

15 

85 

*LUDATF 

16 

22 

CINMSl 

1 

54 

*EQZTF 

15 

86 

*LUELMF 

16 

23 

CINRE2 

1 

55 

*EQZVF 

15 

87 

*LUREFF 

16 . 

24 

CINRE3 

1 

56 

EULTRN 

9 

88 

MAINSV 

1 

25 

COMPUT 

11 

57 

+FFLOAD 

3 

89 

MATDl 

16 

26 

CONINP 

1 

58 

FFREAD 

12 

90 

MATMUL 

1 

27 

COOED 

9 

59 

FORGFl 

16 

91 

NAMBFl 

10 

28 

+CORE 

10 

60 

FREAD 

3 

92 

NAMBF2 

10 

29 

CP 

13 

61 

+GVPLOT 

18 

93 

NAMBMl 

10 

30 

CPBFl 

13 

62 

HBFRE3 

17 

94 

NAMCNl 

10 

31 

CPBF2 

13 

63 

HJRE3 

17 

95 

+NAMES 

10 

32 

CPBMl 

13 

64 

HTRE3 

17 

96 

NAMFAl 

10 

*IMSL routines to be 
+Fortran routines to 

replaced by the contractor, 
be revised or replaced for 

CDC : 

system operation. 



TABLE 7 CONCLUDED. 


No. 

Routine 

Seg. 

No. 

Routine 

Seg. 

No. 

Routine 

Seg. 


Name 

No. 


Name 

No. 


Name 

No. 

97 

NAMGFl 

10 

129 

RECCNl 

1 

161 

SOLEGl 

15 

98 

NAMISl 

10 

130 

RECFAl 

1 

162 

SOLEG2 

15 

99 

NAMMSl 

10 

131 

RECISl 

1 

163 

SOLFRT 

16 

100 

NMARKl 

16 

132 

RECKLl 

1 

164 

SOLFRl 

16 

101 

NMARK2 

16 

133 

RECMSl 

1 

165 

SOLTHl 

16 

102 

NM2RE2 

10 

134 

RECOVR 

8 

166 

SOLVE 

■14 ■ 

103 

NM2RE3 

10 

135 

RECRE2 

1 

167 

S02EG1 

15 

104 

NSPACE 

3 

136 

RECTHl 

8 

168 

S02EG2 

15 

105 

PINBFl 

4 

137 

RENBFl 

9 

169 

S02FR1 

16 

106 

PINBMl 

4 

138 

RENCNl 

9 

170 

S02TH1 

16 

107 

PINCNl 

4 

139 

RENFAl 

9 

171 

STRRE3 

17 

108 

PINEGl 

4 

140 

RENISl 

9 

172 

SV2EG1 

18 

109 

PINFAl 

4 

141 

RENKLl • 

9 

173 

SV2EG2 

18 

110 

PINFRl 

4 

142 

RENMSl 

9 

174 

SV2FR1 

18 

111 

P INGEN 

4 

143 

RENRE2 

9 

175 

SV3TH1 

18 

112 

PINGFl 

4 

144 

RESTHl 

9 

176 

TDPRE3 

17 

113 

PINISl 

4 

145 

RESEGl- 

17 

177 

TRNSFH 

9 

114 

PINKLl 

4 

146 

RESEG2 

17 

178 

*UERTST 

. 14 

115 

PINMSl 

4 

147 

RESFRl 

17 

179 

*UGETIO 

14 

116 

P INPUT 

4 

148 

RESTHl 

17 

180 

UPDATE 

12 

117 

PINPVl 

4 

149 

RESTOR 

17 

181 

UPDMDl 

12 

118 

PINRE2 

4 

150 

REITHI 

9 

182 

UPDPVl 

12 

119 

PINRE3 

4 

151 

RE2EG1 

17 

183 

UPDTHl 

•12 

120 

PINTHl 

4 

152 

RE2EG2 

17 

184 

VDEP 

17 

121 

PRINT 

18 

.153 

RE2FR1 

■ 17 

185 

*VHSH2C 

15 

122 

PRINTl 

5 

154 

RE2TH1 

17 

186 

*VHSH2R 

15 ■ 

123 

PRINT2 

6 

155 

SAVE 

18 

187 

*VHSH3R 

15 

124 

PRTHHC 

17 

156 

SAVEGl 

18 

188 

*VXADD 

16 

125 

PRTRE3 

7 

157 

SAVEG2 

18 

189 

*VXMUL 

16 

126 

PRTSTR 

17 

158 

SAVFRl 

18 

190 

*VXSTO 

16 

127 

PSTRSS 

17 

159 

SAVTHl 

18 

191 

WZDGFl 

17 

128 

RECBFl 

1 

160 

SOLEGN 

15 




*ISML routines 

to be 

replaced 

by the 

contractor, 




+Fortran routines to 

be revised or replaced for 

CDC system operation. 



TABLE 8. FUNCTIONS OF FORTRAN ROUTINES IN BASE PROGRAM 


Routine 

No. Name Brief Description 


+ 1 

ACON 

IBM routine, called by ACONDP, which converts 
floating point or integer numbers to alpha 
numeric characters. Not used for CDC appli- 
cations . 

+ 2 

ACONDP 

IBM routine to be replaced by equivalent CDC 
routine ENCODE. It is called by subroutines 
NAMES, LBLEGl, and LBLEG2 . It provides inter- 
face between a global double precision program 
module and single precision system routines. 

3 

ASCND 

Sorts a vector of integers into ascending order 
Ex: 3, 3, 3, 2, 1,1, 4, 4^ to 1,1, 2, 3, 3, 3, 4, 4. 

4 

ASMALL 

Sets computer storage needed for assembling 
mass (M), damping (C), and stiffness (K) 
matrices and force vector. 

5 

ASMMCK 

Assembles the global M,C, and K matrices from 
the local matrices and forms the global force 
vector. 

6 

ASMH2 

Assembles the rotor impedance matrix H. 

7 

ASMZ2 

Assembles the rotor impedance matrix (Z) and 
the force vector (F) for component RE3 and 
forced response solution FRl . 

8 

ASSMBL 

Controls assembly of M,C, and K matrices for 
the different solution methods. 

Q 

BET 

Saves non-zero elements of the BETA matrix 
(product of BETA), BETAl , and BETA2 below) 
in vector form for fast assembly. 

10 

BETO 

Forms BETAO transformation matrix relating 
local element coordinates to the global element 
coordinates. 

11 

BETl 

Forms BETAl transformation matrix relating de- 
pendent and independent coordinates in the 
global system. 

12 

BET2 

Forms BETA2 transformatin matrix relating 
coordinates defined by routines INDCl and 
INDC2. 

13 

BLOCK 

Initializes data in COMMON blocks. 

NOTE : Subroutines CINBFl (no. 13) through 

CINRE3 (no. 23) perform the same func- 
tion, which is to convert component 
input data to internal program units 
(slug, foot, sec, and radian). The 
component name is identified by the 
last 3 characters in the routine's 
Fortran name. 

14 

CINBFl 

See note above (this routine handles both in- 
plane and vertical bifilars). 

15 

CINBMl 

See note above . 
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TABLE 8. CONTINUED. 


Routine 

No. Name Brief Description 


16 

CINFAl 

See note above . 

17 

CINFRl 

See note above. 

18 

C INGEN 

See note above. 

19 

CINGFl 

See note above. 

20 

CINISl 

See note above. 

21 . 

CINKLl 

See note above. 

22 

CINMSl 

See note above. 

23 

CINRE2 

See note above . 

24 

CINRE3 

See note above. 

25 

COMPUT 

Controls lower level routines processing the 
data and solving system equations of motion at 
each step (parametric variation or time step). 

26 

CONINP 

Controls routines converting input units to in- 
ternal units. 

27 

COORD 

Defines the number of coordinates to be used 
for each component. 

+28 

CORE 

IBM routine, called by ACON, which provides 
computer storage. Not used or CDC applications 

29 

• CP 

Controls component dedicated routines forming 


local M,C, and K matrices and local force 
vector. 


NOTE : Subroutines CPBFl (no. 29) through CP2RE3 
(no. 37) perform the same function, which 
is to calculate the local M,C, and K 
matrices and local force vector for a system 
component. The component name is identi- 
fied by the last 3 characters in the rou- 
tine name . 


30 

CPBFl 

See 

note 

above 

31 

CPBF2 

See 

note 

above 

32 

CPBMl 

See 

note 

above 

33 

CPFAl 

See 

note 

above 

34 

CPISl 

See 

note 

above 

35 

CPKLl 

See 

note 

above 

36 

CPMSl 

See 

note 

above 

37 

CP2RE2 

See 

note 

above 

38 

CP2RE3 

See 

note 

above 



Routine 

Name 


TABLE 8. CONTINUED. 


No. 


Brief Description 


NOTE : Subroutines CRDBFl (no. 38) through CRDRE3 
(no. 47) perform the same function, which 
is to specify the number of coordinates, 
connection node numbers, and Euler angles' 
for each component. The component name 
is identified by the last 3 characters in 
the routine name . 


39 

CRDBFl 

See note above (this routine handles both in- 
plane and vertical bifilars). 

40 

CRDBMl 

See note above . 

41 

CRDCNl 

See note above. 

42 

CRDFAl 

See note above . 

43 

CRDGFl 

See note above . 

44 

CRDISl 

See note above . 

45 

CRDKLl 

See note above. 

46 

CRDMSl 

See note above . 

47 

CRDRE2 

See note above . 

48 

CRDRE3 

See note above. 

49 

DEPCl 

Forms attributes of the dependent coordinates 
of the assembly. 

50 

DMTRX 

Displays non-zero elements of a matrix or a 
vector. 

51 

DSUB 

Writes out the name of the subroutine being 
called for debug printout 

* 52 

EIGZF 

IMSL routine called by S02EG1 and S02EG2 for 
eigensolutions . It calculates eigenvalues 
and eigenvectors (optionally) of the system 
A*X= *B*X where A and B are real 

matrices. 

* 53 

EQZQF 

IMSL routine called by EIGZF for eigensolu- 
tions. It reduces two input matrices A and 
B, simultaneously, A to upper Hessenberg and ■ 
B to upper triangular form. 

* 54 

EQZTF 

IMSL routine callee by EIGZF for eigensolu- 


tions . It reduces an upper Hessenberg matrix 
A to quasi-upper triangular form while keep- 
ing matrix B triangular. 



TABLE 8. CONTINUED. 

Routine 

No. Name Brief Description 

* 55 EQZVF IMSL routine called by EIGZF for eigensolu- 

tions . It calculates the eigenvalues and 
eigenvectors (optionally) of the system A * Z 
= * B * Z where A is quasi-upper triangular 

and B is upper triangular. 


56 

EULTRN 

Forms 3x3 rotation transformation matrix 
of direction cosines from input Euler angles 
(THETA, PHI, XSI ) . 

+ 57 

FFLOAD 

Reads component name and element number, prints 
out run heading and image of input data, con- 
trols reading of input data for successive 
cases and stores input data into 2 vectors, 
ILOAD for integer values and DLOAD for real 
values. The WRITE statements are computer . 
dependent (CDC or IBM). 

58 

FFREAD 

Updates input data; called by UPDATE in the 
execution of a parametric variation study. 

59 

FORGFl 

Forms the force vector from the input excita- 
tion forces for component GFl . 

60 

FREAD 

Updates input data called by INPUT at the 
start of the run. This routine is identical 
to FFREAD; a duplicate routine is used to re- 
duce the size of the first overlay segment. 

+ 61 

GVPLOT 

Writes titles, output results and labels to a 
file on unit 2 for plotting purposes. The 
WRITE statements are computer dependent (CDC 
or IBM). 

62 

HBFRE3 

Restores hub interface forces and moments 
resolved to the hub local axis from a solution 
for rotor type RE3 . 

63 

HJRE3 

Computes performance index J for HHC application 

64 

HTRE3 

Computes HHC optimum control angles . 

65 

HXRE3 

Computes dependent coordinate response with HHC. 

66 

lASKSP 

FUNCTION lASKSP returns information defining 
computer storage allocations for a specific 


variable. 

EX: LOCEXP=IASKSP ( lA, A, 4HEXPL , KIND, KLASS , 

NROWS , NCOLS ) . The function subroutine lASKSP 
provides the storage location, LOCEXP, for 
variable EXPL in the lA vector which contains 
integer values (the A vector contains real 
values) and the storage locations for the 
attributes KIND, KLASS, NROWS, NCOLS, where 

KIND = 1 for real, 2 for complex values 
KLASS = 1 for all cases 
NROWS = number of rows in EXPL 
NCOLS = number of columns in EXPL 
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f 


TABLE 8. CONTINUED. 

Routine 

No. Name Brief Description 

If the variable name EXPL is not found in the 
catalog entry list of the lA vector, a 
message is printed to that effect and pro- 
gram execution is terminated. The programmer 
should then use the function ISETSP'to reserve 
storage for the variable EXPL (see writeup for 
ISETSP, no. 68). 

67 INDCl Stores the element number, connection node 

number, displacements and velocities of the 
independent coordinates of the assembled 
system. 

68 INDC2 Calculates the size of the independent coor- 

dinates vector and assigns global element 
numbers and node numbers to the set of inde- 
pedent coordinates by removing those coordinates 
expressible in terms of modal coordinates for 
elastic structures. 

69 INPCDC Controls reading, printing, and converting to 

internal units of the component input data. 

+ 70 INPUT Reads the component input data. It also reads 

and writes up to 4 title cards. CDC format 
is 4A10, IBM format is 10A4. 

71 IRMVSP Frees computer storage allocations no longer 

needed for the variable specified. 

EX: CALL IRMVSP ( lA, A, 4HEXPL ) . The catalog 

entry list is first searched to find the 
variable name EXPL. If found, then storage 
is freed and the variable EXPL does not exist 
any longer. If the name EXPL is not found, 
than space was not previously allocated and 
therefore none can be freed; program yields 
a message to that effect and proceeds with 
execution. 

72 ISETSP FUNCTION ISETSP reserves computer storage for 

the variable specified. 

EX: LOCEXP=ISETSP (IA,A,4H EXPL, KIND, KLASS, 

NROWS, NCOLS ) . The routine first checks to see 
if storage has already been reserved for 
variable EXPL. If it has, it returns the 
value LOCEXP (location in the lA vector) 
and program execution continues. If it hasn't, 
it reserves storage locations in the lA and 
A vectors if enough space is available. If 
the space requested is not available, a 
message is printed to that effect and pro- 
gram execution is stopped. The programmer 
should then either reduce the storage re- 
quested for EXP or increase the size of the 
lA and A allocation vectors. 
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TABLE 8. CONTINUED. 


No. 

Routine 

Name 

Brief Description 

73 

LABEL 

Forms labels for output system coordinates 
and methods of solution. 

+ 74 

LBLEGl 

Forms the labels for eigenvalue solution EGl 
with no system damping. CDC uses ENCODE, IBM 
uses ACONDP. 

+ 75 

LBLEG2 . 

Forms the labels for eigenvalue solution EG2 
with system damping. CDC uses ENCODE, IBM 
uses ACONDP . 

76 

LBLFRl 

Forms the labels for the forced response 
solution FRl. 

77 

LBLTHl 

Forms the labels for the time history solution 
THl. 

* 78 

LEQT2F 

IMSL routine called by IMSL routine LINV2F and 
by the Base Program subroutines NMARKl and 
NMARK2 . It solves a set of linear equations 
(highly accurate). 

79 

LINK 

Controls the reading of input data trans- 
mitted by external programs E927 and G400. 

* 80 

LINV2F 

IMSL routine called by subroutine NMARKl . 

It calculates the inverse of a matrix (highly 
accurate - uses full storage). 

* 81 

LINV3F 

IMSL routine called by S02FR1 and HTRE3 . 

It calculates the inverse of a matrix (not as 
accurate as LINV2F but requires less storage). 

82 

LK2RE2 

Reads rotor matrices (M,C, and K) from program 
E927 stored in unit 10. 

83 

LK2RE3 

Reads rotor impedance matrix and force vector 
from program G400 stored in unit 11 or as 
specified by the user on unit 5 (standard 
read unit). 

+ 84 

LOAD 

Reads input data. The READ statements are 
computer dependent (CDC or IBM). 

* 85 

LUDATF 

IMSL routine called by IMSL routine LEQT2F. 
It performs L-U decomposition by the CROUT 
algorithm with optional accuracy test. 

* 86 

LUELMF 

IMSL routine called by IMSL routines LEQT2F 
and LUREFF. It eliminates part of the solu- 
tion of AX = B. 

* 87 

LUREFF 

IMSL routine called by IMSL routine LEQT2F. 
It refines the solution of a set of linear 
equations . 

+ 88 

MAINSV 

Controls the Base Program logic flow, sets 


size of the lA and A vectors, initializes 
both vectors to zero, and sets the first 5 
locations of the integer vector lA, At the 
end of the computer run, MAINSV prints out 
the maximum size of the working storage. 

For CDC code, add PROGRAM SIMVIB (...) card. 

MATDl Solves for the displacement response of the 

system using the Newmark Beta method to 
find the time history solution. 
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TABLE 8. CONTINUED. 


Routine 

No. Name Brief Description 


90 MATMUL Multiplies two matrices together. 

NOTE : Subroutines NAMBFl (no. 86) through 
NAMCNl (no. 89) all perform the same 
function, which is to assign names to 
the different degrees df freedom of 
each system component. 


91 

NAMBFl 

See note above . 

92 

NAMBF2 

See note above. 

93 

NAMBMl 

See note above. 

94 

NAMCNl 

See note above. 

+ 95 

NAMES 

Stores element number for each component in its 
coordinate labels. CDC uses ENCODE, IBM- 
uses ACONDP . 

96 

NAMFAl 

See note above before NAMBFl. 

97 

NAMGFl 

See note above before NAMBFl . 

98 

NAMISL 

See note above before NAMBFl . 

99 

NAMMSl 

See note above before NAMBFl. 

100 

NMARKl 

Calculates parameters required to solve the time 
history equvations of motion by the Newmark 
Beta Method for initial points. 

101 

NMARK2 

As above but for intermediate points. 

102 

NM2RE2 

See note above before NAMBFl . 

103 

NM2RE3 

See note above before NAMBFl . . 

104 

NSPACE 

Defines number' of input loader locations re- 
quired for system components . 

NOTE : Subroutines PINBFl (no. 100) through PINMSl 
(no. 110) perform the same functions, 
which is to print the input for the speci- 
fied component. 

105 

PINBFl 

See note above (this routine handles both in- 
plane and vertical bifilars). 

106 

PINBMl 

See note above. 

107 

PINCNl 

See note above. 

108 

PINEGl 

See note above. 

109 

PINFAl 

See note above. 

110 

PINFRl 

See note above. 

111 

P INGEN 

See note above . 

112 

PINGFl 

See note above. 

113 

PINISl 

See note above. 

114 

PINKLl 

See note above. 



TABLE 8. CONTINUED. 


No. 

Routine 

Name 

Brief Description 

115 

PINMSl 

See note above . 

116 

P INPUT 

Controls the routines which print the com- 
ponent input data. 

117 

PINPVl 

See note above before PINBFl. 

118 

PINRE2 

See note above before PINBFl. 

119 

PINRE3 

See note above before PINBFl. 

120 

PINTHl 

See note above before PINBFl. 

121 

PRINT 

Prints output results. 

122 

PR INTI 

Controls the printing of component input data 
except for I SI and RE2 . 

123 

PRINT2 

Controls the printing of component input data 
for ISl and RE3 . 

124 

PRTHHC 

Prints summary of HHC results . 

125 

PRTRE3 

Prints rotor impedance matrix and force vector 
from G400 program as stored on unit 11 or from 
user input on unit 12. 

126 

PRTSTR 

Prints stress matrices from G400 program. 

127 

PSTRSS 

Prints summary of stresses. 



NOTE: Subroutine RECBFl (no. 118) through RECMSl 
(no. 123) perform the same function, which 
is to recover initial values of displace- 
ments and velocities for time history 
restart conditions . 

128 

RECBFl 

See note above (this routine handles both in- 
plane and vertical bifilars). 

129 

RECNl 

See note above. 

130 

RECFAl 

See note above . 

131 

RECISl 

See note above . 

132 

RECKLl 

See note above . . 

133 

RECMSl 

See note above . 

134 

RECOVR 

Controls the recover of initial values of com- 
ponent displacements and velocities for a time 
history restart condition. 

135 

RECRE2 . 

See note above before RECBFl. 

136 

RECTHl 

See note above before RECBFl . 



TABLE 8. CONTINUED. 

Routine 

No. Name Brief Description 

NOTE : subroutine RENBFl (no. 127) through RENRE2 
(no. 133) perform the same function, which is 
to initialize the values of displacements and 
velocities in the globaldependent coordinate 
system and/or recover the component coordinates 
number. 


137 

RENBFl 

See note above. 

138 

RENCNl 

See note above. 

139 

RENFAl 

See note above. 

140 

RENISl 

See note above. 

141 

RENKLl 

See note above. 

142 

RENMSl 

See note above. 

143 

RENRE2 

See note above. 

144 

RENTHl 

Controls the recovery of the initial displace- 
ments and velocities for the various components. 

145 

RESEGl 

Controls restoring of component displacements 
from the global to the local axis system for 
eigensolution EGl . 

146 

RESEG2 

As above but for eigensolution EG2 . 

147 

RESFRl 

As above but for forced response solution FRl . 
In addition, it controls the restoring of the 
hub interface forces for component RE3 . 

148 

RESTHl 

As above for RESEGl but for the time history 
solution THl. 

149 

RESTOR 

Controls the restoring of component displace- 
ments from the global to the local axis system 
for the four methods of solution ( EGl , EG2 , FRl , 
and THl ) . 

150 

REITHI 

Resolves component displacements, velocities, 
and accelerations from local to global axis 
system. 

151 

RE2EG1 

Restores component displacements from the global 
to the local axis system for eigensolution EGl. 

152 

RE2EG2 

As above but for eigensolution EG2 . 

153 

RE2FR1 

As above but for forced response solution FRl. 

154 

RE2TH1 

As above but for time history solution THl . 

155 

SAVE 

Controls the saving of output results for 
printing and/or plotting. 

156 

SAVEGl 

Controls saving of output results for eigen- 
solution EGl. 

157 

SAVEG2 

As above but for eigensolution EG2 . 



TABLE 8. CONTINUED. 


No. 

Routine 

Name 

Brief Description 

158 

SAVFRl 

As above but for forced response solution FRl. 

159 

SAVTKl 

Saves output results for time history solution 
THl . 

160 

SOLEGN 

Controls choice o.f eigensolution EGl for systems 
without damping or eigensolution EG2 for systems 
with damping. 

161 

SOLEGl 

Controls eigenvalue solution EGl. 

162 

SOLEG2 

Controls eigenvalue solution EG2 . 

163 

SOLFRT 

Controls choice of forced response solution 
FRl or time history solution THl. 

164 

SOLFRl 

Controls forced response solution FRl. 

165 

SOLTHl 

Controls time history solution THl. 

166- 

SOLVE 

Controls type of solution requested: eigen- 

solution (EGl or EG2 ) or forced response/ 
time history (FRl or THl). 

167 

S02EG1 

Calculates the eigenvalues and eigenvectors 
for the system defined by the mass and stiff- 
ness matrices . 

168 

S02EG2 

Calculates the eigenvalues and eigenvectors 
for the system defined by the mass, damping, 
and stiffness matrices. 

169 

S02FR1 

Calculates the forced response solution. 
Assembles rotor impedance matrix for component 
RE3 with impedances from other components. 

170 

S02TH1 

Controls the calculations of the time history 
response of the system using the Newmark Beta 
method to find the solution. 

171 

STRRE3 

Calculates harmonics of stresses and pushrod 
loads and one-half peak-to-peak stresses and 
loads. 

172 

SV2EG1 

Saves output results for eigensolution EGl. 

173 

SV2EG2 

Saves output results for eigensolution EG2 . 

174 

SV2FR1 

Saves output for forced response solution FRl. 

175 

SV3TH1 

Saves output results for time history restart 
cases on unit 8. 

176 

TDPRE3 

Assembles T matrix in the dependent coordinates 
domain. 

177 

TRNSFM 

Forms the matrix transformations relating the 
system coordinates. 

*178 

UERTST 

IMSL routine called by IMSL routines EIG2F, 


EQZTF, LINV2F, LEQT2F, LUDATF, and LUREFF. 
It generates error messages. 



TABLE 8 . CONCLUDED . 

Routine 


No. 

Name 

Brief Description 

*179 

UGETIO 

IMSL routine called by IMSL routine UERTST. 

180 

UPDATE 

Updates time or parameter dependent coordinates 

181 

UPDMDl 

Defines displacements, velocities, and accel- 
erations for the independent coordinates at 
the next time increment. 

182 

UPDPVl 

Updates the independent parameter requested 
by component PVl . 

183 

UPDTHl 

Updates the time for a time history solution. . • 

184 

VDEP 

Calculates dependent coordinates for HHC and/, 
or stress options for the RE3 component. 

*185 

VHSH2C 

IMSL routine called by IMSL routine EQZVF. 
Zeroes out a single complex element of a 
matrix. 

*186 

VHSH2R 

IMSL routine called by IMSL routines EQZQF, 
EQZTF, and EQZVF. 

*187 

VHSH3R 

IMSL routine called- by IMSL routine EQZTF. 
Zeroes out 2 elements of a matrix. 

*188 

VXADD 

IMSL routine called by IMSL routines LUREFF 
and VXMUL. It adds a double precision float- 
ing point number to the extended precision 
accumulator. 

*189 

VXMUL 

IMSL routine called by IMSL routine LUREFF. 
It multiplies two double precision floating 
point numbers and adds the extended pre- 
cision result to the extended precision 
accumulator. 

*190 

VXSTO 

IMSL routine called by IMSL routine LUREFF. 

It moves a double precision floating point 
number contained in the most significant half 
of the extended precision accumulator into 
the location specified in the VXSTO argument 
list. 

191 

WZDGFl 

Recovers weighting factors for HHC option 
from the A storage vector. 


*ISML routines to be replaced by the contractor. 

+Fortran routines to be revised or replaced for CDC system operation. 
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TABLE 9- BASE PKOCRAM COMMON B1X)CKS/SUBR0UTINES CKOSS-UEFEKENCE LIST 

4 ♦ 4 


1 2 3 4 5 6 7 a 9 10 11 12 13 14 IS 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 


Subroutine 

Ndfne/ 

Common 

Block 




<1 §u V. 

'j''' VI Vic^w vO" ui 


o f! P S S £ ?: s:'| I S4‘ sS S leg 0 , gg 
&FAk ,3 ?. S S =1 2 2 Sll s & So- 


UI 


pj Subroutine 
H ;/ Name/ 

! Common 
Block 


1 

ACONVR 









ACONVR 

1 

2 

a:iTA 

X 








CliETA 

2 

3 

CDATA 

X 


X 

X 


X 



n)ATA 

3 

4 

CUELUG 

X 

X 

X X 

X 

X 

X 



CDEBUG 

4 

5 

crki 




X 

X 




CFRl 

5 

6 

CCEN 




X 

X 




CCEN 

6 

7 

cc;fi 




X 



X 


CGFl 

7 

8 

CINOIJT 




X 





CINOUT 

0 

9 

Cl TEH 






X 



CITER 

9 

10 

C1.AL3EL 









CLABEL 

10 

IJ 

CMAXA 




X 





CHAXA 

11 

12 

CNAHES 



X 

X 


X 

X 

X 

CNAMES 

12 

13 

CI'LOT 




X 





CFLOT 

13 

14 

CKE2 




X 




X 

CKE2 

14 

15 

CHE3 


X 

X 

X 


X 


X 

CRE3 

15 

10 

CHE3H1 




X 





CRE3H1 

16 

17 

CKE3H2 









CRE3H2 

17 

18 

ensTRT 




X 


X 



CRSTRT 

18 


NOTE; Denotes computer dependent subroutines (IBM or CDC). 


* Denotes IMSL routines which must be replaced by the Contractor. 
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TABLE 9- CONTINUED. 


32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 




. . I -rr O O , L- Subroutine 

“i Name/ 

' * fc ^1 P O UJ o ^ .J Common 
tj) Kl ^ □ -i n: g U) ft- Block 

VliluJUJllluJl^U-lLLicB 


1 

ACONVR 


X 




ACONVR 

1 

2 

CBLTA 






CBETA ' 

2 

3 

CDATA 


X 

X 

X 


CDATA 

3 

4 

COLDUC 

X 





CDEBUG 

4 

5 

CfRl 






CFRl 

5 

6 

CGEN 





X 

CGEN 

6 

7 

CGFl 






CCFl 

7 

8 

CINOUT 

X X 

X 



X 

CINOUT 

6 

9 

Cl TER 






CITER 

9 

10 

CLADEL 





X 

CLABEL 

10 

11 

CNAXA 


X 

X 

X 


CMAXa\ 

11 

12 

CNAMES 






CNAMES 

12 

13 

CPLOT 





X 

CPLOT 

13 

14 

CRE2 X 






CRE2 

14 

15 

CKE3 






CRE3 

15 

16 

CRE3H1 






CRE3H1 

16 

17 

CRE3H2 






CRC3I12 

17 

18 

CRSTRT 


• X 



X 

CRSTRT 

18 
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TABLE 9. CONTINUED. 


62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 76 79 80 01 62 63 84 65 66 67 88 89 90 91 92 93 94 95 


Subroutine 

Name/ UJ V/' 

Coinnon X* i'* 

Block Vi. U- 2; 

ei n >7. p 

i i: ^ — 


'p 

\v 

r 

u- 

H 

-J 

D 


a. 


;y -5 

P 

s: 



lU i] □ s: CJ UJ (J ^Common 

^ I- x; X.I s i: ^.Biock 

Dd.c£ 


1 

ACONVR 










X 


ACONVR 

1 

2 

CBETA 







XXX 

X 




CBLTA 

2 

3 

COATA 




X 



X 


X 



CDATA 

3 

4 

CUEBUG 

X 

X X 








X 


CDECUG 

4 

e 

CFRl 


X 










CFRl 

5 

6 

CGEN 




X X 








CGEN 

6 

7 

ccn 












CG FI 

7 

8 

CINOUT 



X 

X 

X 

X 




X 

X 

CINOUT 

8 

9 

Cl TEH 












CITER 

9 

10 

CIJ\BEL 




X 








CLABEL 

10 

11 

CMAXA 



X 

X 

X 

X 





X 

CMAXA 

11 

12 

CNAMES 




X 



X 


X 



CNAMES 

12 

13 

CP LOT 












CPI.OT 

13 

14 

CRE2 







X 


X 

X 


CRE2 

14 

15 

CRE3 









X 

X 


CRE3 

15 

16 

CRE3H1 

X 








X 

X 


CRE3H1 

16 

17 

CRE3H2 










X 


CRE3H2 

17 

18 

CRSTRT 











X 

CRSTRT 

18 



TABLE 9. CONTINUED. 


96 97 

98 

99 

100 

101 

102 

103 

104 

lOS 

106 

107 

108 

109 

110 

Ill 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

Subroutine] 

Name/ 

Common 

< U- 
U- 

\n 

<-» 

i: 

ii 

or 

<"d 

e 

n! 

u> 

rO 

I 

I’l: 

v*» 

K) 

< » 

u- 

V 

£• 

< » 
c 1 

! 

• ti ' 

a 

u. 

—4 

f . 

0 

tl. 

vO 

t4 

V-) 

:? 

vr 

Z3 

z: 

o_ 

i 

-cS 

D- 

f ^ 

d 

i 

H 

"> 


(V 

0 Subroutine 
CE Name/ 
Common 

Block 

r 

a. ^ 
2 ^ 

X 

<c 

2 

<:L 

2. 

cC 

■SL 

2 

CL 

5 

2 : 

Zl 

rr 

C- 

2 

Cl. 

"> 

(S- 

dl 

qT 

">■ 

u. 

• > 

2 

a. 


a- 

C 

2 

0- 

2 

a 

2 ! 

a,' 

5 

D- 

5; 

0 - 

(V 

D. 

t-- Block 
d 

a 


1 ACOMVR 






ACONVR 

1 

2 CDFTA 






CBETA 

2 

3 CPATA 



X 



CDATA 

3 

4 cde:bug 






CDLBUG 

4 

5 CFIU 

X 





CFRl 

5 

6 CCEJN 






CGEN 

6 

7 CCFl 






CGFl 

7 

8 CINOUT 



X 


X 

CINOUT 

8 

9 CITER 






CITER 

9 

10 CLABEL 






CLABEL 

10 

11 CM/\XA 






CMAXA 

11 

12 CNAJIES 


X 

X 


X 

X CNAMES 

12 

13 CPLOT 






CFLOT 

13 

14 CRE2 


X 




CRE2 

14 

15 CRE3 

X 

X 




CRE3 

15 

16 CKE3H1 


X 


X 


CRE3H1 

16 

17 CRE3H2 






CRE3H2 

17 

18 CRSTRT 






CRSTRT 

18 


M 

U1 
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TABLE 9. CONTINUED. 


SvjU- 

zoutiiie 

Ndiiie/ 

Coftuncn 

Block 


125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 140 

141 

1 1 

fo 

ixi 

»- 

</) 

1 

iX 

1 

u. 

CO 

(J 

2 

tH 

U- 

*J 

T » 

< ) 

V* 

_i 

u. 

u 

Vj 

VI 

> 

a 

<J 

OJ 

Is: 

11 

t 

x-4 

h 

a 

7> 

'•1 

•> 

V> 

2> 

vA 

▼4 

5 

»-4 

2 

V- 



•b 


i'J 

'U 

ii» 


iD 

•ii 




'.U 

iU 

•j) 

lA- 

; 

. 1 -. 

v^ 

a; 


lU 

q: 


u: 

Gi 


tK 

0: 

u: 



Cil 

L 

(V 


to 

^1$ 



n: 




r-» 

x: 

Ul 

DC 


O 

LU 

ti) 

Uj 


:r. 

h 

V"} 



d) 

10 

vO 

LU 

VO 

IU 

i/) 

IU 

lO 

LU 


>-t 

i'* 

0: 

fr 

rw 

u: 

CkC 

0: 


C/' 





] 

ACONVR 






2 

CBFJTA 



X 

X 

X 

3 

CDATA 

X 

X 

X 

x‘ 

X 

4 

CM-ibUC 

X 

X 



X 

5 

CFRl 





X 

C 

CGEN 






7 

CGFJ 





X 

8 

CINOUT 

X 




X 

9 

LITER 






10 

CU^BEL 






J 1 

CriAXA 



X 

X 

X 

12 

CNAMES 

X ' 

X 




13 

Cl-LOT 






14 

CRE2 






15 

CKE3 





X 

16 

CRE3H1 X 





X 

17 

CRE3II2 





X 

16 

CRSTRT 







Glib- 

routine 

Nsme/ 

Conimon 

Block 


ACONVR 

1 

CBETA 

2 

CL'ATA 

3 

CUFBflG 

4 

CFRl 

S 

cgi:n 

6 

CGFl 

7 

CINOUT 

8 

Cl TER 

9 

CLAtiEL 

10 

CMAXA 

11 

CNAMES 

12 

CPLOT 

13 

CRE2 

14 

CRE3 

15 

CRE3H1 

16 

CRE3H2 

17 

CRSTRT 

18 
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TABLE 9. CONTINUED 


Subroutine 

Nome/ 

Coiomon 

Block 


152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

Ot 





rt 




1 



"TJ 

"i 

\u! 

-ji 

o' 


rf{ 





T6 

-H 


O 

ID 

rt 

ra 

lu 

3:1 

Uj' 

\ai 

O 

ViJ 

> 

<C 

o 

<1. 

U. 

< 

1 

(,!!)! 

i 

VU 

tu 

vj 

a 

el 

o' 

lO 

Ui 

n:. 

VO’ 

Vu 

S 

r» 

a 

H 

ri 

-o 

Uj 

2 

o 

e 

VL 

r-t 

N 

X 

•5>; 

w 


Cd 





■ ^ 


v>. 


'A 

v-J 



J2i 

Jl 

tr 




\T) 

'■) 



1^ cn loi 
^ ' ’ 


Subroutine 
Name/ 


W ^ Common 
— Block 


'll; 


1 ACONVR 

2 CBETA 

3 CDATA 

4 CDE13UG 

5 CFHl 

6 CCEN 

7 CGFl 

8 CINOUT 

9 CITER 

10 CLABEL 

11 CMAXA 

12 CNAMES 

13 CPLOT 

14 CRE2 

15 CRE3 

16 CRE3H1 

17 CRE3H2 

18 CRSTRT 



X 


X 





ACONVR 

1 

CEETA 

2 

CDATA 

3 

CDEEUG 

4 

CFRl 

5 

CGIJN 

e 

CGFl 

7 

CINOUT 

6 

CITER 

9 

CLADEL 

10 

C^L^:<A 

11 

CNAflES 

12 

CTLOT 

13 

CHE2 

14 

CRE3 

15 

CRE3H1 

16. 

CRE3H2 

17 

CRSTRT 

18 
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TABLE 9. 


* 


* * * A 





Subroutine 

Name/ 

Common 

Block 


ACONVR 1 
CEETA 2 
COATA 3 
COHDUG 4 
CKRi 5 
CGEN 6 
CCFl 7 
CIIJOUT 8 
Cl TER 9 
CLAEEL 10 
CMAXA 1 1 
CNAMES 12 
CPLOT 13 
CRE2 14 
CRES 15 
CRE3H1 16 
CRE3H2 17 
CRSTRT 18 



TABLE 10, 


BASE PROGRAM INPUT/OUTPUT DATA FILES 


File 

Number 

Fortran 

Name* 

Subroutine(s ) for 
Referencing File 
(Segment No. of 
Subroutine) 

& 

Input/Output Function of File 
and Character of. Data 

1 

KINPUT 

1. INPUT (3) 

Read title cards (up to 4). 



2. LOAD (3) 

Read component input data 
(standard read file). 

2 

KSAVE 

GVPLOT (18) 

Store titles, number of parametric variations, total 
number of output coordinates, coordinate labels and 
final results for plotting. 


5 lOUNIT LK2RE3 (7) Read G400 rotor inipedance matrix and force vector as 

specified by the user. 


6 KPRINT (As Needed) Write titles, input data image, input for each 

component intermediate calculations for debugging and 
final results (standard write unit). 


7 KPUNCH (As Needed) Standard punch unit. 


8 


KRSTRT 


1. RECOVR (8) 

2. RECTUl (8) 

3. SV3TH1 (18) 


Read coordinate dimensions for time history restart cast:. 
Read data (displacements and derivatives) for restart case 
Write coordinate dimensions and data (displacements and 
derivatives) for restart case. 


10 IOXRE2 LK2RE2 (7) Read E927 stiffness, damping, and mass matrices. 

11 IOXRE3 LK2RE3 (7) • Read G400 rotor impedance matrix and force vector from 

program. 

12 IOXRE3 LK2RE3 (7) Read user input of rotor impedance matrix and force vector 

*Specificd in subroutine BLOCK 



TABLE 11 . NEW SUBROUTINES REQUIRED FOR ADDITION OF A SUBSTRUCTURE 


Number 

Routine 

Name 

Segment 

Number 

Called by Subroutine 
(in Segment Number) 

Description of New Routine 

1 

CINBF4 

1 

CONIKP (1) 

Converts input to internal units. 

2 

RECBF4 

1 

RECOVR (6) and 
RESTHl (17) 

Recovers displacements and velocities 
for restart cases. 

3 

PINBF4 

4 

PINPUT (4) 

Prints component input data. 

4 

CRDBF4 

9 

COORD (9) 

Forms coordi::=ite attributes. 

5 

RENBF4 

9 

RENTHl (9) 

Initializes displacements and velocities 
for restart cases. 

6 

NAMBF4 

10 

LABEL (10) 

Assigns * uixs to the. output coordinates. 

7 

CPBF4 

13 

CP (13) 

Foims component mass, damping and stiff- 
ness matrices. 


180 



TABLE 12. REVISIONS TO SUBROUTINES REQUIRED FOR ADDITION OF A SUBSTRUCTURE 


Number 

Subroutine 

Name 

Segment 

Number 


Description of 

Revisions 

1 

BLOCK 

. 1 

1. 

COMT'ION/CNAMES/ 1 Y ( 2 1 ) 





2. 

DIMENSION iX(21) 





3. 

DATA IX/4HFA1 , . . . , 

4HBF4 / 




4. 

DO 10 I = 1, 21 


2 

CONINP 

1 

1. 

COMMON/CNAMES/XFAl, . . 

. , XBF4 




2. 

INTEGER XBF4 





3. 

1F(ANAME.EQ.XBF4 ) CALL 

CINBF4(A(LOCDAT) ) 

3 

NSPACE 

3 

1. 

COMMON/CNAMES/ I X ( 2 1 ) 





2. 

DIMENSION IY(21) 



3. DATA IY/12, . . . , NDLOAD/ 

where NDLOAD = number of input items needed to 
specify new element BF4 

4. DO 10 I - 1, 21 

5. Define NWORK = number of work variables needed for 
BF4 ( 0 ) 

‘ 6. Revise NTOTAL = NDLOAD *■ NWORK as appropriate for 
BF4 


4 

PINPUT 

4 

1. 

2. 

3. 

COMWON/CNAMES/XFAl , . . . , XBF4 

INTEGER XBF4 

IF(ANAME.EQ.XBF4) CALL PINBF4 ( A( LOCDAT ) , I ELEM, KPRINT ) 

5 

LINK 

7 


Same as revisions 1 and 2 for PINPUT (seg. 

4) 

6 

RECOVR 

a 

1. 

2. 

3. 

Same as revisions 1 and 2 for PINPUT (seg. 4) 
IF{ANAME.EQ.XBF4) CALL RECBF4 (A(LOCDAT), 

NNG.A( IXDEP). A( IDXDEP) ) 

If no initial displacements and velocities are to 
be recovered for BF4, then code 
IF(ANAME.EQ.XBF4) NNC » NNG ♦ NCBF4 • 
where NCBF4 » number of coordinates associated 
with BF4.' 
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TABLE 12. CONCLUDED. 



Subroutine 

Segment 



Number 

Name 

Number 


Description of Revisions 

7 

COORD 

9 

1. 

Same as revisions 1 and 2 for PINPUT <seg. 4) 




2. 

IF(ANAME.EQ.XBF4) CALL CRDDF4 ( A(LOCDAT) , NL0CAL« 
NNODE'L, NCNL, EULERL ) 

8 

REN’nil 

9 

1. 

Same as revisions 1 and 2 for PINPUT (seg. 4) 




2. 

IF(ANAM£.EQ.XBF4) CALL RENBF4(A(LOCDAT) ,NNG,X0EP, 
DXDEP ) 




3. 

If no displacements and velocities are to be 
initialized for BF4, then code 





IF (ANAME.EQ.XDF4) NNG NNG * NCBF4 
(same as RECOVR above) 

9 

LABEL 

10 

1. 

Same as revisions 1 and 2 for PINPUT (seg. 4) 




2. 

IF(ANAME.EQ.XBF4) CALL NAriBF4 ( A( LNAME) , NLOCAL) 

10 

UPDATE 

12 


Same as revisions 1 and 2 for PINPUT (seg. 4) 

11 

ASSMBL 

13 


Same as revisions 1 and 2 for PINPUT (seg. 4) ' 

12 

CP 

13 

1. 

Same as revisions 1 and 2 for PINPUT (seg. 4) 




2. 

IF(ANAME.EQ.XBF4) CALL CPBF4 ( A( LOCDAT) , NLOGAL, MO, 
CO. KO, FO.IDIM) 

13 

SOLVE 

14 


Same asrevisions 1 and 2 for PINPUT (seg. 4) 

14 

SOIFRI 

16 


Same as above 

IS 

S02TH1 

16 


Same as above 

16 

RESTOR 

17 


Same as above 

17 

RESTHl 

17 


Same as revisions for RECOVR (seg. 8) 

16 

SAVE 

la 


Same as revisions 1 and 2 for PINPUT (seg. 4) 
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TABLE 13. BASE PROGRAM TEST CASES SUMMARY 





BASE PROGRAM 
SUBSTRUCTURES 



OPERATIONAL AND 
SOLUTION MODES 




Case 

No. 

Case Description 

BFl 

BFZ BHl FAl ISl RE2 RE3 

CN.l CFi 

MSI 

GEN 

PVl EGl EG2 

FRl 

TIU 

Size of 
work 
Vector 

CPU 

Time 

(Sec) 

i 

Inplnne bifilars with 5 fuselage 
modes using a forced 
renponse solution with longi- 
tudinal dfid lateral excitation 
forces . 

V 


V 

V 

V 


V 


2161 

o.s 

2 

Vertical bifilars with 5 fuselage 
nointal modes using a forced 
response solution with a vertical 
excitation force. 


V 

V 

V 

V 


V 


2307 

0.4 

3 

Elastic beam model comprising 
5 segments constrainted to 
vertical displacement using 
eigensolution 1. 


V 

V 


V 

V 



7865 

0.0 

4 

Fi.xed system absorber with mass 
added to base using a forced 
response solution with a 
longitudinal force at various 
ficqucncies. 


V 

V 

V 

V 

V 

V 


106S 

0.4 

*S 

As case 4 above using a time 
history solution calculated up 
to 1 second-line printer 
output suppressed. 


V 

V 

V 

V 

X 


V 

1032 

4.5 

6 

As case ^ above using the time 
history restart option from 0.05 
to 0.10 second-line printer output 
suppressed. 

V 

V 

V 

V 

X 


V 

1131 

0.5 

NOTES : 

1) V indicates substructures and solution method used. 

2) X indicates operational modes and solution methods which cannot be used. 

3) * input and output are provided in this manual. Plots are provided except 

for Case 12. 
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TABLE 13. CONTINUED. 






BASE PROGRAM 
SUBSTRUCTURES 



OPERATIONAL AND 
SOLUTION MODES 



Cifse 

No. 

Case Description 

BFl 

BF2 

BMl FAl 

ISl RE2 RE3 

CNl CFl 

MSI 

GEN PVl EGl EG2 

FRl 

Size of 
Work 

THl Vector 

CPU 

Time 

(Sec) 

*7 

One d.o.f. (Vertical Motion) 
nodal isolator witli 2 fuselage 
and one transmission normal mode 
using a forced response solution 
with a vertical excitation 
force at various frequen- 
cies (see Figure 44a). 




V 

V 

V 

V V 

V 

30S5 

1.0 

8 

Two d.o.f. (Vertical and Lateral 
Motion) nodal isolator with one 
transmission normal mode with 
the vertical force acting on 
node 1 (:see Figure 44b). 



• 

V 

V 

V 

V V 

V 

2794 

0.4 

9 

Two d.o.f. (Vertical and Lateral 
Motion) nodal fuselage and one 
transmisjion tiurmal mode with the 
vertical force acting on node 3 
(see Figure 44c). 




V 

V 

V 

V V 

V 

3090 

0.4 

10 

As Case 9 above but with 6 
normal modes for the fuse- 
lage and six normal modes for 
the transmission. 




V 

V 

V 

V V 

V 

5915 

0.7 

*11 

Six d.o.f. nodal isolator 
assembled from 4 isolator bars# 
positioned 90* to each other. 




V 

V 

V 

V V 

V 

+14831 

2.30 

NOTES ; 

1) V indicates substructures and solution method used. 

2) X indicates operational modes and solution methods which cannot be 

3) * input and output are provided in this manual. Plots are provided 

4) * case cannot be run in an interactive mode since size is greater 

used. 

except for Case 12. 
than 6400 words. 
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TABLE 13. CONCLUDED. 


BASE PROGRAM OPERATIONAL AND 

SUBSTRUCTURES SOLUTION MODES 


Size of CPU 

Case Work Time 

No. Case Description BFI BF2 BNl FAl ISl RE2 RE3 CNI GFl MSI GEN PVl EGl EG2 FRl THl Vector (Sec) 

*12 E927 elastic rotor with 12 rotor V V V V ♦14529 2.2 

coordinates and S fuselage inodes 
using elgensolution 2. 


*13 User input of rotor impedance V VVVVXXVX 1985 1.0 

matrix with 6 hub displacements 
and 4 fuselage normal modes using 
forced response solution. 

14 C400 program rotor impedance V.VVV XXVX 2190 0.6 

matrix with C hub displacements 

and 3 ARES model normal inodes. 

No IIHC and no stress options used. 

15 As Case 14 with the V VVV XXVX 5397 1.4 

stress option invoked. 

16 As Case 14 with both HHC V VVV XXVX 5538 1.5 

and stress options invoked. 

NOTES : 1) V indicates substructures and solution method used. 

2) X indicates operational modes and solution methods which cannot be used. 

3) * input and output are provided in this manual. Plots are provided except for Case 12. 

4) ♦ case cannot be run in an interactive mode since size is greater than '8400 words. 



TABLE 14. BASE PROGRAM TEST CASES TEKTRONIX PLOTS SUMMARY 


Test Figure Abscissa Ordinates 

Case No. Name Units Name Units Plot Description 


5 

45 

TIME 

Sec 

1. DELT 

Ft 





2. X 

II 


7 

46 

FREQUENCY Hz 

1. 2T 

G*b 




2. ZB 

II 


11 

47a 

FREQUENCY 

Hz 

1. XBF 

2. YBF 

3. 2BF 

G»s 

II 

II 


47b 

FREQUENCY 

Hz 

1-4 ZBF 

G*S 


47c 

FREQUENCY 

Hz 

1-4 

ZBF 

G' 

' s 

47d 

FREQUENCY 

Hz 

1-4 

ZBF 

G' 

' s 

47e 

FREQUENCY 

Hz 

1-4 

ZBF 

G' 

's 

47f 

FREQUENCY 

Hz 

1-4 

ZBF 

G' 

's 

47g 

FREQUENCY 

Hz 

1-4 

ZBF 

G' 

's 

47h 

FREQUENCY 

Hz 

1. 

ZBF 

G' 

's 


Time history response of fixed system absorber 
and modal structure longitudinal displacements. 
Response increases with time because system is 
not restrained . 

Forced response of transmission and fuselage 
vertical displacements with one 1-0 isolator 
system for a forcing frequency range up to 
50 hz. Fuselage motion is a minimum near fre- 
quencies of 17 and 27 hz and a maximum near 21 hz. 

Forced response of fuselage displacements at the 
location of the first 2-D isolator for a system 
made up by four 2-D isolators. Frequency range 
is up to SO hz . 

Forced response of fuselage vertical motion at 
the four isolator locations. Positions 1 and 3 
show almost the same response; similarly for 
positions 2 and 4. 

Figure 45c is a perspective plot of the results 
of Figure 45b with zero density value. 

Same as Figure 45c with the default density value. 
Same as Figure 45d with depth reversed. 

Same as Figure 4Sd with abscissa reversed. 

Same as Figure 45d with ordinate reversed. 

Forced response of fuselage vertical motion at 
the location of the first 2-D isolator- as a function 
of isolator mass and frequency. Test case 11 was 
run for 5 values of isolator mass: 24 to 64 pounds 
in increments of 10 pounds. The response is plotted 
for cases 1, 3 and 5 to show the plot package capa- 
bility of plotting data from different cases (up 
to 5 cases can be used). 
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TABLE 14. CONCLUDED. 


Test 

Fi gure 

Abscissa 


Ordinates 



Case 

No. 

Name 

Units 

Name 

Units 

Plot Description 


13 

48a 

3RD MODE 

Hz 

1. 

. RX 

Lbs 

Forced response results of hub forces 

using 



FREQ 


2. 

. RY 

•' 

user input of rotor impedance matrix < 

and force 





3. 

. RZ 

II 

vector for a frequency up to 30 hz. 



48b 

3RD NODE 

Hz 

1. 

. X 

G*s 

As Figure 46a but with display of hub 

accelera-> 



FREQ 


2. 

, y 

II 

tions. 






3. 

. z 

II 




187 



APPENDIX A: G400/F389 COUPLED PROGRAM 

A.l PROGRAM STRUCTURE AND USAGE 

The computer programs G400 and F389 were coupled by means of a 
Fortran driver program to allow blade circulations and rotor 
induced velocity fields to be calculated iteratively to achieve a 
consistent rotor aeroelastic response and variable inflow state 
among the programs. The G400 program' analysis and user's manual 
are presented in References 2 and 3. . The F389 program analysis 
and user's manual are described in References 8. 

The driver program is called VIBIT. The coupled program seg- 
mentation structure is illustrated in Figure A.l. From this 
figure, it is seen that the G400 program is made up of 3 segments: 
G400PG, NIAM and MOTION, while the F389 program consists of 6 
segments: F389PG, SECONE, SECTWO, SECTRE, SECFOR and SECFIV. A 

description of the 10 segments making up the G400/F389 program 
appears in Table A.l. CDC dependent routines are: RESETQ (seg- 

ment 2), NIAM (segment 3), SREC and SKPFIL (both in segment 9). 

The logic of the coupled program allows 3 modes of operation: 
G400 alone, F389 alone, and G400/F389 coupled operation. These 
options are illustrated in the logic flow sketch from Figure A. 2. 
The control parameter, N, is the first item of input loaded in the 
program. It controls the different modes of operation as shown in 
the table attached to Figure A. 2. 

The data transfer files required for the execution of the G400 and 
F389 programs in both the uncoupled and coupled modes are shown in 
Figure A. 3. A summary and a description of the various input/ 
output data files is presented in Table A. 2. 

A sample input runstream is presented in Table A. 3 for a 2-cycle 
operation. To run G400 alone, data blocks 1 and 2 only are 
needed. To run F389 alone, data blocks 1 and 3 are needed; 
however, in this case, the user must also provide the linearized 
airfoil characteristics and the axial velocities (if desired) as 
part of input block 3. To rune one complete cycle, then data 
blocks 1 through 4 are needed. Four input locations for thie G400 
data block 2 must be given the values listed below: 

location 48 = 1.0 - calculate harmonics of blade responses. 

location 76 = 1.0 - write initial conditions to unit 24. 

location 998 - DELTA PSI - azimuthal increment for F389 data. 

location 999 = 1.0 - transfer of data to F389 through units 

13, 16 and 23. 
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The F389 input locations 88 and 206 in data block 3 must be set to 
1.0 for coupling with the G400 program. Data block 4 contains 
only new or revised input locations as compared to the G400 data 
already provided in data block 2 which were stored in unit 24 for 
later use. Examples of revised input parameters are: tip loss 
factor should be set to 1.0 (location 7), variable inflow option 
is activated (locations 53 and 54 are set to 1.0), new initial 
conditions should be read from unit 26 (location 76 equals 2.0). 
Other parameters may be changed as desired by the user. To 
execute more than one cycle, the input data blocks 5 and 6 are 
repeated. The F389 data in block 5 contain only new or revised 
locations because the previous loader data from block 3 have been 
stored in unit 25 for use in subsequent F389 runs. The hew 
character of block 5 data for F389 is analogous to the block 4 
data for G400. A precautionary observation is that new or revised 
input data provided in data blocks 4 and 5 are not written to 
units 24 and 25 respectively. Thus, when more than one G400/F389 
cycle is requested, the new or revised input parameters must be 
relocated. For example, in Table A. 3, G400 data block 6 must 
duplicate the revisions provided in data block 4 if they are. still 
required. 

A. 2 PROGRAM EXECUTION 

The execution of the G400/F389 coupled program for one cycle is 
presented in Figure A. 4 and Appendices H and I. Figure A. 4 shows 
the Job Control Language (JCL) used to run the program on the IBM 
370/68 computer system. The input data are contained in the file 
specified under unit 5. A listing of the complete input data is 
provided in Appendix H. Referring to Table A. 3, the input data 
contain the cycle control option (data block 1), G400 data for the 
first pass (data block 2), F389 data for the first pass (data 
block 3 ) , and G400 data for the second and final ( for this exam- 
ple) pass (data block 4). The output for the coupled run is shown 
in Appendix I. Additional results can be printed out for both 
programs but are not shown for this sample run to minimize the 
size of this report. 

A. 3 CPC PROGRAM INSTALLATION AND OPERATION 

This section provides the Job Control Language (JCL) for the 
installation and operation of the G400/F389 program on the NASA- 
Langley CDC NOS computer facility. The installation and operation 
procedures are based on protocol and system library facilities 
associated with the NASA-Langley CDC NOS computer facility. 

The JCL for installation of the G400/F389 program is presented in 
Figure A5 . 

The JCL for execution of the G400/F389 program is presented in 
Figure A6 . 

The JCL instructions apply to remote job entry operation (such as 
a telephone linkup). Revisions needed to run with card images are 
indicated in Figures A5 and A6. 


189 







r 




F389 INTERNAL STORAGE /WORK FILES 


Figure A3. G400/F389 Input/Output Files. 
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START 

COL 


- — ♦— -- 1-- — — —4 ♦ 


2 . 

3 


I //ET473RUN JOB (4045*STAL»15»15)»*CASSARXNO - X735d* 
1 // CLASSsO>TIHEs0l»RSGIONs3000K> 

1 // ttSGCUSSsT.HSSLEVEL=l,N0TIFYsET473 

1 
1 
1 

“1 
1 
1 
1 


//• 

//* 

■7/* 

//** 
//STCPOX 
//STEPLIB 


G400/F3d9 COUPLED PROGRAM 
JCL FROM •BT473.G400F3d9.CMTL(REPORT)' 
"E5?KcOTEiTNnEYCrErS) - 6-400/F-3S9/G-460 


EXEC 

OD 


PGMsCOUPLEOX,TIME»10 

DSNsET473.6<i00F3e9«LOAD,DlSPsSHR 


0001 

0002 

0003 

0004 

0005 

0006 
0007 

'oooe 

0009 

0010 
_oon 
>0012 

0013 

0014 

0015 

’ m6 

0017 

0018 
>0019 

-0021 

0022 

-0023 

'b024 

-0025 

0026 

-0027 

-0029 

0030 

-0031 

mu 

-0033 

0034 

-0035 

0037 

-0038 

0039 

-0D4ff 

0041 

0042 
-0043 


1 

1 

1 

■ I" 
1 
1 
I 


//* 

//» 1 . 
//* 2. 
77*"“3r 
//* 4. 
//* 5. 
//*— — 


INPUT DATA SET IN UNIT 5 (FT05F001) CONTAINS: 

===> RUN OPTION INPUT CARD - 12 FORMAT 
srs> S-400 INPUT - AIRFOIL* BLADE DATA* AND MODE SHAPES 
-sss5 F=1OT INPUT “LOXDTJNCY DITA">40T PISSED THRO BY 
S5s> 6-400 INPUT - NEU LOADER DATA IF HEEDED 
S5S> REPEAT CARDS 3 A 4 ABOVE FOR F-389 I G-400 NEU DATA 


1 //* STAND AR D " ' REH TU HI T 

4- 1 //FT05F001 00 0SN=ET473.G400F389.DATAfCTCLESl)*DISP=SHR 

1 //*» 

1 “7/R STANDARD WRITE UNIT 

1 //* 

r 1 //FTO6FOOI DO SYSOUT«T 
1 //»- 


77» STA ND A RrT 0N C H ~ UNI 1 

//FT07F001 DO SYSOUT=T 

77«W^*»»^^“6=50trPRDGff»^"aSE5 FILES IT*Z?*ZS 

//* 

//* UNIT 11- G-400 DATA ONLY - IMPEDANCE FOR RDYNE 


FORMATTED 


// FT IIF O OI ' 0 D ‘ D SNg«HLEll*UNI T^ feVH DH R*5 P AClsr T HK " *(2 07nnT; 

// 0CBs(RECFMsFB>LRECLs80*BLKSIZEsl9040)»DZSPs(NEU«DELETE) 

//m 

//« — 


1 

1 

I 

-X 

1 

1 

1 


UNIT 24- 6-400 DATA ONLY - LOADER DATA 


FORMATTED 


//FT24F001 DO DSNsAAFILE24tUNITsVUQRK>SPACEs(TRK. (20*10) )* 

// DCB3(RECFM«FB»LRECL«120.8LXSIZEsl8960),DXSPs(NEUiOELETE) 

//*-- 


^ 1 //FT26F001 DD DSNsUFILE26*UNZTxVUORK.SPACEK(TRK>(20>10n. 0046 

I // DCBx(RECFMsFB»LRECL:a0»BLKSIZExl9040)»DXSPs(NEU. DELETE) 0047 

1 //iHHHHHHf f-389 program USES FILES 2*3*15»25*28*29*30 inMni w mm ii M i i i ni »ii nf 0049 

1 //* UNIT 15- F-389 DATA ONLY - UABAT HARMONICS - FORMATTED 0051 

1 : D0 5 2 

\P 1 //FT15F001 DD 0SN=ET473.F389.OATA(UABAT2).OI5PxSHR 0053 


Figure A4. G400/F389 IBM Job Control Language 



1 //* LfNIT 2S- F*389 LOADER DATA - FORMATTED 0055 

1 //• 0056 

M I //FT25F001 DO DSNs*iFILE25.UNITsVMORK*SPACEsrn?K*( 20 »10 H » 0057 

X~^? 0CBsrR£CFttiFBTlRECtsl2o;BUCSt2£5l696bV;blSPs(N " OOSS 

I //»•••.— 0059 

1 //*•—.. f-389 program USES FILES 2,3»2a*29»30 FOR TEMPORARY STORAGE 0060 

1 — -0061 

‘ n* ~X “7/PT02P0M fib bSNnAPlLCB2*lNitsVKORK»SPACEsftRK;(20 0062 

1 // 0CBs(RECFM:VeS.LRECL8l9060»DLKSIZEsl9069). 0063 

1 // 01SP3(NtU»DELETE) 006A 

l> 1 //FT03F001 00 0SN:|lFILE03.UHIT:VUORK»SPAC£s(TRK»(20>10))» 0065 

1 — 7 / TCe=<mrh'^BsnjiECLmb 6 i>;BLK§m=^^ 0066 

1 // OISPs(NEW,OELETE) 0067 

14' I //FT28F001 DO DSNsUFlLE28>UNlT:VUORK»SPACE:(TRK> (20 »10 ) ) » 0068 

1 // DC6=(RECFM=V8S.LRECLsl9060»BLKSIZE-19069)» 0069 

I’ 7/ OISP=rWEH.OELETEl 0070 

IS" 1 //FT29FC01 00 OSN=ttFILE29,UNITsVMORK,SPACEs(TRK» ( 20»10 ) ) » 0071 

1 // OCB=(RECFM:VBS>LRECL=19060»BLKSIZEsl9069)» 0072 

1 // OISP=(NEW*OELETEI 0073 

" ’ TCI "' //F » AOKU OI~DO~P^mnLr307UMnsWDRKTSPACE=(TRKTT2b;iOin 0079 

1 // OCB5(RECFM=VBS.LRECL=19060,BLKSIZEs19069), 0075 

X // OISPsCNEWtOELETEI 0076 

1 //» -———^•....,-0077 

• — “I “7/ir--=-=- G-AOO/F=3B9 COUPUITPROGRAM^ALSO OSES ■FItES“13‘;i4Tl6723” ■“ "0078 

1 //* — .—.—...0079 

1 //* — . G-AOO OATA FOR THE F-389 PROGRAM - FILES I3»16>23 0060 

1 //*.— — ...... .........—0081 

1 //u - W IT 13- G- gO(y"PATA' FOgT-389~':~TO’^5T^^^ FORMAnED'' — 0082 

1 //» 0083 

n 1 //FT13F001 00 DSN=UFILE13#UNITsVWORK,SPACEs(TRK,(20»10l)t 008A 

1 // DCB:( RECFMsFB* LRECL«80 >BLXSIZE:190A0 ) »DISPs( NEU »DELETE ) 0085 

1~ 77M-~-:r=-— -.-0ofe6 

1 UNIT 16- 6-400 OATA FOR F-389 - BUOE INPUT - FORMATTEO 0087 

1 //* — ........ 0088 

IP 1 //FT16F001 DO DSNsttFILE16»UNITsW!ORK>SPACE:(TRK,(20»10))t 0089 

r~V7 TJCBSTRECFM=FB-,"tRECt2I20-,"BlKSIZr='I6960J'; 0090 

1 // DISP=(NEM»DELETE) 0091 

1 //*- ..... ... — .............. 0092 

1 //* UNIT 23- G-400 OATA FOR F-389 - AIRFOIL DATA - FORMATTEO 0093 

In 1 //FT23F001 DO OSNsAAFlLE23»UNITsvuORK»SPACEs(TRK»(20>10n> 0095 

1 // DCBs(RECFMsFB>LRECLsl20»BLKSIZEsl8960)> 0096 

1 // OISP=(NEW.OELETE) 0097 

X ' ■ 08 98 

1 //* F-389 DATA FOR THE G-400 - FILE 14 0099 

1 //*— ...... 0100 

X //» UNIT 14- F-389 OATA FOR 6-400 - HARMONICS V.I.- FORMATTED 0101 

10 1 //FT14F001 DO DSN=iAFILE14*UNlT=VM0RK,SPACE=( TRK, ( 20,10 )) . 0103 

1 // 0CB:(RECFM=F6.LRECL:d0»BLKSIZEsl9040)» 0104 

1 // DISP=(NEM.DELETE1 ^ 0105 

1 /» <==s EHD OF JCL CONTROL CARDS 0107 

1 //* 0108 

11 1 // 0109 


Figure A4. Concluded. 
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/JOB 

* /NOSEQ 

VIBITL,T1400,CM300000. 

USER , . 

CHARGE, ,LRC. 

RFL(300000> 

GET.FILE3. 

GET, FILEA. 

GET,FILE5X. 

GET, FILES. 

GET,FILE7. 

FTH(I=FILE3,R=0,L=TAPE10, B=VIBGO, A) 
FTNtI=FILEA,R=0,L=TAPEll,B=VIBlGO,A) 

FTN( I=FILE5X, R=0,L=TAPE12 ,B=VIB2G0, A) 

FTN( I=FILE6,R=0,L=TAPE13, B=VIB3G0, A) 

FTH( I=FILE7,R=0,L=TAPE1A, B=VIBAGO, A) 

REPLACE ( VIBGO) 

REPLACEC VIBIGO) 

REPLACE ( VIB2G0) 

REPLACE(VIB3G0) 

REPLACE( VIB4G0) 

ATTACHC ALTMLIB/UN=LIBRARY) 

SEGLOAD(B=VIBABS) 

LDSET( LI B=ALTMLIB, PRESET =2ER0,MAP=SBEX/LMAP ) 
LOAD( VIBGO) 

LOAD(VIBIGO) 

LOAD(VIB2GO) 

LOAD( VIB3G0) 

LOAD( VIBAGO) 

NOGO. 

REPLACE(VIBABS) 

DAYFILE-, JCLOUT. 

REPLACE, JCLOUT. 

REPLACE(TAPE10=FORT3) 

REPLACE (TAPEll=FORT A) 

REP LACE < TAP E12 = FORT5X) 

REPLACE(TAPE13=F0RT6> 

REPLACE(TAPE1A=F0RT7) 

REPLACE , LMAP . 

EXIT . 

DAYFILE, JCLOUT. 

REPLACE, JCLOUT. 

REPLACE ( TAPEl 0=FORT3) 

REPLACE(TAPEll=FORTA) 

REPLACE(TAPE12=F0RT5X) 

REPLACE(TAPE13=F0RT6) 

REPLACE(TAPE1A=F0RT7) 

REPLACE, LMAP. 

/EOR 


FIGURE A5. CDC JCL FOR INSTALLATION OF THE GAOO/F3S9 PROGRAM. 



F3S9PG GLOBAL BET DAT , BLKDAT , BLOCKl , CONST 1 , C0NST2 , C0NST3 , CONTR L-SAVE 
F3S9PG GLOBAL CWAKES , GR P L, LAFDAT , P S I VAR , STOP I T , TEMP ST , XAFDAT-S AVE 
GAOOPG GLOBAL ADER I V , B LDDOF , BL OCK , BLSTR S , CBR V AR , CLDMTB , DAT AL V-S AVE 
GAOOPG GLOBAL DATEET , DOFVAR , DSINTS , FLXBTC, GLAURT , HARR AY , HRMINF-SAVE 
GAO OP G GLOBAL IMISCl , IMP DAT , I NRTHB , I OF 3S9 , L I NE AR , MI SCI , MI SC2-S AVE 
G4 0 0PG GLOBAL MI SC3 , MI SCA.MI SC5 , NVDAMP , P NDULM , P HDVAR , S I NT S-S AVE 
GAOOPG GLOBAL SP NVAR , TR IM, TRHSHT , TTBNDG-SAVE 
SUB TREE F3S9PG-(SECONE,SECTWO,SECTRE-(SECFOR ,SECFIV> ) 

TREE VIBIT-(GAOOPG-<NIAM,MOTION> ,SUB) 

END VIBIT 

/EOF 

* For card images operation: 

1. delete /JOB and /NOSEQ instructions (-first 2 cards). 

2. replace /EOR with a 6/7/S multiple punch card. 

3. replace /EOF with a 6/7/S/9 multiple punch card. 


FIGURE A5. CONCLUDED. 
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it /JOB 
it /MOSEQ 

VIBRUN, T1 400, CM30 0000 . 

USER,— T- . 

CHARGE, ,LRC. 

RFL(300000) 

GET,VIBABS. 

GET,TAPE5=G400R3. 

GET,TAPE15=F3 9WAB2 <==== HARMONICS OF FUSELAGE INDUCED 

VIBABS, TAPES. VELOCITIES FROM WABAT PROGRAM 

DAYFILE.JCLOUT. 

REPLACE, JCLOUT. 

R EP LACE ( OUTPUT: OUT VI B) 

REPLACE(TAPE'1.1 = G400M3) < = = = = I.NPEDANCE MATRICES STORED IN UNIT 11 

EXIT. 

DAYFILE.JCLOUT. 

REPLACE, JCLOUT. 

R EP L ACE ( OUTPUT: OUT VI B) 

REPLACE(TAPE11:G400M3) 

DMD(200) 
it /EOF 

it For card images operation; 

1. delete /JOB and /NOSEQ instructions (-first 2 cards). 

2. replace /EOF with a 6/7/S/9 multiple punch card. 


FIGURE A6. CDC JCL FOR EXECUTION OF THE GA00/F3S9 PROGRAM. 
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TABLE A.l. G400/F389 SEGMENTATION STRUCTURE DESCRIPTION 


Segment 

Leading 

Segment Description 

Number 

Routine 

1 

VIBIT 

Controls execution of the uncoupled and' coupled program versions. 

2 

G400PG 

Controls execution of the G400' program. Reads and writes airfoil data and 
initializes variables. 

3 

NIAM 

Reads and prints input data for blade characteristics mode shapes, and vari- 
able inflow, forms various papameters, sets up deflection arrays, and calcul- 
ates pitch coupling terms. 

4 

. MOTION 

Initializes variables, sets up data for hub motion and pendulum absorbers, 
calculates trim required quantities, solves for the time history response, 
calculates rotor performance, blade stresses, and hub loads, and performs 

• ■ 


harmonic analysis of final results. 

5 

F389PG 

Controls execution of the F389 program. Initializes the loader vector. 

6 

SECONE 

Reads and prints input loader data. 

7 

SECTWO 

Reads airfoil data, calculates rotor attitude, blade flapping, and blade coordi- 
nates, sets up wake model, calculates aerodynamic parameters, reads axial 
velocities, and prints out intermediate results. 

8 

SECTRE 

Controls input of the wake geometry and solution of the rotor induced velocity 
field. 

9 

SECFOR 

Sets up the rotor wake geometry. 

10 

SECFIV 

Solves for the rotor induced velocity field and provides output of the cir- 
culations and the variable inflow velocities and harmonics. 
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TABLE A. 2. G400/F389 I NPUT/OUTPUT FILES 


Fila 

No. 

1 

2 

3 

5 

6 
7 

11 

13 

14 

15 

16 

23 

24 

25 

26 


F389 G400 G400/F389 

Mode Mode Mode Description of Input/Ontput Files 


X 

X 

X 

X X 

X X 


1. G400 output of initial conditions if loc 76 = 0 

2. G400 output of blade and airfoil data and axial velocities 
for F389 

X Temporary storage of circulation solution geometric influence 

coefficients in block form 

X Temporary storage of intermediate solution quantities 

X Standard read file 

X Standard write file 


X 

X X 

X 
X 
X 

X 
X . 

X 


standard punch file - F389 output of variable inflow 

G400 save unit for rotor impedance matrix and force vector 

G400 save unit for axial velocities transmitted to F389 program 

F389 save unit for variable inflow transmitted to 6400 program 

F389 input of induced velocities from the WABAT program 

G400 save unit for blade data transmitted to F389 progreim 

G400 save unit for linearized airfoil data transmitted to F389 
program 

G400 loader data storage for coupled operation 


X F389 loader data storage for coupled operation 


X X 


G400 output of initial conditions if loc 76 = 1 or 2. 
Initial conditions are read from unit 26 only if loc 76 = 2 
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TABLE A. 2. CONCLUDED. 


File 

No. 

F389 

Mode 

G400 

Mode 

G400/F389 

Mode 

Description of Input/Output Files 

28 

X 


X 

Temporary storage of circulation solution geometric influence 
coefficients 

29 

X 


X 

Temporary storage of far field geometric influence coefficients 

30 

X 


X 

Storage of wake coordinates 



TABLE A. 3. G400/F389 INPUT OUTLINE FOR TWO CYCLES 


Data 

Block 

Program 
G400 or F389 

Input Description 

1 

G400 and F389 

Load N in 12 format for all runs: 

N = -2 . - run F389 alone 

N = -1,0. - run G400 alone 

N 1. - run N cycles of G400/F389 

2 

G400 

G400 program input data for: 

a) airfoil characteristics 

b) blade data 

c) mode shapes 

d) variable inflow (if desired) 
Data stored in unit 24. 

3 

F389 

F389 program input data. Data stored in unit 25. 

4 

G400 

For the 2nd pass on the G400 program, only new or 
revised locations need to be loaded. The G400 data 
of block 2 above are stored in unit 24 and read in 
subsequent calls to G400. 

5 

F389 

If more than one cycle is requested, then only new 
or revised locations need to be loaded for the 2nd 
pass on the F389 program. The F389 data of block 3 
above are stored in unit 25. 

6 

G400 

This is the 3rd pass on the G400 program for 2 or 
more cycles. These data could be the same as that c 
block 4 above. Repeat data blocks 5 and 6 for more 
. cycles . 


O 



APPENDIX B: E927 PROGRAM 


B.l PROGRAM EXECUTION 

A brief description of the Rotor Aeroelastic Stability Analysis 
(Program E927) is provided in section 5.2 of this report. The 
instructions for executing the E927 program are provided in 
Reference 5. For this contract, location 110, which controls the 
coupling of the rotor matrices with the bifilar analysis, is not 
applicable. The input value is ignored in the revised program 
logic. 

The IBM job control language (JCL) needed to execute the E927 
program is presented in Figure Bl. Unit 10 is used to store the 
rotor matrices (mass, damping, and stiffness), as indicated by the 
control card no. 8. 

A typical input case is shown in Appendix J. The generated output 
is presented in Appendix K. The rotor matrices printed out at the 
end of the case are stored in unit 10 and are used subsequently in 
the execution of the Base Program (SIMVIB). Case 12 from Table 13 
employs the rotor matrices from this sample run. ’ 

B.2 CPC PROGRAM INSTALLATION AND OPERATION 

This section provides the Job Control Language (JCL) . for the 
installation and operation of the G400/F389 program on the NASA- 
Langley CPC NOS computer facility. The installation and operation 
procedures are based on protocol and system library facilities 
associated with the NASA-Langley CPC NOS computer facility. 

The JCL for installation of the E927 program is presented in 
Figure B2 . 

The JCL for execution of the E927 program is presented in Figure 
B3. 

The JCL instructions apply to remote job entry operation (such as 
a telephone linkup). Revisions needed to run with card images are 
indicated in Figures B2 and B3 . 
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//ET473RUr< JOB ( 4045,STAL,19,153 , ‘CASSARINO - X7358' , 

// CLASS=0 ,TIME=01 ,MSGLEVEL=( 1,0) ,REGIOM=3000K, 

// MSGCLASS=T,N0TIFY=ET<+73 

EXECUTIOM OF E-927 PROGRAM REVISED FOR SIMVIB PROGRAM 

//STEPOl EXEC PGM=E90BCSV1 

//STEPLIB OD DSM=ET473.E90BCHAS.LOAD,DISP=SHR 

-------- — _--—---—- — — — — — - — 

»** READ IHPUT DATA FROM UNIT 5 



//FT05F001 DD 0SN=ET473. SIMVIB. DATAC E927RUN ) ,DISP=SHR 

»,(*-_ -; . 

*** k'RITE OUTPUT TO UNIT 6 

/y* 

//FT06F001 DD SYSCUT=T 

— , — — — 

*** PROGRAM USES WORK FILES 8 AND 11 

yyy 

//FT08F001 DD UNIT=NORK,DISP=(NEW,DELETE ) ,SPACE=(TRK, ( 5,1 ) ) , 

// _.DCB=(RECFM=VSB,LRECL=680,BLKSIZE=2050) 

//FTllFOOl DD UNIT=WORK,DISP=(NEW, DELETE), SPACE=(TRK, (13,1)), 

// DCB=(RECFM=VSB,LRECL=40<f,BLKSIZE=2050) 

yyy — 

*»» OUTPUT MATRICES ARE. KHITTEM TO. UNIT. . 10 FOR SIMVIB PROGRAM.' 

yii«. — 

//FTlOFOOl DD DSN=ET473.E927CUT.DATA,DISP=(NEM,CATLG), 


// UNIT=TS0,SPACE=(TRK,(5,5)), 

// DCB=(RECFM=FB,LRECL=30,BLKSIZE = 3120) 

; 

END OF JCL SETUP 

yyy 




Figure Bl. E927 Program IBM JCL 



/JOB 

=< /NOSEQ 

ES27LHK, T1 400, CM30 0000. 

USER , . 

CHARGE, ,LRC. 

RFLCSOOOOO) 

GET,E927BF. 

GET,E927SV. 

FTN<I = E927BF, ■? = 0 , L = T AP E 1 , B=E9BFG0 , A) 

FTM( I=E927SV,R=0,L=TAPE2,B=E9SVGO, A) 

REPLACE(E9BFGO) 

REPLACE (E9SVGO) 

ATTACH( ALTMLIB/UN=LIBRARY) 

SEGLOAD<B=ESVABS) 

LDSET(LIB=ALTMLIB,PRESET=ZERO,MAP=SBEX/LMAP) 

LOAD(E9BFGO) 

LOAD(E9SVGO) 

HOGO. 

REPLACE(ESVABS) 

DAYFILE, JCLOUT. 

REPLACE, JCLOUT, 

REPLACE(TAPE1=F0RTBF) 

REPLACE(TAPE2=FORTSV) 

REPLACE, LMAP . 

EXIT. 

DAYFILE, JCLOUT. 

REPLACE, JCLOUT. 

REPLACE(TAPE1=F0RTBF) 

REPLACE(TAPE2=F0RTSV) 

REP LACE, LMAP. 

/EOR 

SHAKIT GLOBAL DYH I NP , 1 HE I GN , 1 HE I G , EOF 6 , P R NT SW , TMDS , N IMI C , L AGD AM 
TREE SHAKIT-(PRELIM,DYMMAT,AERMAT, EIGER) 

END SHAKIT 

^ /EOF 

* For card images operation: 

1. delete /JOB and /NOSEQ instructions (-first 2 cards). 

2. replace /EOR with a 6/7/S multiple punch card. 

3. replace /EOF with a 6/7/S/9 multiple punch card. 


FIGURE B2. CDC JCL FOR INSTALLATION OF THE E927 PROGRAM 



^ ./Job 

-> /NOSEQ 

ESa/RUN.Tl'+OO.CMSOOOOO. 

USER, . 

CHARGE, ,LRC. 

RFL(SOOOOO) 

GET,ESVABS. 

GET,TAPE5=ES27R1. 

ESVABS, TAPES. 

DAYFILE, JCLOUT. 

REPLACE, JCLOUT. 

REP LACE ( TAP E6 = E9270UT) 

REPLACE(TAPE10=E927MAT> <=== ROTOR & HUB MATRICES STORED IH UNIT 10 
EXIT. 

DAYFILE, JCLOUT. 

REPLACE,JCLOUTi' 

REPLACE(TAPE6=E9270UT) 

REPLACE(TAPE10=E927MAT) 

DMD(200) 

* /EOF 

^ For card images operation: 

1. delete /JOB and /NOSEQ instructions (-first 2 cards). 

2. replace /EOF with a. 6/7/S/9 multiple punch card. 


FIGURE B3. CDC JCL FOR EXECUTION OF THE E927 PROGRAM 
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